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A Reaction Motor 


2,000 YEARS AGO... 


AND NOW 


@ Hero was a city boy from Alexandria (Egypt) who was 
always pyramid-lengths ahead of his class when it came to 
things mechanical. Tinkering around one day in the year 
120 B.C. he came up with a gadget that would revolve 
indefinitely as long as he kept a brisk fire under it. Called 
the “Aeolipile,” this device was a primitive reaction motor. 
But alas, Hero’s engine ended up as a mere curiosity — 
for lack of an application. 


Today, challenging new requirements are constantly creat- 
ing the need for power sources which only rocket engines— 
modern reaction motors—can satisfy. RMI engineers have 
successfully harnessed this reaction principle for a wide 
variety of important airframe propulsion requirements 
and are extending the application of rocket power to other 
equally important fields. 


4 Static testing a new RMI rocket 
engine of the most advanced design. 


Career opportunities available for experienced mechanical, 
aeronautical, electrical and chemical engineers, physicists, 
chemists. Send complete resume to employment manager. 


REACTION MOTORS. INC. 


Denville, New Jersey 
Affiliated with OLIN MATHIESON CHEMICAL CORP. 
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Scope of JET PROPULSION 


Jet Proputsion, the Journal of the American Rocket Society, is 
devoted to the advancement of the field of jet propulsion through the 
publication of original papers disclosing new knowledge and new de- 
yelopments. The term “‘jet propulsion”’ as used herein is understood 
to embrace all engines that develop thrust by rearward discharge of a 
jet through a nozzle or duct; and thus it includes systems utilizing 
atmospheric air and underwater systems, as well as rocket engines. 
Jet PROPULSION is open to contributions, either fundamental or ap- 
plied, dealing with specialized aspects of jet and rocket propulsion, 
such as fuels and propellants, combustion, heat transfer, high tem- 
perature materials, mechanical design analyses, flight mechanics of 
jet-propelled vehicles, astronautics, and so forth. Jet PropuLsion 
endeavors, also, to keep its subscribers informed of the affairs of the 
Society and of outstanding events in the rocket and jet propulsion 


Limitation of Responsibility 


Statements and opinions expressed in JET PROPULSION are to be 
understood as the individual expressions of the authors and do not 
necessarily reflect the views of the Editors or the Society. 


Subscription Rates 


One year (twelve monthly issues).................02.00-- $12.50 

Foreign countries, additional postage............... add .50 


Change of Address 


Notices of change of address should be sent to the Secretary of the 
Society at least 30 days prior to the date of publication. 


Information for Authors 


Preparation of Manuscripts 


Manuscripts must be double spaced on one side of paper only with 
wide margins to allow for instructions to printer. Submit two copies: 
original and first carbon. Include a 100-200 word abstract of paper. 
The title of the paper should be brief to simplify indexing. The 
author’s name should be given without title, degree, or honor. A 
footnote on the first page should indicate the author’s position and 
affiliation. Include only essential illustrations, tables, and mathe- 
matics. References should be grouped at the end of the manuscript; 
footnotes are reserved for comments on the text. Use American 
Standard symbols and abbreviations published by the American 
Standards Association. Greek letters should be identified clearly for 
the printer. References should be given as follows: For Journal 
Articles: Authors, Title, Journal, Volume, Year, Page Numbers. 
For Books: Author, Title, Publisher, City, Edition, Year, Page 
Numbers. Line drawings must be made with India ink on white 
paper or tracing cloth. Lettering on drawings should be large 
enough to permit reduction to standard one-column width, except 
for unusually complex drawings where such reduction would be pro- 
hibitive. Photographs should be clear, glossy prints. Legends must 
accompany each illustration submitted and should be listed in order 
on a separate sheet of paper. 


Security Clearance 

Manuscripts must be accompanied by written assurance as to 
security clearance in the event the subject matter of the manuscript 
is considered to lie in a classified area. Alternatively, written assur- 
ance that clearance is unnecessary should be submitted. Full respon- 
sibility for obtaining authoritative clearance rests with the author. 


Submission of Manuscripts 

Manuscripts should be submitted in duplicate to the Editor-in- 
Chief, Martin Summerfield, Professor of Aeronautical Engineering, 
Princeton University, Princeton, N. J 


Manuscripts Presented at ARS Meetings 


A manuscript submitted to the ARS Program Chairman and 
accepted for presentation at a national meeting will automatically 
be referred to the Editors for consideration for publication in Jer 
Propusion, unless a contrary request is made by the author. 


To Order Reprints 


Prices for reprints will be sent to the author with the galley proof, 
and orders should accompany the corrected galley when it is returned 
to the Assistant Editor. 
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FORD INSTRUMENT 
QUANTIZES 
LIGHT 


Colored lines or spots on a piece of paper can become a means 
of conveying information — Rorschach charts, impressionist 
paintings and survey maps are all visual message carriers. Ford 
Instrument engineers found it necessary to translate such color 
information into electrical or mechanical quantities, (quantiza- 
tion) with less distortion than is inherent in the usual photo- 
graphic techniques. Such quantities can in turn be used as signals 
that actuate computers, make offset plates, and generally put to 
use the information implied by the difference between the colors 
or the distribution of the colors. 


REO ORANGE YELLOW GREEN BLU 
QUANTIZED OUTPUT 
DISTINCT COLOR SEPARATION OR 


The quantization performed by Ford is not restricted to color alone. 
For example, a black and white photograph represents an aggre- 
gate of light and dark areas of varying shades, and this display 
must frequently be converted into continuous or discrete electrical 
quantities for various purposes and uses. Ford engineers recently 
developed equipment which can quantize and record the various 
degrees of color, or gray areas in photographic negatives, and to 
correlate this information into usable data. This equipment was 
developed for a classified project — the equipment is unavailable 
for general use — however the technical know-how gained by 
Ford — combined with Ford’s superior pruduction and engineer- 
ing facilities —is available in the creation of light quantizing 
equipment for you. 

Light quantizing is but one of the many facets of Ford Instru- 
ment design and development. For more information about Ford's 
products, services and facilities, write for an illustrated folder. 
Ford engineers will be happy to discuss your problems of control 
with you. 


88 
FORD INSTRUMENT COMPANY 
DIVISION OF SPERRY RAND CORPORATION 
31-10 Thomson Avenue, Long Island City 1, New York 
Beverly Hills, Cal. * Dayton, Ohio 


Ford’s capabilities are among 
the finest in the country 


Three dimensional cams are used in 
elaborate computing devices to charac- 
terize shell ballistics, magnetic variation, 
or to solve some basic mathematical 
function. Precision in 3-D cams is of 
vital importance. Ford Instrument de- 
signed and built a unique machine that 
can produce extremely accurate cams 
from a skillfully made master. As many 
as two thousand data points are end- 
milled to set precisely the contours of 
the handcut masters. 


Equipment used for defense must un- 
dergo rigorous tests for accuracy and 
dependability in combat. At Ford, en- 
vironmental testing laboratories repro- 
duce extremes of desert or arctic battle, 
shock of warship broadside, salt fogs and 
heavy seas. When flaws have been 
detected and corrected, equipment is 
okayed for volume production and use 
throughout the armed services. 


| 
Typical of Ford Instrument’s 40 years 
of experience in precision control is its 
work in the field of nuclear power. The 
Company, for example, is building the 
control rod drive mechanism for the Sea- 
wolf, second atomic submarine. Reactor 
designs, sensing mechanisms, control 
equipment and systems, nuclear calcula- 
tions, and other specialized equipment 
and abilities are oftered by the Company 
to this expanding industry. 
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the formation of a separate 


Missiles Engineering Department 


Growing importance of missiles 
in the nation’s defense has led to the 
separation of missiles engineering 
from aircraft engineering functions 
at Douglas Aircraft Company. 


Leadership in this important field 
has been won by Douglas in 14 years 
of development and design of 
guided missiles for the Armed Forces. 
Douglas is currently engaged in 
eight major missiles projects, under 
contracts from the Air Force, 

Army and Navy. 7 


Formation of the new department 

at Douglas opens new opportunities 
for engineers and scientists interested 
in the missiles field. Write to: 

E. C. Kaliher, Engineering Personnel 
Manager, Missiles, Douglas Aircraft 
Co., Santa Monica, California. 
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Even temperatures, 
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RCA 


ANNOUNCES OPPORTUNITIES IN 
MISSILE TESTING INSTRUMENTATION | 
AT THE AIR FORCE MISSILE TEST CENTER q 


The World's largest testing laboratory and range 
extending from the central east coast of Florida 
to the Mid-South Atlantic 


TRACKING TECHNIQUES INCLUDE 
VARIOUS PHASES OF 


Radar 

Optics 

Timing 
Telemetry 
Communications 


1. Instrumentation Planning 

2. Systems Engineering— Military Applications 

3. Installation and Maintenance of 
Instrumentation Equipment 


Benefits to you include: prestige of association with the world leader in electronics « 
interesting assignments in missile testing « professional advancement « pleasant 
Florida living « liberal Company benefits « relocation assistance. 


Our new programs provide opportunities for electronic engineers and scientists, 
physicists, mathematicians, and opto-mechanical engineers and scientists. You should 
have bachelor’s or advanced degree. Experience in data acquisition, transmission 
or processing desired. 


Send resume to: 

Mr. H. N. Ashby, Technical Employment 
Missile Test Project 

RCA Service Co., Inc. 

P.O. Box 1226, Melbourne, Florida 


RADIO CORPORATION OF AMERICA 
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Military coding equipment takes one pulse and 
inserts it into a delay line and in effect sends it 
over a number of paths, each of different lengths. 
Combining the output of the paths gives a pulse 
train with pulses spaced in accordance with arti- 
ficial length of the path. Ordinarily the flexibility 
of the equipment is limited by the fixed taps in the 
delay line and the accuracy is established by 
auxiliary circuitry. 

Now Admiral research has developed a com- 
pletely new type of delay line which is infinitely 


variable within its over-all capacity. It is adjust- 
able with the greatest facility for any desired 
interval. The accuracy of this line is limited only 
by the accuracy of the measuring equipment. 
Moreover, the Admiral delay line requires less 
complicated switching apparatus. Weight and bulk 
are reduced. Fewer components permit faster 
production at lower cost. Here is one more exam- 
ple of Admiral’s capabilities in the field of military 
electronics. Address inquiries to: 


Admiral 


OR P 


RATION 


Government Division 


Chicago 47, Illinois 


Look 10 Admiral ror 

@ RESEARCH e DEVELOPMENT e PRODUCTION 
in the fields of: 7 
COMMUNICATIONS, UHF and VHF, airborne and ground. 
MILITARY TELEVISION, receiving and transmitting, airborne 
ond ground. 

RADAR, airborne, ship and ground. 

RADIAC « MISSILE GUIDANCE « CODERS and DECODERS 
TELEMETERING e DISTANCE MEASURING ¢ TEST EQUIPMENT 


FACILITIES BROCHURE 
describing Admiral plants, equipe 
ment and experience sent on 
request. 


ENGINEERS! The wide scope of work in progress at Admiral 
creates challenging opportunities in the field of your choice. 
Write to Director of Engineering and Research, Admiral 
Corporation, Chicago 47, Illinois. 
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can you qualify for 


Rocket Engineering 


careers are being born right now at Rocketdyne* 


It may surprise you to know you can do this 
with or without specific rocket engine experi- 
ence! Engineering experience in heating and 
ventilating, hydraulics, pumps, turbines, 
combustion devices, controls and engine 
instrumentation are just a few of the related 
fields that could open your future at 
ROCKETDYNE. 


At its new Canoga Park headquarters, 
ROCKETDYNE is continuing and expanding 
its pioneering work in the design, develop- 
ment and manufacturing phases of the rocket 
engine business. And at its test laboratory, 
RockKETDYNE has complete research and test 
operations. ROCKETDYNE thus offers oppor- 
tunities in the complete range of rocket 
engine development from preliminary design 
to field testing. 


ROCKETDYNE not only offers the advan- 
tages of a constantly growing and expand- 
ing company, but also the stability that 
comes from its broad range and volume of 
projects...the stability of rocket engine con- 
tracts with all branches of the Armed Serv- 
ices and the guided missile industry. Though 
it is currently emphasizing the production 


of the most powerful large liquid-propellant 
rocket engines in the Western World, Rock- 
ETDYNE is also designing, building, and test- 
ing many other types and sizes. For the 
future, it is carrying on research for even 
more powerful and startling rocket engines. 


So this is your big opportunity to build a 
career that really has a future. A well-paid 
career that assures prestige, stability, and 
selectivity of fascinating work. If you feel 
you can qualify, you owe it to yourself to 
contact ROCKETDYNE today. 


*ROCKETDYNE is North American’s rocket 
engine division. It has just moved into new 
ultra-modern headquarters in Canoga Park, 
located in the beautiful West San Fernando 
Valley of Los Angeles. This area is famous 
for its fine residential sections, modern 
shopping-center convenience, varied recrea- 
tional and entertainment facilities. Any 
point in the San Fernando Valley is just 
minutes drive from the beaches, and the 
weather is pleasant all year ’round. Many 
engineers are interested in advance school- 
ing offered by UCLA, only 4 hour drive 
from our laboratory. 


THESE POSITIONS NOW OPEN AT ROCKETDYNE: 


DESIGN & DEVELOPMENT ENGINEERS 
Mechanical, Chemical, Electrical, Stand- 
ards, Structural and Stress. For rocket 
engine systems design or development. 
Turbine, pump, controls and combustion 
device experience preferred. 

DYNAMICS ENGINEERS 
To analyze rocket engine control systems 
utilizing electronic analog and digital 
computers, B.S., M.E., or B.S.E.E. neces- 
sary. Prefer advanced degree. Experience 
in servomechanisms, systems analysis 
desired. 


COMPUTER APPLICATION ENGINEER 
Application of automatic computers to 


Write Mr. Joe Roberts, Rocketdyne Engi- 
neering Personnel, Dept. 596-JP, 6633 Canoga 
Ave., Canoga Park, California. 


Fesruary 1956 


investigate new methods of numerical 
analysis. 

TEST ENGINEERS 
Experienced on engine systems, combus- 
tion devices, turbines, pumps and engine 
instrumentation. 

EQUIPMENT DESIGN ENGINEERS 
Electrical, mechanical, structural, indus- 
trial. For design of facilities, specialized 
test, and handling equipment. 


WEIGHT ENGINEERS 

ELECTRONICS TECHNICIANS 
SPECIFICATIONS ENGINEERS 
ENGINEERING DRAWING CHECKERS 


ROCKETDYNE 


A DIVISION OF NORTH AMERICAN AVIATION, INC. 
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NEW FAIRCHILD SHIPBOARD 
RADAR SYSTEM FOR AUTOMATIC 
SEARCH AND TRACKING 


Fairchild radar systems have gone to sea with missile launching cruisers 
of the U. S. Navy. 

Now operating with the fleet, the Fairchild SPQ-2 Shipboard Radar System 
shown here was developed to search out and track either missiles or 
aircraft completely automatically. And, controls have been “human- 
engineered” to facilitate operation under battle conditions. 


Color and shape coding in this new Fairchild radar assure rapid, sure 
identification of all controls. Rugged, shock-resistant construction 
protects equipment. 

Here again is proof of Fairchild Guided Missiles Division’s continuing 
leadership in design, research and production of vital electronic equipment. 
A Division of Fairchild Engine and Airplane Corporation 


FAIRCHILD 


GUIDED MISSILES DIVISION + WYANDANCH, N. ¥ 
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tig equipment designed for your particular needs. 


ccompanying the rapid strides being mgde in the 
development of rocket powered mifSiles is the 
increasing demand for greater quantities f liquid oxy- 
gen. To keep pace with this deman@ Air Products 
“on-location” oxygen generators are cgpable of fulfilling 
the vital need with unlimited quantgffies of LOX when 


and where required. 


Air Products, Incorporated, haf/a leading position in 
the design and construction offxygen generators in the 


United States. We design ang’manufacture: 


Large Capacity Tonnag 
ties of oxygen regard 


Generators for unlimited quanti- 
iss of size, purity or cycle 

and 
“Packaged” H Purity Generators, producing high 
purity oxygey’ delivered safely under pressure by Air 
Products nted liquid oxygen pumps. 


Let us Aeview your requirements and give you actual 
costsfand engineering information for oxygen genera- 
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AIRBORNE 
SYSTEMS 


DIGITAL DATA 
RECORDING SYSTEMS 


ACCELEROMETERS 
PRECISION 
POTENTIOMETERS 


GYRO 
INSTRUMENTS 


Scientific achievement in the field of aircraft and + 

missile development is literally attaining new horizons © e ee 

... Keeping pace with those milestones of progress, | n rn in 


Giannini, too, is advancing—both in theory and in devel 


practice. the pe 
Reducing the art of high speed high altitude instru- FOREIGN SALES OFFICES: . « 
ment design to a science, Giannini engineers and MILANO, ITALY, Scient 
scientists apply knowledge acquired through years of ppb appease ideas, 
= extensive experience with commercial and military Giannini Limited, 31/45 Gresham Street Amor 
aircraft instrumentation. When flight becomes higher of hig 
; and faster, the intensive research policy and progres- REGIONAL SALES OFFICES: dmul 
ve sive design practice at Giannini will continue to NEW YORK 1, N.Y., Empire State Bldg., CHickering 4-4700 e 
° ° CHICAGO, III., 8 So. Michigan Ave., ANdover 3-5272 optim 

| provide solutions to the ever increasing problems of PASADENA, Calif. 918 E. Green St., RYan 1-7152 

a the aircraft industry. MONROVIA, CALIF., 1307 So. Myrtle Ave., ELliott 9-5381 The wv 
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Engineering positions are open at fi 
several Giannini locations for career- G. M. GIANNINI & CO., INC., PASADENA, CALIFORNIA é 
: minded young men—write for details. i 
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SCIENCE AND ENGINEERING 
AT LOCKHEED MISSILE SYSTEMS DIVISION 


THE CREATIVE APPROACH TO MISSILE SYSTEMS FLIGHT TEST ANALYSIS 


The work of the Aerodynamics Department of Lockheed 
Missile Systems Division requires a creative approach to flight 
test planning and analysis. Consequently, research as well as 
development studies are carried on in flight determination of 
the performance, stability, control, flutter, aeroelastic, 

and aero-thermodynamic characteristics of missile airframes. 


Scientists and engineers are given full scope to explore new 
ideas, develop new experimental and evaluation techniques. 
Among projects of present interest are the development 

of high-performance free-flight models and other advanced 
simulation techniques and the accompanying determination of 
optimum flight plan and instrumentation system characteristics. 


The whole spectrum of flight test activities is covered : 


1 Application and improvement of experimental planning 
techniques, including use of probability and statistical 
theory to improve test result, accuracy, reliability, and use- 
fulness, and to decrease the’expenditure of time and money. 


2 Determination of range and precision requirements of 
systems for optimizing results in obtaining aerodynamic, 
structural and thermal parameters. 


3 The development of mathematical and physical analogs for 
prediction and data for missile performance, control, 
flutter, aeroelastic and thermal studies. 


Those possessing a high order of ability applicable to this 
field of endeaver are invited to write : 


MISSILE SYSTEMS DIVISION 


research and engineering staff 


LOCKHEED AIRCRAFT CORPORATION 
VAN NUYS, CALIFORNIA 


| 
| 
| 
A 
fee 


this new 


Aa the priltie producer of this advanced rocket fuel — 
| Westvaco has compiled the data you will need for the 
“Storage and Handling of \sym-Dimethylhydrazine” 
and will be pleased to sen copy. 
“Commercial quantities of | MH are immediately available. 
‘Your orders and requests {or information of a more 
fundamental nature wil) reccive prompt attention. 


DPEVELOPMENT DEPARTMENT 


Westvaco Chior-Alkali Division 
FOOD MACHINERY AND CHEMICAL CORPORATION 
" General Offices : 161 East 42nd Street, New York 17 

CHARLOTTE, N.C. * CHICAGO, ILL. « DENVER, COLO. + PHILADELPHIA, PA. » ST. LOUIS, MO. 


FOOD 
AND CHEMICAL 
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VOLUME 26 NUMBER 2 


Editor-in-Chief 


MARTIN SUMMERFIELD 


AMERICAN ROCKET SOCIETY 


A Comparative Investigation of a Homogeneous Combustion 
Chamber With a Two-Stage Combustion Chamber 


E. A. DEZUBAY! 


Westinghouse Research Laboratories, East Pittsburgh, Pa. 


This analysis compares the performance of two idealized 
combustion systems. In the first system, the fuel and air 
react chemically in a single thoroughly mixed, constant 
flow, combustion chamber. In the second system, the 
fuel reacts initially in one chamber with only a fraction of 
the total air, with the addition of the remainder of the air 
completing the reaction in a separate chamber. If the air 
is introduced in two such definite “‘stages,”” the chemical 
kinetics of the reaction permit a more complete utilization 
of the available chemical energy than if the fuel and air are 
permitted to react in one homogeneous volume. It is 
shown that the two-stage reaction is preferred to the 
single-stage reaction. 


Introduction 


- ees analysis presented in this report contains a partial 
solution of the idealized problem: Can a combustion 
chamber or burner composed of two homogeneous reactors, 
operating sequentially, exceed the performance of a single 
homogeneous reactor operating at the same over-all air fuel 
ratio? 

This problem is not limited to a constant reactor tempera- 
ture and fixed composition as the one treated by Avery and 
Hart (1).2. The pertinent equations are developed for a 
“staged” reactor; that is, one in which fresh air and partially 
burned products from a previous burning zone are mixed and 
permitted to undergo further reaction. An analysis is pre- 
sented for a two-stage reactor in which the first stage or pilot- 
ing zone is assumed to operate at several air fuel ratios, the 
products of which, on exit from this stage, are mixed with 
sufficient fresh air to produce an over-all air fuel ratio of 
twice stoichiometric in the second stage. The reaction con- 
tinues in the second stage. A calculated comparison of the 
heat release rates and combustion efficiencies is made for 
various air flow distributions of a burner having equal volumes 
in the first and second stage. Detailed sample calculations 
are given in the appendix. 


Analysis 


The model used in comparing a single homogeneous re- 
actor with a two-stage reactor is shown in Fig. 1. In the 
single homogeneous reactor, flow conditions are assumed to be 
such that the mixing of incoming gases with the contents of the 
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reaction chamber is instantaneous and, as a result, the com- 
position of the reactor is invariant with time or position. 
In the two-stage reactor, the first stage also has instantaneous 
mixing of the contents and incoming gases. The partially 
reacted gases, emerging from the first reactor, are assumed to 
mix without reaction in the mixing zone and then enter the 
second homogeneous reactor, again, with instantaneous mix- 
ing. This intermediate mixing zone will be assumed to be 
infinitesimally small. 


m, Moles Fuell 
m, Moles Air | 
Initio! Cond. 
To, Po, of Uniformity 
at Temp. T, Po 


Volume 


L Single Homogeneous Reactor. 


m, Moles Fuell | 
ma Moles Air | | 
(primary) | 
Reactor 
eactor Reaction 
Volume Ve 


2. Two Stage Reactor. 


Fig. 1 Analytic methods 


The performance of these two models will be judged by 
heat release and combustion efficiency, for the same volume 
of combustion space and the same over-all fuel air ratio, 
that is V = Vi + V2 and (m,/m,) = (ma + m,2)/m,;. Air 
fuel ratios symbolized by ¢ will denote the molal air fuel 
ratio. Specifically, ¢ without subscript will be the air fuel 
ratio in the first model of the homogeneous reactor; ¢, will 
denote the stoichiometric air fuel ratio; ¢ will be the air fuel 
ratio of the first stage; while ¢» will be the over-all air fuel ratio 
of the burner which is also the air fuel ratio of the second 
stage. For all the conditions considered in this discussion, the 
following air fuel relations will be fixed 

= = 24, 

In order to simplify the general approach to the problem, a 
discussion of the temperatures attainable from chemical 
reactions will be made. In Fig. 2, the temperature rise for 


various fuel air ratios is shown. The dotted linear relation 
labeled as ideal represents the temperature which would be 
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Fig. 2 Temperature fuel-to-air ratio relations 


attainable if all the chemical energy available in the fuel 
could be utilized. In actuality this condition is approached 
only when the fuel concentration is quite low relative to the 
stoichiometric fuel concentration. In the rich or near stoi- 
chiometric range, the flame temperatures that are attainable 
are considerably less, because the chemical energy of the fuel 
cannot be utilized completely due to the dissociation losses 
which prevent the complete conversion of the fuel to waste 
products. 

In order to permit the inclusion of the complicated dis- 
sociation losses, especially at fuel rich mixtures, without pro- 
hibitively tedious algebraic manipulations, a fuel availability 
coefficient is introduced. The fuel availability factor y will 
be defined as the ratio of the adiabatic flame temperature rise 
to the ideal temperature rise (the ideal temperature rise being 
defined as the enthalpy of combustion (H) divided by the 
product of the specific heat (c) and mixture-fuel ratio (¢ + 1)); 
so 


T; — To _ + — To) 
(T; — To): H 


In an actual reactor a temperature rise less than the adiabatic 
flame temperature rise (7', — To) isto beexpected. To obtain 
a temperature rise less than the adiabatic one, only a certain 
fraction y of the available fuel (ym,) is considered to react with 
the oxidant, and any of the unreacted fuel and oxidant are 
assumed to act only as diluents in any homogeneous region. 
This does not preclude the possibility of the unreacted fuel 
and oxidant burning in any subsequent reaction zone. 

The basic equation for a flow reaction implies that the 
product of the fraction consumed y times the mole rate of 
supply of the available constituent ym, is equal to the reaction 
rate r, of the constituent times the volume V of the reactor, or 


y= 


where the reaction rate 
E 
where 
A = collision frequency 
E = activation energy 
R = gas const 
T = gas temperature 
C, = fuel concentration 


Co = oxygen concentration 


The order of reaction would be the sum of the exponents, 
(2b), which for the sake of simplicity will be considered to be 
exactly equal to two, although values between 1.4 and 1.8 
(Refs. 2, 3, 4) would probably be more realistic. 

In considering the concentration of both the air and the 
fuel the following assumptions will be considered. 

1. In the process of combustion the number of moles does 
not change. 

2a. The fuel and air combine at the prevailing available 
air fuel ratio. This assumption is probably quite good at 
rich and near stoichiometric fuel-air ratios. 

2b. The fuel consumed, also consumes a stoichiometric 
number of moles of oxygen (moles of air + , where 7 = 4.76 
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moles of air/mole of oxygen). This condition is possible 
only on the lean side and would be valid only in the lean 
region of ¢/¢s greater than about 1.5. 

With these assumptions the concentrations in the same 
units as the inlet density yo become 


+4) 


is the fuel concentration available for reaction at reactor 
conditions implied by 7. 


C; 


Tovobi 
Tin(1 + 


is the oxygen available for reaction in accordance with 
assumption 2a. 


( bs ) 


is the oxygen concentration available for reaction with the 
fuel in the burner using assumption 2b. 

By combining the reaction flow Equation [I] with the re- 
action rate Equation [II] and substituting the appropriate 
concentration terms, the fuel supplied can be expressed as 


Vry yo? 1 — ( 
yw y T ( exp 


Co (1 — yivi) 


in the lean range of fuel air ratios, and as 


1 — ( =) 
= xp ——)}....[I 
my (1 — yy)\ exp RT [IV] 


for the rich range of fuel air ratios, where in each case the 
energy equation would yield 


_T = + IMT 
= Hy 


From Equations [III] or [IV] and [V] the fuel consump- 
tion at a fixed air fuel ratio of a given reactor can be deter- 
mined for either a desired temperature T or a given conversion 
fraction y. It is now necessary to extend this treatment to 
include the second stage. At the exit of the second stage the 
maximum fuel conversion possible is Yom, where y-2 is the 
availability fraction previously defined and applicable to the 
over-all air fuel ratio ¢2. The fuel converted at the exit of the 
first stage or entrance to the second stage is y:¥,:my; hence 
the fuel that could be burned in the second stage would be 


(yo — 


Letting y2 represent the fraction of available fuel consumed in 
the second stage, the flow reaction equation becomes 


E 
yl ve — = Where ry = AC Coe (cxp 


The fuel concentration is 


yoT'(1 — ye) 
) 

the oxygen concentration, now assuming a lean fuel air ratio 
(over-all), is 


Cre 


then 
nyxX1 + ¢2)? RT, 
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where an energy balance would yield 


+ 1)T2 — — — (1 + 
— wr) 


The condition of compatibility relating the first and second 
stage can be obtained from the fuel supply rate (m;). The 
fuel supplied to the first reactor is calculated from Equation 
{III} or [IV] depending on whether the mixture is lean or 
rich, and represents a fuel consumption fraction changing 
from zero to y:y1 in the first stage. Similarly, Equation [VI] 
represents the fuel initially necessary for a reactor with a 
change in the fuel consumption fraction from y:y; to y2(Y2 — 
yi) in the second stage. Since the amount of fuel (m,) 
which is supplied to the reactor is a fixed quantity, (m,) 
calculated from [III] or [IV] and from [VI] must be equal. 
As a result, the conditions of the first and second stage are 
related as 


+ gu)? wo (7) {ew RT, 


Hence the operating condition for the second stage specified 
by o and 7’, or ye is determined by the operating condition 
of the first stage ¢; and 7’; or y. 

The quantity of greatest interest to engine designers is the 
quantity defined as the specific heat release (S), defined at this 
point as the heat released per unit time per unit volume unit 
density squared; that is 


[VET] 


2 


=m 
Hence the dimensionless heat release (c) for a single-stage 
reactor being supplied with less than a stoichiometric amount 
of fuel would be 


Sn To\? _ os 


and for a two-stage reactor 


Vi + Vol + To y 
Yobs E 


Method of Calculations 


Calculations were made for a series of two-stage reactors, 
each stage of which was assumed to have the same volume. 
The over-all air fuel ratio was arbitrarily specified as twice 
stoichiometric (that is, a molal air fuel ratio of 47.8 or a 
mass air fuel ratio of 31.4 for air propane mixture). The 
primary molal air fuel ratio was permitted to vary from 47:8 
to 15.3. The former collapsed the two-stage reactor into a 
single stage, while the latter represented quite a fuel-rich 
condition in the first stage. 

The method of calculation consisted of the following steps: 

1. For a given air fuel ratio in the first stage, evaluate the 
maximum fuel rate possible by maximizing Equations [IV] 
and [V]. This procedure will produce a value of temperature 
T,,, and combustion fraction y1,,. 

2. Evaluate 7; and y2 from Equations [VII] and [VIII]. 
This involves a trial-and-error solution or graphical solution. 

3. With the values of 7. and y2 obtained in step 2, solve 
for ¢ in Equation [XI]. This procedure will yield the heat 
release for a two-stage reactor, the first stage operating at an 
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air ratio of ¢ and the second at ¢. For the steps 1 through 3 
the first stage is operating at the verge of blowout since the 
maximum possible fuel rate is being used. These maximum 
points are shown on Fig. 3 with circled points. The com- 
bustion efficiency ¢ is defined as 


e = H/c(T: — ToX1 + o2).......... [XII] 


and is used as the abscissa on Fig. 3. 

4. For points at other than maximum fuel rate conditions, 
the temperature of the first reactor is chosen such that Ti, > 
T; > Tim,that is, the first reactor temperature is chosen in the 
range between the adiabatic flame temperature and the 
blowout temperature. Once 7, is chosen, y; can be calcu- 
lated and steps 2 and 3 can be performed. As many points 
as are desired can be evaluated in the range of 7,7 > 7; > Tim. 


is CURVE 421520 
oe Two Stage Reactor Ne 
Overall Air Fuel Ratio 
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!st Stoge Air Fuel Ratio (%,) 
Moles Fuel as Shown 
2 8}+— 
| $,* 23.9 (Stoichiometric 
i \ \ 
c 
© 5x10 
\ 
4 
$, 731.9 
3 
2 \ 
= 
30 40 50 80 90 100 


60 7 
Combustion Efficiency « 
Fig. 3 Behavior of two-stage reactor 


On Fig. 3, variation of heat release with combustion 
efficiency is shown. The lean air fuel ratios of 47.8 and 31.9 
show a consistent reduction in heat release with increasing 
combustion efficiency. It should be pointed out that the air 
fuel ratio of 47.8 is actually a homogeneous reactor visualized 
as a two-stage reactor with no air admission to the second 
stage. The richer air fuel ratios of 23.9 and 15.3 show an 
inflection point in their operating range. At efficiencies 
poorer than the inflection point, the air admission to the 
second stage effectively chills the mixture, while efficiencies 
higher than the inflection point correspond to appreciable 
reaction occurring in the second stage. 

On all the curves, the heat release at the lowest efficiency 
indicates blowout in the first stage (circled points on Fig. 3). 
If the curves are compared in order of descending air fuel 
ratio, it can be seen that greater heat releases are obtainable, 
until the stoichiometric air fuel ratio is obtained in the first 
reaction stage. Decreasing the primary air fuel ratio re- 
sults in appreciable chilling as depicted by the operational 
curve of 15.3 molal air fuel ratio. From these calculations, it 
would seem that for a two-stage reactor the primary stage 
should operate at stoichiometric. 

Equation [VIII], repeated here for convenience, contains 
the inherent explanation for chilling of the reaction. If 
V, = V2, then 


— 1) 
(1 + Ti? 


(1 RT, = B = 
E 


Ye) Yrbs ~ RT2 
If the left-hand side of Equation [VIII] with the appropriate 
value of y,; (Equation [V]) is plotted, the left side of Fig. 4 is 
obtained. Similarly plotting the right side of Equation 
[VII] with the appropriate value of 7, (Equation [VIT]}) 
results in a family of curves with parametric values of 4. 
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Fig. 4 Reaction chilling in two-stage reactor as explained by 
Equation [VIII] 


The temperatures shown on Fig. 4 are defined in this man- 
ner: 1's; is the adiabatic flame temperature attainable with 
air fuel ratio ¢2. T'x is the temperature which the mixture of 
gases discharged from the mixing volume would have just 
at the entrance to the second homogeneous reactor stage. 
Tx would also be the temperature at which the total mixture 
would be discharged if no reaction occurred in the second 
stage. If the first reactor would be operating at its maxi- 
mum capacity of 70.8, the exit temperature T, would be 
determined by point @ on Fig. 4, since this is the only value 
of the right side of Equation [VIII] which would satisfy the 
left side for a given value of y; (0.8 for this point). This 
temperature is very close to 7'» and indicates that only a tiny 
bit of reaction has occurred in the second stage. Such a con- 
dition would result in chilling the reaction. If the first stage 
is loaded less strenuously with y,; in the neighborhood nine 
tenths, three solutions are possible depicted by points @)a, @b, 
and ()c. Such a triple solution gives rise to the inflection 
points mentioned previously in conjunction with Fig. 3. 
This type of solution could easily be the driving force of os- 
cillations since a relatively large variation in temperature 
could occur at random. Such a random temperature fluctua- 
tion would cause a marked change in the pressure drop across 
the combustion chamber due to the variation in the momen- 
tum change. Such a triggering mechanism could easily sus- 
tain pressure oscillations in the cavity represented by the 
burner. 

Point (8) of Fig. 4 represents a relatively modest operating 
condition for which the first and second-stage reactions are 
well completed. 


Conclusions and Recommendations 


From this analysis it is indicated that, by staging, heat re- 
leases higher than those of a homogeneous reactor can be 
obtained. The flow conditions must be regulated carefully 
so that chilling in the second stage does not occur. It is 
hoped that this two-stage reactor treatment can be expanded 
to many stages, and the optimum air flow per stage can be 
found. Experimental verification of the results indicated 
in this analysis is planned in the near future. 


APPENDIX 
Sample Calculation 
Part 1 Constants to Be Used 


Activation energy, EF = 33 kcal/mole (Ref. 5) 
E/R = 29,900°Rankine 
Heating value (lower) of propane = 878,000 Btu/lb-mole 
Average specific heat = 7.34 Btu/lb-mole 
Initial temperature Tp = 600°R 


¢ moles 
T, (Ref. 6) T; —To air/mole fuel v 
4090 3490 23.9 0.728 
3015 2415 47.8 0.978 
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Part 2 Calculations 


In step 1 of the description of the calculations, Equation 
[IV] must be maximized. Rewriting this as 


E 
log. (1 — — 2 log, T, — RT, 


will be a maximum when 


(tos myn (1 + 


or VA yo?T 
which will be when 
( 1 v1 1 ) 2 E 
— yim 1 — Yim Tia Bie 
but 
> = : since Te — 
aT, T;—T 
so that 


T; - — W(Tin — T) 
E 


Rr.” 


the solution of which is 7;,, = 3310 and y:,, = 0.776. Equa- 
tion [IV] then becomes 


me _ Vis Te me, _ 
A + 1)? Yim ( 
which with the appropriate numerics yield 
23.9 ( 600 \?0.224 
= 1.895-10-V; 


which by Equation [VIII] is related to the second stage by 


E 
Vode —pr 

1.895-10*V, = = 

(oe + DAT Ye 

(: = {v2 ws} ) 

which with the appropriate numerics and with V; = V; 
becomes 

(48.8)? (2 21 — 
.895-10-* ———- = 
1.895-10 78 e 2 Ts 


(: {0.978— 0.728 X 0.76 


29,900 
0.944:10% =e (@) 
ye 


and from Equation [VII] we have 
48.8T, — (47.8 — 23.9)600 — 24.9(3310) 


(1 — 0.207y:)....... [1] 


878,000 
734 (0.978 — 0.565) 
T: — 1985 
2 
1030 (2) 
Substituting [2] in [1] we have 
29,900 /600\2/3015 — Ts 
- 
(Continued on page 97) 
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Flame Spreading From Baffles 


RICHARD J. PETREIN,' JOHN P. LONGWELL,’? and MALCOLM A. WEISS® 


Esso Research and Engineering Company, Linden, N. J. 


Flame spreading from 0.5 to 2.25-inch V-gutter baffles 
in two-dimensional turbulent ducted flow was studied 
experimentally. Measurements of combustion efficiency 
were made at distances of 10 to 43 inches downstream of 
the flame stabilizing baffle(s) in a 9 X 5-inch duct. The 
effects on efficiency of several geometric and operating 
variables were observed. Little effect of baffle size was 
noted, but efficiency increased when reducing the spacing 
between baffles (or between baffle and wall) by adding more 
baffles. Pressure dependence was variable, Efficiency 
increased with pressure (over the range of !/; to ‘/; atm abs 
tested) when using geometries giving high flame spreading 
rates, but there was little change when using geometries 
giving low rates. Smaller effects were found on changing 
equivalence ratio (faster burning occurred nearer stoichio- 
metric mixtures) and velocity (efficiency decreased when 
velocity increased). Miscellaneous data are noted for 
changes in tailpipe length, fuel vaporization, and nature of 
turbulence. An empirical correlation is presented which 
groups together in a simple way most of the data taken. 
Many variables interact complexly in flame spreading and 
no accurate simple description seems likely. Thus, 
neither the concept of a wrinkled laminar flame, nor of an 
extended reaction zone, nor of homogeneous combustion 
fits the data generally—although any one concept may be 
acceptable over a restricted range. Burning is probably 
discontinuous in micro-time making inapplicable theories 
assuming continuous burning. 


Introduction 


OME ramjet engines use baffles as flame stabilizers. 

These baffles are aerodynamically bluff bodies. In the 
sheltered eddying region behind the baffle, a flame can be 
stabilized, thus providing an ignition source for the total 
stream. A ramjet engine can be characterized by its combus- 
tion stability and efficiency. In any baffle engine, com- 
bustion stability means simply the range of conditions over 
which a flame will be stabilized in the baffle eddy region; 
combustion efficiency, on the other hand, depends on the rate 
at which the stabilized flame spreads through the remainder 
of the air and fuel passing through the combustor. There is 
relatively much information available on the stability of 
baffles; stability studies have been made in which both 
baffle geometry and operating conditions were widely varied. 
Critical surveys of these data exist, e.g. (1, 2, 3).4 However, 
there is little quantitative information available on the rate of 
flame spreading from stabilized ignition sources in ducted 
streams. Miscellaneous measurements in relatively small 
equipment have been published by Mullen, et al. (4), Berl, 
et al. (5, 6), Wilkerson and Fenn (7), and Williams, et al. (8). 
The purpose of this paper is to report on some experiments 
on flame spreading from baffles in larger equipment with 
several geometries and operating conditions of practical 
interest. Some discussion of the data is included. 
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Il Equipment and Operation 


In all cases, the flames investigated were stabilized by two- 
dimensional V-gutter baffles, 5 inches long, with total apex 
angles of about 35 to 55 degrees and with bases ranging from 
0.5 to 2.25 inches wide. These baffles were mounted normal 
to the stream and parallel to the short side of a 9-inch by 5- 
inch rectangular duct sketched in Fig. 1. Single baffles 
were centered in the duct and multiple baffles mounted sym- 
metrically about the center in a single plane normal to flow 
such that the distance between baffle centers was twice the 
distance between the wall and the centers of the baffles nearest 
the walls. 


BAFFLE RECTANGULAR 
STABILIZERS 
HOMOGENEOUS QUENCH 
AIR-FUEL 9" VALVE," 
MIXTURE =| AND 
EJECTOR 
CHOKING MOVABLE COOLED 
PERFORATED PROBE ANNULAR 
PLATE BAFFLE 


Fig. 1 Sketch of flame spreading apparatus 


About 21 inches upstream of the baffles, a choking per- 
forated plate was used to control velocity distribution. The 
fuel was solvent naphtha, a paraffinic volatile petroleum 
fraction with properties tabulated in (9). Fuel was injected 
upstream of the perforated plate into the 300° F air stream. 
This procedure gave a homogeneous air-fuel mixture at the 
baffle; mixture temperature was about 235° F, varying 
somewhat with equivalence ratio. Except when equivalence 
ratio was the specific variable under study, stoichiometric 
mixtures were always used. Air humidity was held constant 
at 0.025 lb water/ib air. 

Downstream of the baffle, the duct extended for about 
five feet; at this point there was an abrupt expansion into a 
water quench section followed by piping through a butterfly 
valve to a steam ejector. The valve was used for static 
pressure control, the control point taken as just upstream of 
the baffle. Static pressure at the end of the burning zone 
was never more than 10 per cent below the control pressure. 
A removable water-cooled annular baffle (41 per cent block- 
age) was located near the end of the duct (just before the 
expansion) during most of the tests. This baffle was used at 
operating pressures greater than about 1 atm abs to reduce 
pressure fluctuations while burning; at lower pressures, it 
could be removed without affecting operability. 

At distances of 10 to 43 inches downstream of the test 
baffle(s), provisions were made both to measure wall static 
pressures and to insert a movable water-cooled gas sampling 
probe in the stream. The probe traverses were used to deter- 
mine local values of gas velocity and oxygen consumption 
efficiency (calculated from measured oxygen concentrations 
before and after completing combustion). By using local 
mass velocity to weight local efficiency, an over-all weighted 
average efficiency at the sample plane was calculated. Effi- 
ciency is defined as the oxygen consumed divided by the 
oxygen that would be consumed were combustion complete. 

The chief dependent variable in the work was the weighted 
average combustion efficiency in a plane normal to the flow 
at various distances of this plane from the stabilizer. As 
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noted above, the average efficiency was calculated by weight- 
ing measured local efficiency by measured local mass veloc- 
ity. However, making complete velocity and efficiency 
traverses in any survey plane is quite tedious. It was found 
that the drop in wall static pressure (relative to the control 
point) could be correlated with average combustion efficiency 
in the plane of the static pressure measurement. Thus, after 
initial correlations, these easily measured wall static pressures 
were used to both interpolate and extrapolate average 
efficiencies determined from velocity and efficiency traverses. 
Efficiency errors probably averaged about +5 per cent in 
absolute efficiency. 

In these tests, Reynolds number (based on baffle width and 
approach stream conditions) varied from about 2 X 10¢ to 
12 X 104 with most runs at a Reynolds number of 4 X 104 or 
higher. These Reynolds numbers are above the value where 
analysis (2, 3) of several sets of experimental data shows a 
transition in the flame stabilizing mechanism. This transi- 
tion value—of about 1.5 X 10‘—is probably sensitive to the 
nature of the approach flow. In the experiments of (2, 3), 
turbulent intensity was quite low—about 1 per cent—whereas 
in these tests, turbulent intensities averaged about 7 per cent. 
It is likely that the transition Reynolds number occurs at 
lower values for a more turbulent flow. Thus, the tests 
reported here were probably safely past this transition region 
and in a regime of fully turbulent stabilization—and propaga- 
tion, if a similar transition for propagation exists. 


III Experimental Results 


Examination of some detailed traverses is helpful before 
considering effects on over-all average efficiencies. Figs. 2 
and 3 present typical efficiency and velocity traverses for the 
conditions listed. The traverses were made in the center of 
the duct in a direction normal to the flow and the baffle. 
As one proceeds downstream, over-all efficiency and velocity 
levels increase. However, it is of particular interest that, 


a 43" DOWNSTREAM 

w 

60 

& 40 

z 

8 20 


-2 ° 2 3 4 
DISTANCE FROM DUCT CENTER LINE -INCHES 


Fig. 2 Oxygen efficiency profiles downstream of a 2.25-in. V- 
gutter. Pressure = 10’ abs, average velocity approaching 
baffle = 170 ft/sec 
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Fig. 3 Velocity profiles downstream of a 2.25-in. V-gutter. Oper- 
ating conditions as in Fig. 2 


for these time-average measurements, no steep gradients of 
velocity or efficiency are found and all levels (between the 
extremes) exist. Over-all average efficiencies at the 10, 25, 
and 43-inch stations are 9, 24, and 58 per cent, respectively. 
The so-called “flame spreading curve” plots these average 
efficiencies as ordinates against distances from the baffle as 
abscissas. The complete curve is almost always S-shaped 
with slope (the rate of flame spreading) increasing gradually, 
as efficiency increases, to a maximum at efficiencies of about 
50 to 80 per cent and then falling off again as the efficiency 
approaches 100 per cent. 

It is sometimes convenient to compare local flame spreading 
rates by converting flame spreading curve slopes to effective 
turbulent flame velocities (S,). Flame velocity keeps its 
usual definition as the volumetric rate at which unburned 
mixture is completely burned divided by the apparent flame 
front area. A flame front area is calculated by assuming the 
existence of distinct regions of completely burned and com- 
pletely unburned mixture in the duct. The interfaces be- 
tween these regions are zero-thickness flame fronts. Such a 
model was used, for example, by Tsien (10). Although this is 
a convenient model for calculating a fictitious flame velocity, 
it bears no resemblance to the physical reality illustrated by 
Fig. 2. For the two-dimensional systems used here, one finds 
that, for the model above 


dz 

in consistent units where uw is the average stream velocity 
approaching the baffle, h is the distance between baffle center 
and duct wall (or plane of symmetry between two baffles), 
E is the average fraction of oxygen consumed in a plane at a 
distance x from the baffle, and @ is the angle between a tangent 
to the flame front and the mean flow direction. Cos @ can 
usually be taken as unity with little error. For the data of 
Figs. 2 and 3, effective turbulent flame velocities calculated by 
Equation [1] were 4, 12, and 16 ft/sec at the 10, 25, and 43- 
inch stations, respectively. 

The effects of specific variables on the flame spreading 
curve can be examined next. Flame spreading is an ex- 
tremely complex process; many interacting variables operate 
and it would be very unusual if the effect of any particular 
variable were simple and independent of other variables over 
any appreciable range of conditions. 


Baffle Size and Blockage 


There appears to be no consistent effect of baffle size or 
blockage (the percentage of cross-sectional area occupied by 
baffles) in the range tested. Illustrative data are shown in 
Fig. 4. Within the accuracy of the data, single baffles 0.5 to 
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Fig. 4 Effect of baffle size for single and multibaffles at several 
pressure levels. Approach velocity of 170 ft/sec unless noted; 
parameter is baffle width in inches 
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9.25 inches in width (5.5 to 25 per cent blockage) perform 
alike at 20 to 22” Hg abs. Results are similar for single 1 
and 2.25-inch baffles at 10 to 11” Hg and for three-baffle 
arrays of 0.75-inch and 1.5-inch gutters (25 and 50 per cent 


blockage) at 25” Hg. The implication of these results is that 
the size of the ignition source—if large enough to cause any 
flame spreading—does not greatly affect the spreading rate. 
It was also concluded in (8) that the size of small (less than 
0.5-inch) stabilizers did not affect flame width downstream of 
the stabilizer in a duct. 


Baffle Spacing and Duct Size 


Spacing between baffles (or between baffle and wall) was 
altered by simply changing the number of baffles tested. 
Some results are shown in Fig. 5. Efficiencies at the 10-inch 
station for one, two, and three 1.5-inch baffles are about 10, 
40, and 60 per cent, respectively. If there were no interaction 
between baffles, one would expect the efficiency at this point 
to double in the case of two baffles or triple in the case of 
three baffles. However, the increase is about fourfold for 
two baffles and about sixfold for three baffles and, therefore, 
an interaction does exist which gives improved rates of com- 
bustion. 
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Fig. 5 Effect of number of baffles of two sizes. Approach 
velocity of 140 to 150 ft/sec 
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There are not enough data for the region in the immediate 
wake of the baffles to shed light on flame spreading at very 
close distances. This effect of number of baffles is apparently 
restricted to regions of lower average efficiency since the rate 
of increase of efficiency (i.e., flame velocity) is about the same 
(at a given efficiency), regardless of number of baffles, at 
efficiencies higher than about 50 per cent. This is surprising, 
since one might expect that a larger number of ignition sources 
would give faster burning throughout the entire combustion 
zone. Examination of the velocity profile for the multi- 
baffle arrangements shows an almost flat velocity distribution, 
whereas the single baffle shows velocity gradients (cf. Fig. 3). 
Since velocity gradients are expected to result in higher mixing 
rates, this effect may counterbalance the advantage of more 
ignition sources. Thus, the same rates of combustion would 
result in the latter stages. 

A question arose as to whether a physical duct wall would 
have the same effect as a plane of symmetry between baffles. 
A divider was thus inserted into the 9 X 5 duct forming two 
5 X 4.5 ducts. In each of these, a 1.5-inch baffle was 
mounted and operated at about 14” Hg and 125 ft/sec ap- 
proach velocity. Traverses 25 inches downstream of the baffle 
showed the average efficiency to be 74 per cent with or without 
the divider in place. Thus, a plane of symmetry and a duct 
wall have about the same effect—within the accuracy of 
measurement. (In most cases, small differences might be ex- 
pected due to heat losses and altered velocity distributions 
at a physical wall.) 
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Pressure 


The effects of pressure are complex. In general, as effec- 
tive turbulent flame velocities increase—due to a multibaffle 
stabilizer, or a smaller duct, or intense turbulence and mixing 
from any cause—the flame velocities become more sensitive 
to pressure. Flame velocities decrease as pressure decreases; 
the pressure exponent of flame velocity variation ranges 
from about 0.1 to 0.6 depending on the particular geometry 
and operating conditions. 

Several sets of data illustrating varying pressure effects are 
shown in Fig. 6. When six 0.75-inch baffles are used—a 
configuration giving rapid burning rates—there is a marked 
drop in performance when operating at 10” Hg rather than at 
23” Hg. Similarly, when a single 1-inch baffle was tested by 
Ber! (5, 6) in a small duct with intense mixing, flame spreading 
rates were again very high; these high rates resulted in the 
same large pressure effect as shown by the Fig. 6 curves com- 
puted for 10” Hg from (6) and for 30” Hg from (5). On the 
other hand, with low flame velocities, the effect of pressure is 
small. Thus, Fig. 6 shows only small change when operating 
a single 2.25-inch baffle between 10 and 40” Hg and a similar 
small change for a single 1-inch baffle between 8 and 20” Hg. 
(The latter baffle showed anomalous behavior at 11” Hg; 
such anoamlies or poor reproducibility are frequent problems 
in making flame spreading measurements.) 
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Fiz. 6 Effect of pressure for several configurations. Approach 
velocity of 170 ft/sec unless noted; parameter is static pressure 
in inches of mercury, absolute 


Equivalence Ratio (¢) 


One series of tests was made in which a single 2.25-inch 
baffle was used as a stabilizer while equivalence ratio (fuel/air 
ratio divided by the stoichiometric fuel/air ratio) was varied 
between 0.75 and 1.26. The results are shown in Fig. 7. 
The fastest burning rates occurred at stoichiometric, where 
¢ = 1.00, and the slowest rates for the leanest mixture, where 
¢@ = 0.75. Curves for @ = 1.26and ¢ = 0.84 coincided. The 
trends shown are certainly in qualitative accord with other 
combustion measurements; quantitatively there is somewhat 
less influence of equivalence ratio than might be expected 
from relative chemical reactivity (as measured by laminar 
flame velocities, for example). As mixing intensities increase 
or pressure decreases, relative chemical reactivity of the fuel 
should become of more importance. For example, in Mullen’s 
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Fig. 7. Effect of equivalence ratio (¢) 


(4) high output 2-inch combustor, turbulent flame velocities 
varied as laminar flame velocity taken between roughly the 
0.6 power and the 1.8 power (averaging 1.4) for various 
conditions and geometries. Laminar flame velocity was 
changed by changing equivalence ratio and inlet temperature. 


Approach Velocity 


For one configuration—a single 2.25-inch baffle operating at 
10 to 11” Hg—approach velocities of 119, 170, and 201 ft/sec 
were tested. The results—in Fig. 8—show that as approach 
velocity increased, combustion efficiency near the baffle in- 
creased but, at longer distances from the baffle, efficiency 
decreased. Differences were not large and, from these 
limited data, a general conclusion about the effect of velocity 
should not be drawn. 
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Fig: 8 Effect of stream velocity approaching baffle. Velocity 
parameter in ft/sec 


Tailpipe Length 


Longwell (11) has presented data showing that the stability 
of baffles mounted in a duct with a short tailpipe (about 12 
inches long) is two to three times as great as the stability of 
baffles mounted in a duct with a long tailpipe (about 5 feet 
long) like that used here. This result is presumably due to 
pressure and velocity fluctuations in the long duct; these 
fluctuations at least temporarily overload the eddy zone 
behind the stabilizer and thus extinguish the flame. How- 
ever, heavier loading of this zone should then show up as an 
increased efficiency at short distances from the baffle. Some 
scattered data confirm this expectation. For example, in 
one typical test the efficiency 10 inches downstream of the 
baffle was 5.7 per cent using a short tailpipe (in equipment 
otherwise the same as described previously) and about 9 
per cent with the (usual) long tailpipe. Therefore, tailpipe 
length does affect at least the initial flame spreading. The 
average amplitude of pressure fluctuation in the long tailpipe 
system was about 50 per cent greater than in the short tail- 


pipe system. 
Turbulence 


It was observed in Fig. 6 that a 1-inch baffle tested by 
Berl in a 2 X 4-inch duct gave very high flame spreading 
rates relative to those obtained in the 9 X 5-inch duct here. 
Part of the difference was attributed to duct size, but there 
was also a difference in the choking perforated plate used for 
velocity distribution. A plate to simulate Berl’s was thus 
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made up for the 9 X 5 duct. This plate had 0.040-inch holes 
and 5 per cent open area compared to the 0.1875-inch 
holes and 10 per cent open area used in all other tests. The 
results are shown in Fig. 9 for a single 1.5-inch baffle operated 
with each of the plates. A clear advantage exists for the 
plate with smaller holes. 
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Fig. 9 Effect of hole size in choking perforated plate 


To attempt to explain the difference, turbulent intensities 
were measured by Westenberg’s helium diffusion method (12) 
just upstream of the baffle. Observed intensities were about 
6.6 to 6.8 per cent with the large-hole plate and about 7.1 to 
7.5 per cent with the small-hole plate. It is difficult to see 
how this relatively small intensity difference at a high inten- 
sity level could account for the observed differences in flame 
spreading rates. One suspects either that scale of turbulence 
differed appreciably for the two plates or that the average 
intensity obtained by helium diffusion does not give suffi- 
ciently detailed information about turbulent fluctuations. 


Fuel Vaporization 


All of the preceding data were obtained with homogeneous 
mixtures of air and fuel in the equipment described earlier. 
In a practical burner, homogeneity is rarely attained, and it is 
thus of interest to observe the effects of incomplete fuel 
vaporization in a typical system. The system tested con- 
sisted of a 1.5-inch D cylindrical pilot burner centrally 
mounted in a 6-inch D duct. A certain part of the total air 
flowing (termed the per cent pilot) was diverted through the 
pilot and burned completely. The products of combustion 
exhausted into the duct through which the remaining air was 
flowing in an unburned combustible mixture. Measurements 
were made, in a plane 18 inches downstream of the pilot exit, 
of the weighted average combustion efficiency. Conditions 
in this case were a pressure of one atmosphere, a mixture 
velocity of 200 ft/sec, and an inlet mixture temperature of 
250° F. The over-all equivalence ratio was 0.83. 

For several conditions, a comparison was made among three 
types of main-duct fuel injection. The first type used a 
homogeneous vaporized naphtha feed. The second type in- 
jected liquid naphtha 10 inches upstream of the pilot outlet. 
and the third injected Diesel oil (a low-volatility fuel with 
properties given in Ref. 9) at the same point. The main dif- 
ference in these types of injection is obviously the degree of 


Table 1 
Per cent main duct fuel burned 
Homo- 
Per Pilot geneous 
cent equiv. vaporized Liquid Diesel 
pilot ratio naphtha naphtha oil 
11 1.06 57 50 35 
11 0.94 53 47 18 
7 1.06 va 42 25 
7 0.83 5 39 22 
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vaporization of fuel attained—from 100 per cent in the first 
type to a negligible amount in the third. Some of the results 
are shown in Table 1. 

Clearly, poor vaporization can have a pronounced effect on 
flame spreading. Allowance for nonhomogeneity must thus 
be made in any practical system. 


Empirical Correlation of Data 


Examination of the preceding flame spreading data shows 
that there are two simple ways in which the flame spreading 
curve seems to be altered as a parameter is altered. The 
first is a change in the rate of flame spreading or slope of the 
curve. The second is a translation of the distance scale; 
although different rates of flame spreading may occur in about 
the first 10 inches behind the baffles, rates frequently follow 
the same curve for distances downstream greater than about 
10 inches. Therefore, a dimensionless correlating parameter, 
B, was developed empirically as the abscissa of a generalized 
flame spreading curve. The parameter consists of ratios of 
the operating variables to reference conditions—each ratio 
raised to an exponent best fitting the data. Thus 


Here, x is the distance downstream of the stabilizer (with 36 
inches as reference), wp is the approach velocity (with 150 ft/ 
sec as reference), P is the static pressure (with 1 atm abs as 
reference), and S, is the laminar flame velocity (with a ref- 
erence of 63 cm/sec, Dugger’s (13) value of the flame 
velocity of stoichiometric propane-air mixtures at 235° F); 
it was convenient to assume the S,; of naphtha to be pressure 
independent and to equal the S, of propane, thus making use 
of Dugger’s values at one atmosphere. Both Az (inches) and 
the pressure exponent, n, are functions of the flame spreading 
distance, i.e., the distance between baffle center and a wall or 
plane of symmetry between baffles. For the data of Figs. 4 
to 8, the best values of Ax and n were found to be those 
shown in Fig. 10. If these values are used to define B, 
curves of the type shown in Fig. 11 result. In Fig. 11 
appears every curve of Figs. 4 to 8 inclusive (except the data 
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Fig. 11 Empirical correlation of flame spreading data 
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of Berl). A mean curve has also been drawn in. Because 
Fig. 11 fixes efficiency as a function of B, then the first terms 
on the right of Equation [2] fix the slope of the usual flame 
spreading curve, and the Az/36 term fixes the translation. 

About a dozen other curves were determined—not re- 
ported in previous figures—for various combinations of 
variables within the ranges already cited. All of these curves 
fall within the scatter of Fig. 11. Although this scatter 
could be reduced by a more complex correlating parameter, 
the increased precision was not deemed worth the increased 
complexity. The curve of Fig. 9 for the small-hole velocity 
distribution plate, and the curves of Berl (Fig. 6) do not fit the 
correlation of Fig. 11. The former curve is translated from 
the mean and the latter curves vary greatly from the mean in 
slope. Thus, some features of the flow—which probably are 
not easy to characterize—can change both the slope and 
translation of the flame spreading curve. 

Clearly, no correlation as simple as that of Fig. 11 could 
be an accurate general description of a process as complex as 
flame spreading; therefore, extension of the correlation to 
other geometries or a wider range of operating conditions is 
not recommended. However, the curve does reasonably 
describe most of the rates observed in this study. 


IV_ Interpretation of Data 


No theory exists that will quantitatively predict the rates 
of flame spreading in turbulent ducted flow. Flame spread- 
ing, like most other combustion phenomena, is affected by 
many interacting variables. Under some conditions, the rate 
may be completely controlled by factors that would be of 
little importance under other conditions. For example, in the 
limiting case of unmixed flow, flame spreading could occur as a 
laminar flame front (or a somewhat wrinkled laminar flame 
front). For such a case, S; ~ S, ~ P®° because pressure 
usually does not greatly affect laminar flame velocity. At 
the other limit, with extremely turbulent flow and violent 
mixing, there would be a profound pressure effect and S; ~~ 
P°8 (Ref. 14) if combustion were occurring homogeneously. 
The pressure exponents observed in this study ranged from 
about 0.1 to 0.6 and thus were intermediate between the 
extremes above. Higher exponents, it will be recalled, were 
obtained with a closer baffle spacing or a smaller duct, either 
of which would give more intense mixing. 

Not only were observed pressure exponents intermediate 
between the limits theorized above but observed flame 
velocities were similarly intermediate. The data presented 
give equivalent values of S, from 5 to 50 ft/sec; if burning 
occurred as a laminar flame front under these conditions, S, 
would be 2 ft/sec, and if it occurred as optimum homogeneous 
combustion, S; would range up to 340 ft/sec. (This “opti- 
mum homogeneous combustion” assumes a theoretical system 
like that proposed by Avery and Hart (15), except that no 
physical wall separates burned and unburned mixtures, and 
uses the kinetic data of Longwell and Weiss (14). Such a 
system gives the highest flame spreading rates.) 

Although “optimum homogeneous burning” does not 
occur, it is possible that homogeneous combustion is occurring 
locally. For example, consider the efficiency profiles of Fig. 
2. One can treat any particular profile by assuming homo- 
geneous combustion in all local regions across the profile. 
Thus, in each increment of gas, local chemical composition and 
local adiabatic temperature rise correspond to the local 
observed efficiency; heat or mass transfer normal to the flow 
direction is neglected. Local reaction rates corresponding 
to the local efficiencies were taken from (14). Summing up 
the local rates (weighted for local mass velocity) over the 
profiles gives one an over-all differential change of efficiency 
with time or distance, i.e., a value of dE/dz at each station. 
These values can be substituted in Equation [1] to obtain S;. 
The results (for the three efficiency profiles of Fig. 2 and 
corresponding velocity profiles of Fig. 3) are shown in Table 2. 
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Table 2 
Si, ft/sec 
Inches Average % Experi- Local 
downstream _ efficiency mental homogeneous 
10 9 4 27 
25 24 12 57 
43 58 16 60 


The over-all flame velocities calculated from local homogene- 
ous reaction rates are much greater than experimental values. 

The computation above can be pursued for one more step. 
Allow combustion to proceed as the gases move downstream, 
again assuming that the burning is locally homogeneous. 
One can then calculate how efficiency profiles and effective 
turbulent flame velocities change while moving downstream 
from any given starting profile. The results of such a cal- 
culation are shown in Figs. 12 and 13; the starting profiles are 
those at the 10-inch station of Figs. 2 and 3. Again, radial 
heat and mass transfer are neglected, and the error so intro- 
duced increases as distance from the reference station in- 
creases. 
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Fig. 12 Calculated efficiency profiles assuming locally homo- 
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Fig. 13 Calculated flame velocities assuming locally homo- 
geneous combustion 


Note first (Fig. 13) that the very high initial value calcu- 
lated for S; rapidly falls off when proceeding downstream. 
At a distance of about 4 to 5 inches downstream, the calcu- 
lated value drops to equality with the experimental value. 
At greater distances, it drops even farther to only about 15 
per cent of the experimental values. 

As the calculated S; falls, Fig. 12 shows that the calculated 
efficiency profiles become steeper and steeper. This results 
because burning rate is a very sensitive function of efficiency. 
If the initial efficiency exceeds about 30 to 40 per cent, burning 
proceeds rapidly toward completion. If the initial efficiency 
is below about 30 to 40 per cent, the rates never do really get 
started. Thus, the rapidly burning parts of the profile 
quickly burn themselves out. Only slow-burning material of 
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low initial efficiency remains. This behavior of the profile 
accounts for the rapid drop in calculated S;. 

The calculated profile 15 inches downstream of the reference 
profile can be compared with the experimental profile at this 
station (25 inches from the baffle in Fig. 2). There is an 
obvious lack of agreement. Although the experimental 
profile is at a higher efficiency level, there are no steep effi- 
ciency gradients. Thus, the concept of local homogeneous 
burning fails to predict not only the over-all effective flame 
velocities but also the general shape of the efficiency profiles, 
Part of the failure of calculated and experimental profiles to 
agree, but certainly not the greatest part, is due to the re- 
strictive assumption of no transfer normal to the flow direc- 
tion. 

No satisfactory general mechanism for flame spreading 
has been suggested. Conceivably, for relatively smooth flow 
and low burning rates, flame does propagate as a laminar front 
with extended surface area. However, the conditions for 
which this mechanism seems applicable are quite limited. 
Observed pressure dependence and over-all burning rates 
make this theory untenable for most conditions of practical 
interest. At the other extreme, combustion does not proceed 
homogeneously even when local time-average compositions 
are taken into account. The instantaneous rates predicted by 
homogeneous local combustion are usually high; the rates 
predicted by optimum homogeneous burning are always much 
too high. In general, rates and pressure dependence for 
experimental flame spreading may have any values between 
the extremes predicted by (wrinkled) laminar flame front 
theory and by homogeneous reaction theory. 

An intermediate mechanism was suggested by Summerfield, 
et al. (16) who proposed that a turbulent flame propagates by 
means of an extended reaction zone; this zone is similar in 
principle to a laminar flame front, but the essential burning 
mechanism is not a laminar flame. An extended reaction 
zone is a useful qualitative concept and does allow a con- 
ceptual bridge between laminar and homogeneous flames. 
However, it is not clear how a quantitative application to 
these data can be made. Summerfield’s treatment suggests a 
close analogy between a turbulent and a corresponding 
laminar flame; the data cited in support of this analogy are 
quite limited in range. Furthermore, in order to compute a 
turbulent flame velocity or space heat release rate—assuming 
that all parameters of the corresponding laminar flame are 
known—one must know the approach stream turbulent 
diffusivity and the thickness of the turbulent reaction zone. 
In general, this latter quantity would be quite difficult to 
define, much less to measure in ducted high-velocity flow if 
it could be defined. 

It seems likely that combustion efficiencies at any point in 
space are not constant in micro-time. Time-average measure- 
ments showing a given efficiency probably really represent 
averages of alternate gas pockets of unburned, partly burned, 
and fully burned material. Conclusions based on a homo- 
geneously burning mixture at the time-average efficiency 
could then be grossly in error. 

If burning does consist of detached masses of burned, 
partially burned, and unburned material, one would expect 
burning rates to be particularly influenced by the nature of the 
flow. The size, nature, and transfer rates of these small 
detached masses would likely depend on stream turbulence 
and on velocity gradients -in the stream. However, the flow 
itself would depend on the rate of burning; burning causes gas 
acceleration and velocity gradients, for example. Thus, a 
feedback mechanism is created with burning rate and flow 
interdependent. The result is that burning becomes par- 
ticularly sensitive to the nature of the original flow. Further- 
more, one might expect such a feedback mechanism to cause 
oscillations. It is well known that flow and pressure oscil- 
lations almost invariably exist in ducted flow when combus- 
tion is permitted to approach completion. 

(Continued on page 97) 
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High-Frequency Combustion Instability in Liquid 


Propellant Rocket With Concentrated Combustion 
and Distributed Time Lag 


SIN-I CHENG! 
Princeton University, Princeton, N. J. 


In view of the turbulent conditions in the combustion 
chamber in a liquid propellant rocket motor, different pro- 
pellant elements will encounter different dynamic and 
heat transfer conditions before they are converted into 
burned gas. This fact is represented in the simplified 
model of analyzing high-frequency combustion instability 
assumed in Ref. (1)? by the nonuniformity of the time lags 
and the space lags of different propellant elements. The 
effect of nonuniform time lag on the high-frequency insta- 
bility of the system is analyzed with the spatial distribution 
of combustion assumed as concentrated at arbitrary axial 
location, i.e., uniform space lag. It is inferred from the 
analytical results that spreading the sensitive time lag 
about a mean value is stabilizing, but its over-all effect, 
especially for small time lag spread, is not very significant 
as compared with the stabilizing effect of the supercritical 
nozzle. 


Introduction 


HE stability of high-frequency small longitudinal dis- 

turbances in the combustion chamber of a liquid pro- 
pellant rocket with constant injection rate has been investi- 
gated in (1,2, 3). The combustion of each propellant element 
is assumed to take place instantaneously after traveling a dis- 
tance called the ‘total space lag” from the injector end during 
a period called “total time lag.” The last part of the time 
lag will vary with the variations of the physical properties of 
the gaseous medium through which the unburned element 
is traveling. This is called the “‘sensitive time lag.”” The 
configuration of the combustion system can then be schemat- 
ically described by two functions defining the distribution of 
propellant elements in different ranges of the sensitive time 
lag and of the total space lag. Owing to the complicated 
mathematical formulation, it is assumed in previous investiga- 
tions that all elements have the same sensitive time lag. 
The spacewise distribution of combustion is found to be 
very important in determining the stability. It is therefore 
necessary to investigate the effect of the fact that different 
elements have different sensitive time lags. 

In the present paper we shall study the system with all 
elements having the same total space lag but different sensi- 
tive time lags; that is, a system with a single concentrated 
combustion front but the propellant elements burning at the 
instant ¢ became sensitive to the physical states of the sur- 
rounding gases at different instants ¢ — 7. 

Objection has been raised as to the feasibility of this as- 
sumption. It should be noted that, with the complicated 
patterns of the injection velocity and the recirculating motion 


This work was supported or Bureau of Aeronautics under 
Contract no. NOas 53-817-c. This pon is condensed from 
Princeton Aero. Engng. Rep. no. 246, 1953. 

Assistant Professor, Engineering Department. 
Mem. ARS. 

2 Numbers in parentheses indicate References at end of paper. 


Fesruary 1956 


of the gas and the unburned elements, two elements which 
burn at the same distance from the injector face will, in 
general, possess different injection velocities, follow different 
geometric paths, and therefore have different time lags. 
This assumption allows greater generality without making 
unnecessary, restrictive, yet uncertain description of the 
flow. Systems with distributed combustion should, of 
course, be treated separately. 

In previous analyses the sensitive time lag (dimensionless) is 
assumed to be of the order of unity so that all propellant ele- 
ments sense the variations of the thermodynamic states of the 
burned gas only in the immediate neighborhood of the combus- 
tion front. When the sensitive time lag of an element is very 
large, the element will sense the chamber oscillations during 
its travel through a definite spatial range § — £; immediately 
upstream of the combustion front. The effect of such very 
large sensitive time lag and the effect of the spread of such large 
sensitive time lag are also treated in the present paper. 
The author is grateful to Prof. L. Crocco for initiating this 
investigation. 


Formulation and General Discussion 


As assumed in previous investigations, the unburned pro- 
pellant elements suspended in, and carried along by, the 
burned gas stream are assumed to occupy negligible volume 
so that the entire cylindrical combustion chamber is filled 
with burned gas. The mean flow velocity of the burned gas 
upstream of the concentrated combustion front (region 1) 
is practically zero and that downstream (region 2) is a con- 
stant determined by the geometry of the deLaval nozzle. 
Let x represent the axial distance downstream of the injector 
end as a fraction of the combustion chamber length and the 
concentrated combustion front be at x = & For the acousti- 
cal oscillations in each of the two uniform flow fields, the time 
independent part of the velocity perturbation v(x) as a frac- 
tion of the stagnation sound speed in the burned gas and that 
of the density perturbation 6(x) as a fraction of the stagnation 
density are given (2) by Equations [1] and [2] 


= cy, exp (+ ax)[1 — exp 
= —c,, exp (+ar) [1 + exp (—2axr)]f 


v(r) = EXP (; {1 = ; »|t 
= —c,, exp {1 + Bew| 
1-@ 1 — # 


where a = \ + ww is the characteristic constant with \ repre- 
senting amplification coefficient and w the angular frequency 
based upon the wave propagation time for the acoustical wave 
to travel the length of the combustion chamber filled with 
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stagnant burned gas. @ is the mean gas velocity as a fraction 
of the sound speed in stagnant burned gas. The constant B 
is defined by 


1 + 1(6, 
I(6, 


with J standing for the specific acoustical admittance ratio of 
oscillating gas flow through a deLaval nozzle. _ The specific 
acoustical admittance ratio is the ratio of the fractional 
velocity variation to the fractional density variation of the 
gas entering the nozzle. This admittance ratio depends on 
the reduced frequency £ of the oscillation, the steady-state 
gas velocity entering the nozzle, and the geometry of the sub- 
sonic portion of the nozzle. The values of J(8, %) for a nozzle 
with linear steady-state velocity in subsonic part are repro- 
duced from (4) as Fig. 1. For this particular noazzle,. the 
reduced frequency £ is equal to the ratio of the angular fre- 
quency of oscillation and the steady-state velocity gradient 
in the nozzle. The energy outflow through the nozzle con- 
stitutes an inevitable damping action even in the absence of 
viscosity and conductivity. 


BB, a) = 


z= a f 


N 


4 


Fig. 1 Specific admittance ratio 1(8,z) = R + iS of super- 
critical nozzle as a function of reduced frequency 6 = w/z and 
dimensionless entering velocity 7 


The solutions [1] and [2] are to be matched together 
through the boundary conditions at the concentrated com- 
bustion front z = &. The first part of the boundary condition 
is the continuity of the density perturbations of burned gases 
across the combustion front at any instant. Accordingly, 
the characteristic constant a = \ + ww is continuous across 
the front. The second part of the boundary condition is con- 
cerned with the velocity discontinuity ». — » across the com- 
bustion front that is prescribed by the variation of the local 
burning rate due to the variation of the sensitive time lags of 
the propellant elements. It is through this boundary con- 
dition that the time lag distribution affects the stability of 
the system and that the present investigation differs from that 
of (2). 

The instantaneous value of the sensitive time lag 7 of a 
given propellant element is related to its steady-state value 
7 in terms of the local gas pressure p at the position x’ at the 
instant ¢’ by the following integral relation as introduced by 
Croeco (1) 


[x’(t’),t’ jdt’ = const = p" [x'(t’),t’] dt’. . [4] 


t 
Here the integrand represents the local instantaneous rate of 


activation and the constant value of the integral represents 
the total amount of activation that is required for the given 
element to start the instantaneous combustion. The ex- 
ponent n is called the interaction index which represents the 
average extent of interaction between the rates ef combustion 
processes and the different oscillations of the correlated phys- 
ical quantities such as pressure, temperature, etc., existing 
in the combustion chamber. 

From physical considerations of the successive steps of the 
processes leading to complete combustion in a liquid rocket, 
the sensitive time lags of different elements are expected to 
scatter in a certain region in the neighborhood of some mean 
value 7m = + fmin)/2. Let be the fractional total 
extent of the spread, 7max — fmin = 2eo7m. Let fie) be the 


- fractional amount of the propellant elements having sensi- 


tive time lags lying between 7 = 7m(1 + €) and fmin = 7m (1 — 


flo) =1, f(-#) =90 f ; 
—«, de 

Since we are considering only instability of the intrinsic 
type where the rate of injection and the pattern of distrilbu- 
tion of propellant elements are independent or insensitive to 
the oscillations in the combustion chamber, f(¢) is independ- 
ent of time. The rate at which propellant elements with 
sensitive time lags lying between 7 and 7 + d are being 
transformed into burned gas is given as 


The total rate of burned gas generation m, is obtained by 
integrating ™;dz7 over the complete range of 7. Differen- 
tiating Equation [2.4] to obtain dr/dt and replacing the 
instantaneous quantities by the sum of mean quantities and 
their respective perturbations of the exponential type, one 
obtains for local isentropic small oscillations the fractional 
burning rate perturbation 


5] 


Both df/de and + = 7m(1 + e¢) of the integrand are functions 
of «. If 4(€) is not zero, Equation [6] can be rewritten as 


= yne—™ [1 — [7] 
mi 
with 


where c and ;, are real quantities. 

Equating (m, — mi)/m; from Equation [7] to the frac- 
tional variation of the net mass outflow rate from the con- 
centrated combustion front, one obtains 


vo — + — yn)d + exp (—az.) = 0.... [9] 


which prescribes the discontinuity of the velocity perturba- 
tion ve Ye 

By substituting Equations [1] and [2] into Equation [9] 
and noting that the Mach number of gas flow is small (except 
for the case of throatless motor) so that 1 — u? & 1, one ob- 
tains the following characteristic equation 


— Bexp [2a(1 — 
a + Bexp [2a(1 — 


— tanh «| = 


(1 — yn) + yne exp (—az-,).. . [10] 


for the determination of the characteristic constant a = 
\ + ww of the system. Equation [10] is closely analogous 
to Equation [4.1.1] in (2). The only difference is that the 
expression exp (—a7) in Equation [4.1.1] is replaced by 
c exp (—a7-) in the present case. 
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For the determination of the stability boundary with a = 
iw, c represents the magnification factor of the effect of the 
time lag spread and ;, represents the effective time lag that 
enters into the stability calculation of the combustion system 
with the particular distribution function f(e) in question. 
Both c and 7, in general depend on the type of the distribu- 
tion function f(e), the extent of spread 27,,¢9, the position & of 
the concentrated combustion front, and especially the fre- 
quency w of the particular mode of neutral oscillation. Once 
the dependence of c and 7, on these parameters has been 
found for a particular distribution of time lag, the deter- 
mination of the stability boundary and the frequency of the 
neutral oscillation can be carried out. 

Since during the sensitive time lag, the propellant elements 
are upstream of the combustion front, that is, in region (1), the 
value of 6(&) and 6(€) should be evaluated using Equation 
{1]. In general, &; is different from £; then the expression 
[8] would encounter some singularity when 6(£) = 0 which 
does occur when the combustion front is located at the node 
of the pressure oscillation in the combustion chamber. This 
is only introduced by the division by 4(¢) which is not at all 
necessary except for the purpose of convenience and com- 
parison with the corresponding expressions in (2). 


Systems With Time Lag Comparable to Wave 
Propagation Time 


We shall first consider the case where the value of the pres- 
sure sensitive time lag = is of the order of unity as was treated 
in (2). Since the velocity of the unburned propellant ele- 
ments approaching the combustion front is small as com- 
pared to the small quantity u which is the velocity of the 
burned gas downstream of the front, it is possible to take 
ti(¢ — 7) = &(t) as was done in (2). Equation [8] for neutral 
oscillations then becomes simply 


d 
0 


The only complex quantity involved in the integral of Equa- 
tion [11] is exp [—tw7m(1 + €)] whose modulus is unity. 
According to Equation [5] the integral of the modulus of the 
complex integrand is unity. Therefore, for neutral oscilla- 
tions and for completely arbitrary distribution of 7, of the 
order of unity, we must have 


If (df/de) is a symmetric or even function of ¢, one has 


The even function of (df/de) corresponds to odd function 
of f(e); i.e., antisymmetric distribution of f(#) with respect 


to 7 a8 was assumed in (1) while dealing with the low-fre-. 


quency case, and 7,, was used in the place of 7, for such anti- 
symmetric f(7). The simplest examples of such symmetric 
(df/de) are 


@ 1 i 
(i) of = — = const with a = —— 
de WE0T m 
de 4a Be = 


as are given in (1). For arbitrary symmetric function 
(df/de), the integral*in Equation [13] can be carried out 
analytically by expanding df/de into Fourier series and inte- 
grating term by term. The coefficients of the Fourier ex- 
pansion should, of course, be properly normalized so that 


° Tf the integral is negative, c should be taken as the absolute 
value of the integral and 7, as r/w — 7m. 
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Fig. 2 Magnification factor C of the effect of the spread of the 
sensitive time lag as a function of the ratio of the extent of spread 
to the period of oscillation 27,,¢)/T 


Equation [5] is satisfied. Both c; and c2 are plotted in Fig. 2 
as damped oscillating function of w7neo/7. (The asymptotic 
behavior of c; is 1/w, and that of cz is 1/w*.) The value of c 
for other forms of df/de has a similar behavior. 

It is important to note that, for particular type of dis- 
tribution function, the value of ¢ depends not only on the 
manner of distribution f(e) but also on the ratio of the total 
extent of the time lag spread 27,,€9 to the period T of the given 
neutral oscillation. The magnitude of c remains less than 
unity and approaches unity when ¢) approaches zero, as it 
should be. Furthermore, the factor c could vanish for par- 
ticular values of the time lag spread as compared to period of 
oscillation for each type of distribution. In Equation [14] 
the value of c, for type (i) distribution vanishes when the total 
extent of time lag spread 27,9 is equal to any integral mui- 
tiples of the period; and c2 for type (ii) distribution vanishes 
when 27,¢9 is equal to any odd integral multiples of half 
period not less than three halves of the period. Since w is ap- 
proximately integral multiples of 7, and fmeo is only sub- 
jected to the restriction that 7me9 < 7m, the situation that 
c = Ois possible. 

The magnitude of c is of considerable importance in the sta- 
bility behavior of a given system as can be seen from the fol- 
lowing qualitative consideration. In the case of intrinsic 
instability, the only excitation agent is the variation of the 
specific rate of burned gas generation which is produced by 
the variation of the sensitive time lag. Equation [6] shows 
that this variation of gas generation rate can be divided into 
two parts, the first part is due to the pressure variation at the 
end of the time lag which is the same for all propellant ele- 
ments burning simultaneously at the same location. The 
second part is due to the pressure variation that each element 
sensed in the beginning of its sensitive time lag. The con- 
tribution of the second part is different for different elements 
having different time lags even if they burn together. With 
optimum timing, the sum of the exciting contribution of the 
second part is added to that of the first part in stimulating 
instability. The quantity c represents the maximum amount 
of exciting contribution with optimum timing of the second 
part that is available to augment the excitation of the first 
part. Therefore, intuitively, we see that smaller values of 
¢ means a more stable system with optimum timing. The 
fact that c is, in general, less than unity (Equation [12]) when 
there is a spread of time lag indicates that the effect of the 
spread of time lag is always stabilizing so far as the uncondi- 
tional stability is concerned. The stabilizing effect is maxi- 
mized when c = 0; that is, when the increase of burning rate 
due to the pressure variations that a group of elements sensed 
in the beginning of their sensitive time lags is just balanced 
by the decrease of the burning rate due to that of another 
group of propellant elements. The qualitative result that 
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a smaller value of c means a more stable system as obtained 
from the prévious physical arguments is clearly verified by 
the calculated results given in later sections of the present 
paper. 

Since c is a function of both w and 7, the determination of 
the stability boundary with given time lag spread involves 
the simultaneous solution of Equations [10] and [11]. How- 
ever, if 27meo/7’ <1 for the fundamental or the first few 
higher modes, the value of c for each mode of neutral oscilla- 
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Fig. 3 Critical values of the interaction index n corresponding to 
neutral oscillations of frequency w = Bx, with c = 0.9 when the 
combustion is concentrated at different position ¢, as a fraction 


of the chamber length 
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Fig. 5 Unstable ranges of the effective sensitive time lag 7, and 

the position ¢ of concentrated combustion front from the injector 

face as a fraction of the chamber length when interaction index 
n=1 


90 


Fig. 4 The lowest critical values (h = 0) of the effective sensi- 

tive time lag 7 corresponding to neutral oscillations of frequency 

w = Bx with c = 0.9 when the combustion is concentrated at dif- 
ferent position ¢ as a fraction of the chamber length 


tion [with kr — w = O(a) and k = a given integer] can be 
taken as approximately constant. Then the critical values 
of interaction index n and time lag 7, corresponding to neutral 
oscillations of frequency w are obtained as 


11-X l+e 
=- 1 


.. [15] 


with 


X+i¥ — tanh |... [16] 


— Bexp [2iw (1 — &)] 


1 + Bexp [2iw (1 — &)] 


For given rocket geometry X + iY is a function of w only. 
Typical results of calculation for the case with 8 = w/z are 
shown in Figs. 3 to 5. The minimum value of n compatible 
with any unstable oscillations in a given system is an index 
of the relative stability of systems with different geometry. 


Cc 
“O 2 4 6 8 10 


Fig. 6 Minimum values of n compatible with any unstable high 
frequency oscillations vs. magnification factor c of time lag spread 
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Figs. 7a, b Critical values of the interaction 1 > < AY 
index n vs. reduced frequency 6 = w/rwith 
cosine-type time lag distribution (a) «) = 1, — rege? = 
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Larger Mmin indicates a more stable system. For a given 
nozzle 


1 — X(w) _ 1 — X(wo) 
2y 
with w) determined by Y(w) = 0. The calculated results 
are shown in Fig. 6. The stabilizing effect of spreading the 
time lag is verified. 

For systems with 27,,¢)/7' comparable to unity, Equations 
[10] and [11] must be solved simultaneously and the function 
d/de must be specified. The solution can be effected 
through the relation 


(without time lag spread). . [17] 


sin wi, 


and the relation between c and w7,¢. Typical results are 
shown in Figs. 7 and 8 and unstable ranges of 7, and é are 
illustrated in Fig. 9, with n = 1 and c = cz. The values of 
min for the same system with two different types of time lag 
distribution are compared in Fig. 10. The system with con- 
stant df/de is found to be more stable than the same system 
with cos df/de. 


Systems With Time Lags Much Larger Than 
Wave Propagation Time 


When the sensitive time lag = is of the order of 1/% or larger, 
t; may be appreciably different from ¢. However, the extent 
of the time lag spread is not likely to be of the order of 1/%, 
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and é; can be considered as the same for all the elements hav- 
ing different time lags. Thus, in Equation [8], 5(&)/5(£) can 
be taken outside of the integral. Since 6(;)/4(€) is a real 
quantity, the effective sensitive time lag 7, for neutral oscilla- 
tions is not affected by the fact that §;< The only effect of 
&; < & appears as a multiplying correction factor 4(&;)/5(é) = 
cos wé;/cos wt on the magnification factor c as defined by 
Equations [11] and [13]; i.e., cs = c(cos wé;/cos wt). The 


nel 


4 


Fig.9 The lowest unstable ranges of the effective sensitive time 
lag 7, and the position ¢ of concentrated combustion front from 


the injector face as a fraction of the chamber length when n = 1 
with cosine-type time lag distribution 
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Fig. 10 Minimum values of n compatible with any unstable high 

frequency oscillations as functions of the fraction extent «o of 

time lag distribution; dotted curve, constant df/de; solid curve, 
cosine df/de 


characteristic Equation [10] remains unchanged and the sta- 
bility boundary can be determined in exactly the same manner 
as described in last section. By comparing corresponding 
curves for given values of & and ¢, in Figs. 11 and 12, the 
effect of time lag spread about the mean value is found to be 
stabilizing. 

Based on the previous physical argument, it can be ex- 
pected that the effect of &; < é is stabilizing or destabilizing 
as compared to the case with ¢; = & depending on whether 
|cos wé:/cos wt| |. Since both cos wé; and cos wt can 
vanish, the effect seems to be quite complicated. However, 
more careful considerations indicate that it simply means a 
progressive displacement of the nodal and the antinodal 
positions downstream with the injector end being always an 
antinode of pressure oscillations. This trend is illustrated in 
Figs. 11 and 12 for the fundamental mode of oscillation with 
the position é for infinite nmin corresponding to the pressure 
node. For higher modes, the figures become complex but the 
same trend persists. As a result of the downstream displace- 
ment of the nodal positions, the whole pattern of stable 
and unstable ranges of 7, ¢ for given n is displaced similarly 
and the magnitude of this displacement increases with in- 
creasing magnitude of ¢ — &; as illustrated in Fig. 13. 
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Fig. 11 Minimum values of n compatible with fundamental mode 
of high frequency oscillations without time lag spread but with 
nonvanishing sensitive space lags ¢ — &;; that is, the fractional 
axial distance traveled by the propellant elements during the 
sensitive time lag 
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Fig.12 Minimum values of n compatible with fundamental mode 

of high frequency oscillations with small time lag spread c = 0.9, 

but with nonvanishing sensitive space lags ¢ — &;; that is, the 

fractional axial distance traveled by the propellant elements 
during the sensitive time lag 
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Fig. 13 Critical positions ¢, separating the stable and the unstable 
ranges of the positions of the combustion front as a function of 
sensitive space lag ¢ — & 
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Thermal Stresses in Solid Propellant Grains 


RICHARD D. GECKLER' 


Aerojet-General Corporation, Azusa, Calif. 


A method is presented for computing the maximum 
thermal stress in solid propellant grains that can result 
from temperature cycling. The method is based upon a 
first approximation obtained from the analytical solution 
for thermal stresses assuming Young’s modulus to be inde- 
pendent of temperature. A more accurate solution is 
then found by numerical integration taking account of 
the variation of Young’s modulus with temperature. 


OLID PROPELLANT rocket motors are subjected to 
varying ambient temperatures in nearly all applications. 
The possiblity that thermal stresses may cause the propellant 
grains ‘to crack therefore makes the computation of the 
thermal stresses a matter of practical importance. 

The present ‘discussion is concerned with hollow, circular 
cylinders often employed in rockets as either unrestricted- 
burning or internal-burning grains. The cylinder of external 
radius R,, internal radius R;, and infinite length is initially at 
the temperature 7. At times ¢>0, heat is transferred to the 
external cylindrical surface from a medium at temperature 7’, 
with a constant heat transfer coefficient h; there is zero heat 
flux at the internal surface. The heat conduction equation is 


oT (27 Ri<r<R, 
with the boundary conditions 
or r=R, 
or t>0 [2] 
oT r=R, 
h(T, — T), [3] 


The thermal diffusivity x =  /pc, where ) is the thermal 
conductivity, p is the density, and c is the specific heat 
capacity of the substance, all assumed constant. The 
solution of these equations giving the temperature at any 
radial position p = r/R, depends upon (a) the rate of heat 
transfer to the outer cylindrical surface, which is specified by 
the Biot number @ = AR,/), (b) the time, which is specified by 
the Fourier number F = x«t/R,*, and (c) the ratio of the 
internal to the external diameters m = R,;/R,. Tables 
giving the temperature, in the form of the unconsumed 
temperature difference 6 = (7, — T)/(T, — T;), as a func- 
tion of these variables have already been presented in (1).? 

During the heating or cooling period the difference between 
the thermal expansion of the material near the internal and 
external surfaces will produce thermal stresses. The radial 
stress o, is zero at both surfaces, but it is a maximum near 
the middle of the web—a tensile stress during heating and a 
compressive stress during cooling. The tangential stress o» 
is zero near the middle of the web, but it is high at both sur- 
faces—a tensile stress at the internal surface and a compressive 
stress at the external surface during heating; a compressive 
stress at the inner surface and a tensile stress at the external 
surface during cooling. 

Since the initial temperature of the propellant grain is uni- 
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Thermo! Stresses ino Hollow Cylinder 
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Fig. 1 Profiles of tangential thermal stress in a hollow cylinder 
(E independent of temperature) 


form, the thermal stresses start out at zero; they then increase 
to a maximum absolute value with time, and decrease to zero 
once more as the temperature again becomes uniform. This 
behavior is illustrated in Fig. 1 showing tangential stress pro- 
files through the web at various times for a cylinder with m = 


0.20. In this figure we use a dimensionless measure of the 
tangential stress that is defined by 
(1 — 


where » is Poisson’s ratio, E is the modulus of elasticity, and a 
is the linear coefficient of thermal expansion. A similar 
dimensionless measure of the radial stress is defined by 


(1 — v)o, 
Ea(T, — T;) 


In both cases positive values of the stresses refer to tension 
and negative values refer to compression. 

In the simplest case where all the thermal and mechanical 
properties of the cylinder are independent of temperature, 
the radial and tangential stresses are given by Timoshenko (2) 


as 
1 — R? Re 


1 +- R? 


The longitudinal stress ¢, depends upon the nature of the con- 
straints applied at the ends of the grain; in the only case we 
shall consider, the resultant of the normal forces over any 
cross section of the cylinder is zero and the longitudinal stress 
is given by 


2, 


The analytical solutions to Equations [1], [6], and [7] were 
given by Jaeger in (3) for somewhat more general boundary 
conditions than we need here. When simplified to relate to 
the boundary conditions expressed by Equations [2] and [3], 
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Jaeger’s results can be written as follows 


= exp (—55;*)fi(m, 8, 5;) 


[ — 5:) — fa poi, may | [10] 


Zo = — exp @, 


$=1 


1 2 


Here, the functions f,, fe, fs, and f, are defined by 


— BJo(6;)] 
— BJo(8;)]* — (@* + 


pdi, mb;) = — Yo( pd: )Ji(ms;) .. . [13] 

Fd = /pd;. . . [15] 
and the 6; are roots of equation 

— BYo(5)] — — = 0. . [16] 


In these equations, Jo and J; are Bessel functions of the first 
kind, and Yo and Y; are Bessel functions of the second kind 
according to Watson’s notation (4). 

Since the roots of Equation [16] had been previously com- 
puted (5) in order to obtain the temperature distribution 6, 
relatively little added labor was needed to also obtain values 
of Zgand 2,. Actually for design purposes it is only necessary 
to consider the tangential stress, since its maximum absolute 
value at both the internal and external surfaces is greater 
than the maximum absolute value of the radial stress which 
exists somewhere near the middle of the web. Furthermore, 
it is unnecessary to know the tangential stress as a function of 
time. In Fig. 1, for example, one is interested in knowing 
only the value of 2» at the peak for p = mor at the valley for 
p = 1.0, depending upon whether the grain is being heated or 
cooled. The results of computations based upon Equation 
[11] are given in Table 1 for the outer surface p = 1.0 and in 
Table 2 for the inner surface p = m. Both the maximum 
value of 2», and the value of the dimensionless time ¥ at which 
it occurs are given. The data for the inner surface are also 
given graphically in Figs. 2 and 3. 

As an example of the use of these tables we may consider 
the following problem: A propellant “grain’’ with an internal 


film, 8, [12] 
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Fig. 2 Maximum tangential thermal stresses at the inner surface 


of a hollow cylinder vs. Biot number (# independent of tempera- 
ture) 
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Fig. 3 Value of Fourier number when tangential thermal stress 


at inner surface of a hollow cylinder is a maximum vs. Biot 
number (# independent of temperature) 


diameter of 2.78 in. and an external diameter of 9.71 in., 
initially at a temperature of —30°F, is suddenly placed in a 
room maintained at 150°F. The metal case and air space 
between it and the propellant, together with the resistance of 
the air film on the outside of the case, result in an over-all con- 
vection coefficient of 3.74 Btu-ft~?-hr-!-°F-!. The thermal 
ductivity of the propellant is 0.181 Btu-ft-?-hr-'-(°F/ 
ft)“, its thermal diffusivity is 0.00578 and _ its 


Table 1 Maximum tangential thermal stresses at the outer surface of a hollow cylinder and the times at which they occur 
(EZ independent of temperature) 
——0.2 —0.8 
— Ze —Ze —Ze —Ze F 
0.01 0.00231 0.35 0.00186 0.22 0.00131 0.11 ai 0.033 
0.05 0.0110 0.25 0.00911 0.17 0.00638 0.086 0.00340 0.026 
0.10 0.0216 0.23 0.0178 0.15 0.0125 0.075 0.00670 0.023 
0.20 0.0408 0.18 0.0340 0.12 0.0243 0.066 0.0130 0.020 
0.50 0.0890 0.13 0.0764 0.094 0.0563 0.052 0.0306 0.016 
1.0 0.151 0.11 0.133 0.072 0.101 0.041 0.0576 0.013 
2.0 0.239 0.068 0.215 0.052 0.171 0.031 0.104 0.011 
5.0 0.383 0.037 0.355 0.029 0.302 0.019 0.205 0.0074 
10 0.498 0.021 0.472 0.017 0.419 0.012 0.310 0.0053 
20 0.596 0.011 0.569 0.0093 0.523 0.0070 0.432 0.0037 
50 0.728 0.0039 0.705 0.0036 0.665 0.0031 0.605 0.0022 
100 0.816 0.0017 0.795 0.0016 0.765 0.0016 0.713 0.0015 
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Table 2 Maximum tangential thermal stresses at the inner surface of a hollow cylinder and the times at which they occur 
(EZ independent of temperature) 
0.2 0.4 0.6 0.8 
8 ze F Ze F ze F Ze F 
0.01 0.00199 0.35 0.00136 0.22 0.00081 0.12 0.00038 0.040 
0.05 0.0096 0.30 0.00663 0.20 0.00402 0.10 0.00183 0.032 
0.10 0.0187 0.27 0.0130 0.18 0.00792 0.087 0.00363 0.028 
0.20 0.0352 0.24 0.0248 0.16 0.0153 0.077 0.00711 0.024 
0.50 0.0766 0.18 0.0554 0.13 0.0353 0.063 0.0170 0.019 
1.0 0.127 0.15 0.0953 0.10 0.0632 0.052 0.0320 0.016 
2.0 0.196 0.12 0.151 0.080 0.106 0.043 0.0574 0.013 
5.0 0.291 0.091 0.238 0.061 0.181 0.033 0.112 0.0099 
10 0.349 0.078 0.297 0.051 0.240 0.028 0.167 0.0087 
20 0.382 0.070 0.336 0.046 0.285 0.024 0.221 0.0077 
50 0.413 0.064 0.369 0.040 0.327 0.020 0.289 0.0070 
100 0.426 0.060 0.380 0.038 0.342 0.019 0.324 0.0067 


linear coefficient of thermal expansion is 5.2 X 1075 °F-1. 
What is the maximum tensile stress in the propellant and at 
what time does it occur? 

Now R, = 0.405 ft; hence the dimensionless parameters 
are 


® = 3.74 X 0.405/0.181 = 8.37 
m = 2.78/9.71 = 0.286 


Since the propellant grain is being heated, the maximum 
tensile stress occurs at the inner surface. Hence Table 2 is 
used, and simple linear interpolation gives (a) for m = 0.2: 
2, = 0.330, § = 0.082; and (b) for m = 0.4: 2 = 0.278, 
§ = 0.054. Linear interpolation between these results gives 
finally for m = 0.286: 2, = 0.308, § = 0.070. The space- 
mean temperature of the propellant grain at this time is 
approximately 40°F, and the modulus of elasticity of the 
propellant will be taken as approximately 5.6 < 10‘ Ibf.in~*. 
Then taking Poisson’s ratio as 0.35, we obtain the tangential 
stress as 


= > — T;)/(1 — v) = (0.308 X 5.6 X 104 X 5.2 X 
X 180)/(1 — 0.35) 
= 249 lbf-in.-? 


and the time at which it occurs is 


t = R25/« = 0.405? X 0.070/0.00578 
= 2.0 hr 


If the grain were being cooled from 150°F to —30°F in- 
stead of being heated, the calculations would be the same 
except that Table 1 is used instead of Table 2. For @® = 8.37 
and m = 0.286, we now find 2, = —0.446 and ¥ = 0.023. 
The space-mean temperature is approximately 120°F and 
the modulus of elasticity of the propellant at this temperature 


can be taken to be 3.3 X 10‘ lbf-in~*. The tangential stress 
is thus 
og = (0.446 X 3.3 X 104 X 5.2 K 10 X 180)/(1 — 0.35) 


= 209 Ibf-in.~? 
and the time at which it occurs is 


t = 0.405? X 0.023/0.00578 
= 0.65 hr 


It is obvious that numerical values of the tensile stress 
obtained as above are not very accurate. Tables 1 and 2 were 
computed on the basis that all the material properties and the 
convection coefficient were independent of temperature. In 
many types of heat conduction problems these assumptions 
are adequate for engineering purposes. However, solid 
propellants are plastic materials and their moduli of elasticity 
may change by as much as a factor of ten over the tempera- 
ture range of interest. Consequently, the assumption that 
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the modulus of elasticity is independent of temperature is not 
a suitable assumption on which to base quantitative calcula- 
tions of thermal stresses. 

Nevertheless, despite the fact that the foregoing results 
must be viewed as crude approximations, they adequately 
demonstrate one important conclusion: the maximum stresses 
occur relatively early in the heating or cooling period while 
the cylinder is near its original temperature; for plastic mate- 
rials the much greater modulus of elasticity at low tempera- 
tures then means that higher stresses will occur during heat- 
ing from a low initial temperature than during cooling from a 
high initial temperature; since failure will occur in tension 
rather than in compression, it is clear that cracks as a result 
of temperature cycling should start most frequently at the 
inner surface during the heating portion of a temperature 
cycle. 

The idealized thermal stresses given in Tables 1 and 2 also 
permit one to estimate in a semiquantitative fashion the 
relative effects of changes in the physical parameters. For 
example, in the heating problem given above, if the effective 
convection coefficient h were reduced by external insulation 
to one tenth its former value, thus reducing @ by the same 
factor to 0.837, then the maximum value of 2» at the inner 
surface would be reduced from 0.308 to 0.094; i.e., by a factor 
of three. A similar reduction would result if the thermal 
conductivity could be increased tenfold or if the external 
radius of the grain could be reduced by the same factor. 

In order to obtain a closer estimate of the thermal stresses 
by taking into account the variation of Young’s modulus with 
temperature, we may make use of Hilton’s solution of the 
equations of equilibrium and compatibility (6). On the 
assumption that the medium is incompressible, i.e., »v = 0.5, 
Hilton derives equations that may be put in the following 
form for the radial and tangential stresses as a function of the 
radial coordinate. 


E 
= ofr) + - ran | [17] 
_ PAR.) 
= FAR.) F\(r) — Fir) [18] 
where 

Ri 

K(r) 


dT 
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= 2[a(r)r*T(r) — a( 
a(r)rT(r)dr — 2 f r*T(r)de(r) 
Ri 


If a is a constant independent of temperature, the last term 
in Equation [21] is zero. 

A complete calculation based on the foregoing equations 
for the example given above is shown in Table 3, and the 
tangential stress distribution is shown in Fig.4. It is assumed 
in this computation that a is independent of temperature, 
thereby permitting the use of only the first three terms on the 
right side of Equation [21]. The temperature dependence of 
Young’s modulus was assumed to be expressed by the em- 
pirical equation ZH = 7.3 X 10‘ exp [—0.0065T). Since Z 
is not constant, the time at which the maximum thermal stress 
occurs is not necessarily that given by Table 2. Trial and 
error computations in the present example show that the 
time is no longer ¢ = 0.070 (equivalent to ¢ = 2.0 hr); it 
has become approximately ¢ = 0.05 (equivalent tot = 1.4 hr). 
The transient temperature distribution at the time computed 
as described in (1) was used as a basis for the computations 
in Table 3. 

Although the absolute values of the stresses computed in 
Table 3 do not possess great accuracy because of the assump- 
tion vy = 0.5, they are adequate for engineering purposes 
to compare with other values computed in the same way. 
The relative effects of changes in the major design parameters 
h, and R, and R; can be defined quite precisely by calculations 
based on Equations [17] to [21]. Examination of the detailed 
calculations in Table 3 reveals that the distribution of 
thermal stresses can be markedly influenced by the manner 
in which Young’s modulus depends upon temperature. For 
this reason the results given in Tables 1 and 2 for the ideal 
case are not usually dependable for precise comparisons and 
are most useful merely to obtain a preliminary estimate at the 
time near which the maximum stress occurs. It is to be re- 


tangential stress og - Ibf - in-~? 
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Fig. 4 Tangential thermal stress vs. radial position in a hollow 
cylinder for § = 0.050, ® = 8.37, m = 0.286, E = 7.3 X 10‘ exp 
( —0.00657') 


marked, of course, that if the propellant is viscoelastic and 
exhibits marked viscous behavior, then it could happen 
that deductions based on Equations [17] to [21] would also be 
inerror. Actually, it appears that for many propellants their 
relaxation times are considerably greater than the time re- 
quired to reach the maximum thermal stress given in Tables 
1 and 2, and consequently their viscoelastic behavior can be 


neglected for engineering design purposes. 


The author desires to acknowledge the aid given him by 
Mr. Myron Lipow in performing many of the tedious cal- 
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culations required for the computation of Tables 1 and 2. 


Table 3 Numerical computation of thermal stresses at § = 0.050, ® = 8.37 with H = 7.3 < 10‘ exp (—0.00657') 
S4arT 2er*T — F; = (2) 
r, rz, r, 7; a, 4arT, (1), 2ar*T, 2ak?T (2), — (1), 
in. in.? °F oF-1 in. in.? in.? in.? in.? 
0.286 1.390 1.932 2.686* —20.4 5.2 X 10-§ —0.00590 Or —0.00410 0 0 
0.3 1.456 2.120 3.087 -—20.4 5.2 X -—0.00618 -—0.00040 -—0.00450 —0.00040 0 
0.4 1.942 3.771 7.324 -—18.3 5.2 X -—0.00739 —0.00369  -—0.00718 —0.00308 0.00061 
0.5 2.428 5.895 14.31 -12.4 5.2 xX 10-* -—0.00626 —-0.00701 -—0.00760  -—0.00350 0.000351 
0.6 2.913 8.486 24.72 -1.0 5.2 X 10-* —0.00061  —0.00868 -—0.00088 0.00322 0.0119 
0.7 3.398 11.55 39.23 16.8 5.2 X 107° 0.0119 —0.00594 0.0202 0.0243 0.0302 
0.8 3.884 15.08 58.59 41.0 5.2 X 10 0.0331 0.00499 0.0643 0.0684 0.0634 
0.9 4.370 19.10 83.45 71.4 5.2 X 10° 0.0649 0.0288 0.142 0.146 0.117 
1.0 4.885 23.57 114.4 106.2 5.2 X 107 0.1072 0.0705 0.260 0.264 0.194 
Table 3 (continued) 
E E 
= E\QXFi, o = E_ 6 = 
E, Ibf- r3}’  Ibf- Ibf- (6), Ibf - = (8), + (8), 
0.286 8.34 X 10‘ 31,000 0 0 0 0 0 0.00645 43,100 278 278 
0.3 8.34 X 10‘ 27,000 1,910 0 0 12.3 12.3 0.00645 39,300 254 266 
0.4 8.21 X 10‘ 11,200 11,360 6.83 0.00584 21,800 127 199 
0.5 7.91 X 10‘ 5,530 15,500 19.4 8.02 100.0 92.0 0.00294 13,400 39.5 131 
0.6 7.33 X 10 2,960 17,600 35.2 21.3 113.5 92.2 —0.0054 8,630 —46.6 45.6 
0.7 6.55 X 10 1,670 18,800 50.5 42.1 121 79 —0.0237 5,670 —134 —55 
0.8 5.59 X 10 953 19,400 60.4 69.0 125 56 —0.0569 3,710 | —211 —155 
0.9 4.59 X 104 550 19,800 64.4 99.2 128 29 —0.111 2,400 —267 —238 
1.0 3.66 X 104 320 20,000(3) 62.1 129(4) 129 0 —0.188 1,550 —292 —292 
Nore: ‘S means integration from R; to r by the trapezoidal rule: Sf(r)dr = Ar[}(fo + fn) +f + fe + + fr-1]Q = 
F(Re)/Fi(Re) = (4)/(3) = 0.00645. 
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A Comparative Investigation of a Homogeneous Combustion 
Chamber With a Two-Stage Combustion Chamber 
(Continued from page 80) 


It is now necessary to solve Equation [3] for J. On Fig. 5 
the solution of Equation [3] is shown in graphic form where 
T: = 2019. This solution corresponds to the type repre- 
sented on Fig. 4 by point @. 


2x10" 


CURVE 421600 


2019 


2000 2030 
T, in°R 


Fig. 5 Graphic solution to Equation [3] 


The heat release o can now be evaluated from Equation 
[XI] as 


Vi + V2 48.8\600 T: 
= 29,900 
(: be Ts 
Y2 


where it should be noted that the bracketted portion from 
Equation [1] is equal 0.944 X 10~, and since V; = V2 


1\ 47.8 
(3) 1) X 0.944 X 10 

= 10.9-:107 

The combustion efficiency as defined by Equation [XII] 
can be computed as 

7.34: 48.8 

«= 878,000 (2019 — 600) = 0.58 

This represents the circled point on Fig. 3 for the curve 
with the first stage operating stoichiometrically. Other 
points on this curve were obtained with arbitrarily chosen 
values of 7; > 3310°R using essentially the same procedure. 
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Technical Notes 


Transient Response 
of the Turbine Flowmeter 


JERRY GREY! 


The Daniel and Florence Guggenheim Jet Propulsion 
Center, Princeton, N. J.. 


N recent years the turbine-type flowmeter has been used 
almost exclusively for steady-state flow measurement in 
rocket and jet engine applications. With the recent growth of 
interest in the measurement of instantaneous values of os- 
cillating flow rates, an examination of the transient response 
of the turbine-type flow sensing element is of some signifi- 
cance. 

Consider a flowmeter which consists of a rotor mounted on 
bearings inside a pipe. The rotor is equipped with N blades, 
each of which is a thin flat plate of tip radius rz, root radius 
r,, and chord c. The blades are inclined at an angle a to the 
axis of the pipe. 

The transient response of this type of sensing element is 
best described by the time constant of the rotor when sub- 
jected to a step change in fluid velocity. That is, if an in- 
stantaneous increase occurs in the flow rate (see Fig. 1), the 
rotor will accelerate from a speed corresponding to the original 
flow rate to a speed corresponding to the new flow rate. The 
time required to accelerate the rotor to its new speed (or 
some fraction of its new speed) is a measure of the time con- 
stant. 


Flow rote 


Rotor speed 


Time 


Fig. 1 Response of sensing element rotor to step increase in 
flow rate 


For the purposes of this analysis, rotor bearing friction 
(which is very small in certain commercial meters of this type) 
will be neglected, and the fluid will be assumed to be an in- 
compressible, nonviscous liquid: The equation of the rotor 
may then be written 


fel. 


where 


T = hydrodynamic torque on the rotor 

I = moment of inertia of the rotor about its axis 

6 = angular displacement 

i = time 

To compute the fluid torque T' on the rotor, it is necessary 
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to find for each blade that component of the hydrodynamic 
foree which produces torque (see Fig. 2). When the fluid 
velocity u is constant, the velocity d@/dt of the blade element 
is such that the effective angle of attack ao is zero, and hence 
the lift force L’ per unit blade length is zero. (Actually, L’ 
and ay are just large enough to balance bearing friction and 
viscous drag effects.) When the fluid velocity wu is increased, 
ay (and hence L’) becomes finite, and the rotor begins to speed 
up. When the blade element velocity r(d@/dt) becomes 


Fig. 2 Hydrodynamic forces on a blade element subjected to a 
step increase in flow rate 


sufficiently large so that a» is again zero, the rotor has reached 
equilibrium, and runs at a steady speed corresponding to the 
new fluid velocity. 

During the transient speeding-up period, the torque ex- 
erted on a given blade element of chord c and length dr is 
(referring again to Fig. 2) 

aT = Lrdr 
= L’ cos é rdr 
where L’ is the force per unit blade length perpendicular to the 
effective velocity V 
or 
aT = Cy’ cos 6 (1/2 pV%c)rdr 


where = lift coefficient = L’(1/.pV*c) 


For a blade of finite length, the effect of induced ‘drag must 
be introduced (1).? 
3 SiN a 
+ (2n/M) 
The torque on a blade element thus becomes 


sin 1 ue \? 
dT = 1 + (r/R) cos » rdr 


CL 


= blade “airfoil’’ efficiency ~0.9 
ao = effective angle of attack = a — 6 
= blade angle with pipe axis 


tan“! r U2 
dt 


fluid density 

final flow velocity (after step) 

initial flow velocity (before step) 

chord of blade element 

radius of blade element (measured from rotor axis) 
(rz — 11)? 


blade area 


= 


blade aspect ratio = 


Intregrating along one blade 
cr sin ap 


* Numbers in parentheses indicate References at end of paper. 


Eprror’s Note: This section of Jer Propusion is open to short manuscripts describing new developments or offering comments on 
— previously published. Such manuscripts are published without editorial review, usually within two months of the date of receipt. 
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Since the blades are generally short and trapezoidal in 
shape, only a small error will be involved if we take a mean 
ll radius R and a mean chord C. With this simplification, and 
taking account of the number of blades N, the equation of 
motion of the rotor becomes F} 
N sin ap CR(r2 — 171) 
nce Making the proper trigonometric substitutions, A 4 
, 
and rnpus*NCR (rz — 11) R do E MY 
eed 
Now, let 
mes 
— 11) 
Fig. 4 Time constant of 1 in. sensing element as a function of 
and flow rate for step change in flow rate 
7 
Ue d t 4 
then 
Oa dz Rd* max. flow rote 
dt wu,dt? B5 
and the equation of motion becomes 
d 
ex- A cosa (tana —s)=— 
x dt 
is Max. roted flow 
Integrating to find the time 7 required for the rotor to 3+ een 
accelerate from a speed corresponding to the initial fluid 
velocity to that of the final velocity ue 
he dz Blade angle (degrees) 
0 Acosadw, tana —z Fig. 5 Time constant of 1 in. sensing element as a function of 
ua blade angle for step increase in flow rate 
It is obvious that the integration cannot be carried out to 14 
its upper limit, since it would, of course, take mathematically ; 
infinite time to reach this limit. In order to avoid this I : 
ist difficulty, it is customary to compute the time required to 
reach a fraction [1 — (1/e)] of the imposed velocity increment, ‘ 
where ¢ is the base of natural logarithms. This required time “a? : 
is called the time constant. 
The upper limit may thus be rewritten z max| flow rate = 
/ Max. rated flow 
Pick-up Coil Alnico Magnet ) 6 7 
Fig. 6 Time constant of Potter sensing elements with 40° blade 
angle for step increase in flow rate 
The time constant K may thus be obtained by 
1 
K 1 [1 = AG a de 
0 A cosa du, tan a —2z 
— tana 
Integrating and introducing the limits 
ae 1 
r. Acosa 
. In order to present numerical results, a commercial meter 
i Fig. 3 Potter electronic flow sensing element manufactured by the Potter’ Aeronautical Company (2) 
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has been examined. A typical sensing element is shown in 
Fig. 3. The results of the calculations are shown in Figs. 4, 
5, and 6. 

Fig. 4 illustrates values of the time constant K as a function 
of water flow rate for the standard one-inch Potter sensing 
element; Fig. 5 presents the variation of K with blade angle; 
and Fig. 6 shows the time constant of the different sizes of 
Potter elements. It is interesting to observe that the size of 
the assumed step increment in velocity has no effect on the 
time constant (i.e., on the time required to reach a fixed per- 
centage of the corresponding rotor speed increment), which is 
in accord with expectations for a meter with linear response to 
flow rate. 

It is especially interesting to note that the '/:-inch Potter 
sensing element has a time constant of the order of 1 milli- 
second at rated flow, and hence is able to respond accurately 
to variations in flow rate up to at least several hundred cycles 
per second. 

‘References 
1 Millikan, C. B., “Aerodynamics of the Airplane,” John 


Wiley & Sons, New York, 1948. 
2 Bulletin §-1, Potter Aeronautical Company, Union, N. J. 


A Missile Air/Gas 
Pressurization System 
GEORGE C. GENTRY! 

Air Force Missile Test Center, Patrick AFB, Fla. 


Developmental flight testing of missiles is in most cases 
the next step from laboratory tests and, as such; is subject 
to a great many variables. In this respect, the supporting 
systems must provide a degree of control approaching 
that found in -the laboratory. The air/gas system 
described here is one of several fluid pressurization systems 
in use at Cape Canaveral, Fla. This system was initially 
designed for use with air, but has since been adapted to 
handle some other gases as well, in supporting certain 
missile programs assigned to this Center. 


Introduction 


HE development of a tactical weapons system is a rela- 

tively complex and arduous process. The elements of 
any given program may be assigned to various agencies 
about the country in order to take advantage of the latest 
advances in the state of the art. Only the air/gas servicing 
system of the flight test phase of a missile program will be 
considered here. This system is designed to support tests 
during the research and development phase only. 


Missile Pressurization 


A missile system may require servicing with high-pressure 
air and/or inert gas in order to support the functions of its 
air-frame system, its guidance and control system, and its pro- 
pulsion system. These functions may vary with respective 
missile programs from pressurized air-frame sections, to null 
control systems or air bearings, to propellants pressurization 
or inert gas blanket, and others. Servicing requirements may 
range from 3000-psi to 6000-psi pressure at controlled flow 
rates, dew point, and capacity for air or gas or both. This 
ground support system may frequently become more com- 
plex than the missile system it is designed to support. ~ One 
type of system is shown in Fig. 1; this system has been de- 
signed to provide the maximum of support system flexibility 


Received November 1, 1955. 
1 Mechanical Engineer. _Mem. ARS. 


100 


/GAS STORAGE 
J 


GAS TRANSFER 


HYDRAULIC 
PUMP 


Fig. 1 Schematic diagram of missile air/gas pressurization 
system 


for some missile test programs at the AFMTC. The major 
components of this system may be described as a compressed 
gas storage and transport trailer, a mobile air/gas compressor, 
and a mobile boost transfer machine. 


Compressed Gas Storage and Transport Trailer 


This gas trailer, Fig. 2, is made up of a tube bank of 30 
tubes mounted on a standard military vehicle. The capacity 
of this unit is between 35,000 and 40,000 cu ft STP (standard 
temperature and pressure) at a working pressure of about 
2400 psi at 75°F. The manifolding at the front of the trailer 
features a rupture disk in each tube with individual gas escape 
lines leading upward and away from operating personnel. 
The manifolding control is located at the rear of the trailer 
in an enclosed cabinet which features temperature and pres- 
sure gages, individual tube shutoff valving, and a main outlet 
line and shutoff valve. This unit has proven quite satis- 


Fig. 2 Compressed gas storage and transport trailer 


factory as a gas supply source in missile operations. 


Mobile Air/Gas Compressor 


This compressor, Fig. 3, is a modified Hardie-Tynes unit, 
mounted in a standard military vehicle. The compressor is 
rated at about 85 SCFM (standard cu ft/min) and has a maxi- 
mum discharge pressure of 3500 psi. Certain alterations 
have been made in the form of additional piping and valving 
to allow the introduction of compressed gas into the com- 
pressor at its rated inlet pressure. To follow the compressor 
system portion of the schematic diagram in Fig. 1, we find 
that compressed gas may enter the system and by-pass 
the compressor entirely and yet pass through the driers and on 
to the boost transfer machine receiver vessels where pressure 
equalization will occur as follows: 

Initial pressure of gas storage trailer, 2400 psi; volume (by 
water measure) of gas storage trailer, 214 cu ft. 

Volume (by water measure) of receivers in boost transfer 
machine, total 37 cu ft. 
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Fig. 3 Mobile air/gas compressor 


Total volume of gas storage trailer and boost transfer 
machine, 251 cu ft. 


P,V, = P2V2 after pressure equalization at constant 


temperature 
2400 (214) = Pz» (215) 
P, = 2046 psi 


With the pressure in the boost transfer machine at 2046 
psi and the desired pressure at 3500 psi, it will be necessary 
to utilize the air/gas compressor to build up this remaining 
pressure. In order to operate this compressor with gas, 
two pressure reducing valves have been installed ahead of the 
compressor inlet. In actual practice, a few pounds of pres- 
sure above atmospheric (about 20 psia) is desirable to prevent 
the diaphragms of these pressure reducing valves from chatter- 
ing and eventually rupturing. Temporary use is being made 
of a Rawinsonde weather balloon as a surge chamber until 
an aluminum vessel can be fabricated. It may be noted 
here that to pump air with this compressor, the weather 
balloon is removed and the gas storage trailer valved off. 

In proceeding with the pressurization of the boost trans- 
fer machine with the air/gas compressor, it is necessary to 
purge the system and take a dew-point reading before intro- 
ducing gas into the boost transfer machine. The dew point 
is normally held to about —80°F. When a need for reactiva- 
tion of the driers is indicated, a reverse flow of heated nitrogen 
gas is introduced which proceeds through the driers and com- 
pressor by-pass line and out through the disconnected gas 
trailer line. 

The pump time to raise the pressure in the boost transfer 
machine from 2046 psi to 3500 psi is found to be about two 
hours by actual practice for helium gas. 


Mobile Boost Transfer Machine 


The boost transfer machine? is a rather complex piece of 
equipment, Fig. 4, and is utilized as a major component of 
this system for the purpose of storing compressed air or gas 
and delivering that air or gas upon demand at precisely 
controlled rates of flow, capacity, and pressure. 

In operation, the electrical system is energized, the hy- 
draulic system is set to by-pass flow back into the reservoir. 


2 See Reference at end of paper. 
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With the air receiver precharged to-3500 psi, the gas pressure 

is equalized between the receiver vessels and the accumulator 

transfer barriers. Pressure equalization is found to be: 
Volume of receiver vessels, 37 cu ft; pressure of receiver 

vessels, 3500 psi. 

Volume of accumulator transfer barriers total, 10.6 cu ft. 

Total volume of receivers and transfer barriers, 47.6 cu ft. 


3500 (37) = Ps (47.6) 
P, = 2846 psi 


The accumulator transfer barriers (accumulator battery 
no. | and no. 2, Fig. 1) are then isolated from the receiver vessels 
and hydraulic pressure is applied to the accumulator bladders 
to provide transfer of gas into the given missile system at the 
desired pressure. The hydraulic pump is of variable dis- 
placement to provide some control over gas delivery time. 
The comparator downstream of the accumulator transfer 
barriers serves as a double check on the dew point and will 
serve to detect possible bladder rupture or some other oil 
leak. The entire system features interlocks and fail-safe 
design. 


Fig. 4 Mobile boost transfer machine 


This typical system may be varied to suit individual missile 
requirements by valving off one of the receiver vessels or by 
adding receiver vessel volume by plumbing in additional 
receivers, or by pressurizing the accumulator transfer barriers 
initially with the receiver vessels, thereby eliminating the 
second step of pressure equalization. Also, one bank of 
accumulator transfer barriers may be isolated or either or 
both banks may be cycled to build up pressures as high as 
6000 psi. 


Operational Support 


In actual support operation, the procedure will be (for a gas 
charging operation) to connect the gas supply trailer line in 
and with an initial supply pressure of 2400 psi to utilize the 
air/gas compressor system to prime the boost transfer machine 
to its storage pressure of 3500 psi. This operation is normally 
performed a day or so in advance of the operational demand. 

The actual missile pressurization is performed within a few 
minutes time as a joint effort of the st transfer machine 
operator and the missile operator. The boost transfer ma- 
chine operator controls the entire marhine from the remote 
master panel in the Blockhouse. ‘The delivered air/gas 
pressure, flow rate, and capacity are monitored and controlled 
by the missile operator from the remote submaster panel in 
the Blockhouse. 
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Times for Interplanetary Trips 
ROBERT T. JONES! 
Ames Laboratory, NACA, Palo Alto, Calif. 


HE writer found it interesting to calculate the times re- 
quired to travel to the various planets at an acceleration 
of one “g,” or 8 X 10‘mph/h. ° 

Of course, with a thrust of such comfortable magnitude one 
can neglect surrounding gravitational fields except for a 
relatively short distance near take-off or landing. The 
“orbit” consists of an essentially straight line with the thrust 
directed toward the destination up to the halfway point, but 
in the opposite direction for the remainder so that the velocity 
is zero on arrival. 

The following table lists the approximate times required, 
and also the maximum velocities acquired in light units v/c, 


in days 
Mercury 2.25 0.0032 
Venus 1.51 0.0021 
Mars 2.08 0.0029 
Jupiter 5.88 0.0083 
Saturn 8.37 0.0118 
Uranus 12.07 0.017 
Neptune 15.45 0.022 
Pluto 17.80 0.025 


In spite of the small values of v/c the energy expended is 
certainly large by present standards. Nevertheless, it is 
interesting to see how quickly even “astronomical” distances 
succumb to the application of a thrust force of such a reasona- 
ble magnitude. 


Received Nov. 14, 1955. 
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A Quasi-Morphological Approach to the 

Geometry of Charges for Solid Pro- 

pellant Rockets; The Family Tree of 
Charge Designs 


J. M. VOGEL! 
The Ramo-Wooldridge Corporation, Los Angeles, Calif. 


One of the principal areas in the design of any solid 
propellant rocket is the selection of a suitable geometry for 
the propellant charge. A morphological classification of 
charge designs is developed, based upon their several dis- 
tinguishing characteristics. From the classification so 
developed it appears that at least 6240 charge designs are 
available, many of which are, of course, impractical for 
various reasons. Using the morphological classification 
as a datum, a quasi-morphological approach is discussed, 
which seems to be peculiarly well adapted to the particular 
field of charge design. It is shown that a family tree inter- 
relating charge configurations can be established by invent- 
ing or identifying a means of compensating for some unde- 
sirable feature of a chosen design, extending the principle 
(i.e., introducing redundancy) or inverting the princi- 
ple. 


Introduction 


FUNDAMENTAL characteristic of solid propellants is 
change in operating chamber pressure P as a function of 
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propellant burning surface (with nogzle area constant) 
according to the relation ; 


P, and P: are pressures corresponding to two burning sur- 
faces A; and Az of a propellant charge, and n is the pressure 
index or burning-rate exponent of the propellant. The pres- 
sure index is characteristic for a particular propellant and is 
usually less than unity, generally ranging in value from about 
0.4 to 0.7. It is apparent that pressure varies as about the 
square of the ratio of any two burning surfaces—say, initial 
and final—and may vary as the cube or higher power. A 
charge design with 10 per cent variation in burning surface 
would display a pressure variation of from 17 per cent (pro- 
pellant with a pressure index of 0.4) to 37 per cent (pressure 
index of 0.7). Thus, because of the magnifying effect of 
the exponent 1/(1 — n), allowable variations in propellant 
burning surface are generally limited to about 10 per cent or 
less since, except for certain special applications, it is usually 
desirable for a solid rocket to burn at approximately constant 
chamber pressure. Designs whose burning surface are ex- 
cessively progressive or regressive are usually not suitable. 


b INTERNAL-BURNING STAR 
Fig. 1 Cross sections of grains illustrating propellant burning in 
parallel layers 
A second fundamental characteristic of solid propellants is 
burning in parallel layers (Fig. 1). This property, together 
with change in chamber pressure as a function of burning 
surface, establishes one of the principal areas of analysis in 
the design of any solid propellant rocket: selection of a suita- 
ble geometry for the propellant charge. 


Scope of the Discussion 


The discussion of charge designs which follows will be 
limited to strictly geometrical considerations. No account 
will be taken of specific applications, if any, for any charge 
design. The effect of physical properties of the charge design, 
such as stress concentrations, resistance to deformation under 
load, or ability to withstand temperature changes, will be 
omitted. No consideration will be given to performance 
factors, such as methods of attaching charges to combustion 
chambers or efficiency of utilization of available chamber 
volume. Certain practical problems, such as methods of 
manufacturing charges or effects of gas flow (especially at 
very high velocity) on propellants or on combustion chambers 
will be purposely overlooked. It is, of course, recognized 
that in actual practice these other considerations are at least 
as important as the selection of a charge geometry. Finally, 
no structurally composite or nonhomogeneous propellants 
will be considered. 
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t) Morphological Approach 


Analysis of the geometry of propellant charges ranges from 
the extremely simple, as with the end-burning grain or the 
internal-burning tubular grain, to the extremely complex, as 
in the case with the internal-burning star or the slotted 

7 cylinder (Fig. 2). Equations have been developed which 
describe the behavior of the burning surface of various 


Me charge designs, and many configurations have been studied 
it and subjected to extensive testing. Most charge designs are 
‘4 based upon an intrinsically neutral burning shape (as the 
l internal-external-burning tubular grain, Fig. 2), or are de- 
A rived from the fundamental internal-burning tubular grain 
: by attempting to compensate in various ways for its intrinsic 
7 progressity. Restrictions (inhibitors) are employed to pre- 
vent certain surfaces from burning. 

f Before discussing the specialized method of interrelating 
t charge designs which is the main purpose of this paper, it is 
e helpful to establish as a convenient datum a pure morphologi- 
eal classification. It is apparent that such a classification 
‘ can be developed on the basis of the several distinguishing 


characteristics of charge designs. A little reflection will 
reveal that one way of classifying these characteristics is the 
following: 

1. ‘Principal directions of burning: 1.1 longitudinally; 
1.2 two longitudinally; 1.3 many longitudinally; 1.4 radially 
out; 1.5 radially in; 1.6 two radially; 1.7 many radially; 


2a. End-burning grain 


2b. Internal-burning tubular grain 


2c. Internal-burning star 
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1.8 laterally; 1.9 two laterally; 1.10 many laterally. 


2. Conirol of burning surface: 2.1 uncontrolled (un- 
restricted); 2.2 controlled (restricted). 

3. Character of main perforation cross section: 3.1 uniform; 
3.2 nonuniform. 

4. Character of ends: 4.1 both flat; 4.2 one conical inside; 
4.3 two conical inside; 4.4 one conical outside; 4.5 two 
conical outside; 4.6 one hemispherical inside; 4.7 two hemi- 
spherical inside; 4.8 one hemispherical outside; 4.9 two 
hemispherical outside; 4.10 one ellipsoidal inside; 4.11 ~ 
two ellipsoidal inside; 4.12 one ellipsoidal outside; 4.13 two 
ellipsoidal outside. 

5. Number of perforations: 5.1 none; 5.2 one; 5.3 two; 
5.4 many. 

6. Number of grains: 6.1 one; 6.2 two; 6.3 many. 

In this classification as with certain primitive tribes, any 
number greater than two is arbitrarily considered as many. 
The matrix in Fig. 3 shows the morphological identification 
of the end-burning grain, the internal-burning star, and the 
cone-and-cylinder charge (Fig. 2). 

It appears that at least 6240 charge designs are available 
from the characteristics selected. In actual fact, even more 
designs are available, since the internal-burning star is 
capable of enormous flexibility by selection of the number of 
star points (two or more), the star diameter, and the angles 
which uniquely define any given star. 


RESTRICTOR 


2d. Slotted cylinder 


RESTRICTOR 


2e. Internal-external-burning tubular grain 


2f. Cone and cylinder 


Fig. 2 Representative charge designs 
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Fig. 3 Morphological identification of three charge designs 


While many of these possibilities must be excluded from 
serious consideration because of their obvious impracticality, 
a large number of feasible designs remain, many of which are 
untested. The morphological approach not only generalizes 
the whole field of charge design, but also is capable of ex- 
pansion in light of new concepts which may be added to this 
area of analysis. 


Quasi-Morphological Approach 

After briefly reviewing the pure morphological classification, 
it is both instructive and interesting to examine a quasi- 
morphological approach which, though specialized, seems to be 
peculiarly well adapted to the particular field of charge 
designs. Remembering that most designs are derived from a 
few basic ones, usually by attempts to compensate for in- 
trinsic progressivity or regressivity, a family tree can be es- 
tablished by means of the following steps: 

1 Invent or identify a means of compensating a chosen 
design. 

2 Extend the principle involved where possible; i.e., 
introduce redundancy. 
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4b, External-burning cruciform grain 


3 Invert the principle involved where possible. 
Examples of means of compensation include the following: 


1 Increase the initial burning surface (internal star). 
2 Arrange for charge to shorten during burning (cone and 


cylinder). 
3 Allow burning to occur on additional surfaces (rod and 
tube, Fig. 4). 


4 Apply inhibiting material to certain surfaces (external- 
burning cruciform, Fig. 4). 

5 Employ combinations of these methods. 

Examples of extension or inversion of a principle include 
burning a grain on two or more ends instead of one, or adding 
additional radial burning surfaces to a rod and tube to 
achieve a multiperforated grain (Fig. 4) or changing an in- 
side conical end to an outside one (Fig. 4). A chart partially 
illustrating these interrelations is shown in Fig. 5. Illus- 
trations of additional charge designs mentioned for the first 
time in Fig. 5 are shown in Fig. 6. 

It is apparent from the chart that many configurations can 
be derived from more than one parent design. For instance, 
the end-burning grain becomes a double end burner and then 
a multiple-disk charge by progressive introduction of re- 
dundancy. The slotted cylinder (longitudinal slots) be- 
comes a multiple disk by inversion. The external-burning 
tube becomes a slab by introducing external longitudinal 


PROPELLANT 


RESTRICTOR 


4d. Multiple disk grain 


PROPELLANT 


RESTRICTOR 


4e. Frustum and cylinder 


Fig. 4 Examples of means of compensating chosen charge designs and extending or inverting a principle 
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Fig.5 Chart showing interrelationship of several charge designs 


6a. Internal-burning screw 


6b. Double end-burning grain 


6c. External-burning star 


ridges; the slab is also derivable from the internal-burning 
star by inversion of internal longitudinal ridges. The slab, 
by redundancy, becomes a multiple slab which, by inversion, 
becomes again a multiple disk. 


Conclusion 


In practice, a particular rocket design problem may be 
rather strictly circumscribed by various dimensional, perform- 
ance, and material limitations. However, the selection of a 
suitable charge geometry for a particular rocket usually pre- 
sents the designer with an opportunity to exercise his in- 
genuity and judgment. This freedom of selection is an 
indication of the intrinsic flexibility of solid propellant charge 
designs. The quasi-morphological approach to geometry of 
solid rockets, while not so generalized nor so powerful as the 
morphological one, is quite useful in developing new charge 
designs by logical extension of existing practice, and may be 
helpful in stimulating re-evaluation or modification of some 
designs widely used at present. 
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Fig. 6 Additional charge designs showing extension or inversion of a principle 
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The Optimum Ratio of Propellants for 
a Liquid Bipropellant Rocket Operating 
Within a Mixture Ratio Tolerance 


J. A. BROUSSEAU, JR.! 
Boeing Airplane Co., Seattle, Wash. 


For any liquid bipropellant rocket system designed for 
propellant exhaustion, close contro] of the burning mix- 
ture ratio and weight ratio of loaded propellants must be 
maintained if sizable quantities of unburned propellant 
are to be avoided at shutdown. This unused propellant re- 
duces the total available impulse from the rocket, and a 
consequent reduction in the over-all performance of the 
flight article in which it is installed. This paper considers 
the problem of unburned propellant resulting from a bi- 
propellant rocket operating over a range of mixture ratios. 
Equations are derived from a mathematical analysis which 
provides the rocket designer with tools for optimizing the 
weight ratio of propellants for any given operating mixture 
ratio tolerance. 


Nomenclature 
W. = available oxidizer propellant in lb 
W; = available fuel propellant in lb 
W,. = total available bipropellant supply in lb = Wo + W, 
R = ratio of available oxidizer to fuel propellant by weight 
W./ W; 
Wo = rate of oxidizer burned in Ib/sec 
wy = rate of fuel burned in Ib/sec 
w: = total bipropellant flow rate in lb/sec = wo + wy 
r = operating mixture ratio of rocket motor = wo/wys 
t = burning time in sec 
L = quantity of unburned propellant at end of any firing in 


Ib 

E(L) = average or expected value of L (Ib) 
g(r) = statistical distribution function for r (nondimensional) 
Definitions 

Available propellant—The total propellant serviced to the sys- 
tem less the quantity of propellant unavailable to the rocket mo- 
tor under any operating condition. 

“Specific impulse’””—The impulse delivered by the rocket motor 
per unit weight of propellant consumed. 


Introduction 


'HE maximum impulse to gross weight for any liquid bi- 
propellant rocket system is obtained when the following 
conditions are satisfied: 

1 The rocket motor is operated at the mixture ratio, oxi- 
dizer to fuel by weight, which will yield the highest “specific 
impulse” from the system propellants. 

2 All of the propellants available for producing impulse 
are burned. 

Condition (1) is easily satisfied, since thermochemical 
calculations will show the operating mixture ratio which gives 
maximum specific impulse or optimum rocket motor per- 
formance for a given chamber pressure and expansion ratio. 
To satisfy condition (2), the operating mixture ratio of the 
rocket motor and the weight ratio of available propellants, 
oxidizer to fuel by weight, must be the same. Any difference 
which exists between the above ratios will result in one pro- 
pellant, oxidizer, or fuel, remaining in one of the tanks when 
the other is exhausted. In practice, the possibility of obtain- 
ing simultaneous exhaustion of propellants from one firing to 
another of the same or similar systems is difficult, if not im- 
possible, due to the variations that exist in those factors 
affecting the operating mixture ratio. Consequently, tol- 
erances are generally specified on the operating mixture ratio 


Presented at ARS 25th Anniversary Meeting, Chicago, IIl., 
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to control these factors in design. In any event, since the 

operating mixture is allowed to vary, it is essential that the 

weight ratio of available propellants be properly fixed such 

that the probable loss of unburned propellant is minimized. 
Unburned Propellant 


The total available propellant supply in the system can be 
expressed as 


Ww, 
Wet = Wy +1) = + 1] 


or W,= WA 


Similarly, the total propellant flow rate to the rocket motor 
can be expressed as 


w,(r + 1) 


r 


or = 


For a given propellant supply and flow rate, the maximum 
burning time available from the fuel is 


=, \R+1 w (R+1) 


and from the oxidizer 


_W,_ WR(rt+1 pW. (1+ 


To exhaust both propellant supplies simultaneously, the burn- 
ing time of the fuel and oxidizer must be the same. Therefore, 
from Equations [5] and [6] 


then R=r 
which verifies the statements noted in the introduction. 

Now, if in a particular run r < R, then from Equations [5] 
and [6], t; < t; therefore, fuel will run out first and oxidizer 
will remain in the system. From Equations [4] and [5] the 
amount of oxidizer consumed within this period is 

wir (pt We 

\R+1 R+1 


= 
From the above and Equations [1] and [2] the amount of un- 
used oxidizer for the run is 


Wir W.R Wer 
R+1 R+1 k+1 


or L = W;,(R — r) whenr<QR............. [7] 


L=Wo- 


In a similar manner, if r > R, then from Equations [5] and 
[6], t, < ty; therefore, oxidizer will run out first and fuel will 
remain in the system. From Equations [3] and [6] the 
amount of fuel consumed within this period is 


BG ( __ WR 
he = Loe \R+1)) 
From the above and Equations [1] and [2], the amount of 
unused fuel for the run is 


WR 
L=W, when r > [8] 
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Minimum Loss for a Finite Operating Mixture Ratio 
Interval 
The loss function represented by Equations [7] and [8] has 
the shape shown in Fig. 1. 


R r 


Fig. 1 


As r approaches R, the amount of unused oxidizer, repre- 
sented by the straight line to the left of R, and the amount of 
unused fuel, represented by the curve to the right of R, ap- 
proaches zero. If the operating mixture ratio r is restricted 
to a finite interval, say rz <r < rg, then to approach the con- 
ditions leading to minimum unused propellant, R must also be 
restricted to the interval r, < R < ry. The loss function for 
different values of R between the interval is shown in Fig. 2. 


UNUSED OXIDIZER 


From Fig. 2, as R approaches rz, the amount of unused fuel 
for any operating mixture ratio between r; and rz is approach- 
ing zero, whereas the amount of unused oxidizer is increasing. 
Conversely, as R approaches rz, the amount of unused oxidizer 
is approaching zero, whereas the amount of unused fuel is in- 
creasing. It follows from the above that the unused propel- 
lant for the range of operating mixture ratios rz to rz will be 
minimized if R is chosen so that the maximum loss in either 
direction from R is the same. In Equation [7], let r = rz. 
In Equation [8], let r = rx. Equating the above, we get 


W,(R — rr) = W;[1 — (R/ru)] 


or (1 + rz) 

R=ryz 
Proportioning the available propellants to the value obtained 
from Equation [9] will minimize the unburned propellant 
losses possible within the normal range of operating mixture 
ratios rz, to ry. The usefulness of Equation [9] can be best 
illustrated from the following assumed example. 

The rocket motor for a certain liquid bipropellant rocket sys- 
tem is allowed to operate between mixture ratios, oxidizer to 
fuel by weight, of 2.90 to 3.10. Determine the unusable pro- 
pellant losses for off-mixture operation within the tolerances 
and the optimum design ratio of available oxidizer to fuel by 
weight for the system. 


Solution 

If we assume that the ratio R of available oxidizer to fuel by 
weight lies between the operating mixture ratio tolerance, 
then the unburned propellant loss for off-mixture ratio opera- 
tion can be best shown by Fig. 3. 

From Fig. 3, the maximum unburned propellant loss for 
any value of R within the tolerances can be represented by 
Fig. 4. This curve shows that the optimum design ratio of 
available oxidizer to fuel by weight lies at 2.948, which is also 
the value obtained from Equation [9]. 
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Fig. 4 
Statistical Considerations 

Referring to the unburned propellant losses represented by 
Equations [7] and [8], these,equations suggest that since r has 
a range of variation, it should be treated as a random variable 
whose behavior is described by some distribution function g(r) 
which can be determined by a statistical analysis of experi- 
mental observations. 

In order to determine the distribution of values of r, and so 
the function g(r), it is necessary to observe many individual 
values of r. These values might be obtained from many runs 
of a typical engine or from runs on many engines. The relia- 
bility of the results on any hypothesis, of course, increases with 
the number of observations. The greater the dispersion of the 
observed values, the greater the number of observations 
needed to yield the reliability desired. 
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The distribution function g(r) is defined by requiring the 
probability that r be in the short interval r; < r <r; + dr; be 
g(r:)dr;. The probability that r lies in an arbitrarily long 
interval r; < r < then 


f g(r)dr 


The graph of this distribution function can be drawn as fol- 
lows. The whole range of values of r is divided into a con- 
venient number of subintervals. On each subinterval g(r) is 
given a constant value equal to the fraction of the whole num- 
ber of observations found in this subinterval. For example, if 
15 out of 150 observations fall between r = 2.90 and r = 2.92, 
then g(r) is given the value 0.1 in this interval. 

Extreme caution should be exercised in identifying a sample 
distribution with the true distribution of r. Statistical 
methods exist for the estimation of the true distribution and its 
parameters in terms of the sample distribution to any degree 
of reliability desired. 

Considering g(r) to be defined for all positive values of r, 
even though it may be zero over a part of this range, we have 


and g(r) >0 
~ Thus the probability of finding r somewhere among the 
positive numbers is one, and probabilities are nonnegative 
real numbers between zero and one. 
Minimization of the Loss 

The average or mathematical expectation of L, designated 
as E(L), is given by summing all of the various values of L for 
given values of r times the probabilities that those values of r 
will occur. In symbols 


B(L) 
Now inserting L in terms of W;, R and r yield 
or 


R must now be chosen so as to minimize E(L). Setting to zero 
the derivative of E(L) with respect to R 


R 
f f g(r) [11] 
0 


then determines R, the weight ratio of available propellants 
which should be used to minimize the expected loss of un- 
burned propellant. It is assumed that the distribution func- 
tion g(r) for the mixture ratios observed for individual engine 
runs is known or has been determined by statistical analysis of 
observations as previously mentioned. Examples of results 
found from Equation [11] for several plausible forms of 
g(r) follow. 


Illustrations 
Example 1 Uniform Distribution 
In this case all values of r on a certain range, say m — d < 


r < m+ d are considered to be equally probable. Then, 


from Fig. 5 
0 r<m-—d 
g(r) = m—-d<r<m+d 
0 m+d<r 


Substitution into Equation [11] gives 


R ors. 
f car -dr = 0 
m-—d R r 
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Integrating these two expressions yields 
=m —d+In(m + d) 


The numerical solution of this equation gives the design tank 
mixture ratio R for the minimum loss for the case of the uni- 
form distribution of r. 

If the range of r is +0.10 about a mean value of 3.00, then 
m = 3.00, d = 0.10, and R is found to be 2.950. 

If d is much smaller than m, the equation above can be 
simplified to read 


m— 


In the above numerical example this yields R = 2.950. It 
should be noted that the above results and the results obtained 
from the minimum loss analysis are equivalent. This condi- 
tion will hold true if the range in mixture ratios is small com- 
pared to the ratios themselves. 


Example 2 Parabolic Distribution 
Again r lies in a certain range of variation but has a higher 
probability of lying in the middle portion of the range than it 


+¢ 


m-d m r 
Fig. 6 
does in regions about the end point of the range (Fig. 6). 
0 r<m-—d 
g(r) = [d? — (r — m—d<r<m+d 
0 m+d<r 


Substituting in Equation [11] gives 
d 
[d? — (r — — = 0 


m—d 
where c has been divided out. 
Carrying out the integration we have 


a(R — m + d) — (R — m)§/3 — + (m? — d?) In 
(m + d)/R + 2m(R — m — d) + (m + d)?/2 — R*/2 = 0 


As in the previous example, we may solve this expression for 
R when numerical values are assigned to m and d. 
Again if d is much smaller than m, the equation may be 
simplified to 
(m + 1)(R — m)* — 3(m + 1)d2(R — m) — 2(m — 1)d? = 0 
Now (R — m) is necessarily of smaller magnitude than d. 
If (R — m)? is much smaller than 3d?, the first term in the 
above may be neglected and the equation solved explicitly for 


R Rem- 
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With the same numerical example as before, m = 3.00 and 
d = 0.10so that R = 2.967. (Since d is here much smaller than 
m and since (R — m)* = 0.0011 is much smaller than 3d? = 
0.0300, the use of the last approximate equation above is 
justified.) 
Other Distributions 

Several other distribution curves such as the “normal” dis- 
tribution, “gamma”’ distribution, etc., have been devised in 
statistical analysis which may more realistically describe the 
probability density function of operating mixture ratio. How- 
ever, the formulas for these more realistic distributions are not 
amenable to mathematical manipulation as readily as the dis- 
tributions noted in the above examples. In any event, for 
any assumed symmetrical distribution of r about the mean m, 
the solution for R would show that 


R< mwhenm > 1 


R>mwhenm< 1 


The above statement is based on the results of the illustrated 
distributions and the fact that R approaches m when the 
density of the probability function at the mean approaches in- 
finity. 

Conclusions 


Results of this paper indicate that, for any liquid bipropel- 
lant rocket operating within a mixture tolerance, the weight 
ratio of available propellants and the mean operating mixture 
ratio should not be the same. The only exception to the above 
is when the mean operating mixture ratio is 1.0. In practice, 
where the mixture ratio range is small compared to the ratios 
themselves, and the ratios are above 1.0, the equations from 
the minimum loss analysis and the statistical distributions 
noted in this paper yield about the same result, namely, that 
weight ratio of available propellants should be somewhat less 
than the mean operating mixture ratio. For operating mix- 
ture ratios less than 1.0, the weight ratio of available propel- 
lants should be somewhat greater than the mean operating 
mixture ratio. 

In the design stages of a new system, where no test data are 
available for determining a suitable distribution of operating 
mixture ratios, it is recommended that the ratio of available 
propellants obtained from Equation [9] be employed. A de- 
sign based on the above may not be optimum, but will insure 
the minimum possible loss of propellant under any condition of 
operation within the tolerance. 
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Longitudinal Forced Vibrations of 
Cylindrical Fuel Tanks 
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HIS note is concerned with the pressure field in the fluid 

and the response of the wall of a cylindrical elastic tank 
due to small longitudinal forced vibrations. It will be seen 
that simple approximate expressions can be obtained, even 
for geometrically complicated configurations of the tank 
bottoms, e.g., spherical or conical bottoms. 
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Fig. 1 shows the cylindrical part of the tank of length L, 
radius R, thickness h, and the longitudinal and radial co- 
ordinates z and r, respectively. The origin of z is selected at 
the surface of the fluid. Assuming that the circular fre- 
quency 2 of the excitation is small vs. c/L—where c is the 
velocity of sound in the fluid—it is sufficient to use the theory 
of incompressible fluids. The velocity u and pressure p can 
be derived from a potential function ¢ 


u = — grad ¢, 


where ¢ satisfies 


As a first step, particular solutions of this equation satis- 
fying the boundary conditions everywhere except at the 
bottom of the tank are obtained. Functions 


F = sn( 4 [3] 


satisfy both, Equation [2] and the boundary condition p = 0 
on the free surface, z = 0, for any real or complex value of yu, 
where J) are modified Bessel functions. The boundary con- 
dition on the cylindrical surface of the tank r = R is ob- 
tained from the fact that the radial velocity of the fluid 


must be equal to the velocity of the tank wall. This bound- 
ary condition could be derived by considering the tank 
as an elastic shell; simpler but approximate expressions will 
be obtained here by neglecting the coupling of longitudinal 
and radial motions. Assuming the cylindrical wall to consist 
of isolated hoops which can displace radially only, the equation 
of motion for the displacement w of the wall of the tank be- 
comes that of a harmonic oscillator 


[5] 


where p; is the mass density of the tank material and w the 
extensional frequency of a ring of a radius R 


Noting that w will be of the form w = C sin [u(z/R) Je 
one finds by substitution that the Expression [3] will satisfy 
the boundary condition for r = R for any value of » which 


satisfies the relation 
To(u) _ hor (2 1) [7] 
uli Rp \Q? 


This equation has an infinite number of roots. If 2 < w, one 
root is real, the rest are imaginary; if 2 > w, all roots are 
imaginary. For the imaginary roots it is convenient to 
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change Equation [3] by using hyperbolic and regular Bessel 
functions. 

In the case of metal tanks of practical dimensions the tank 
frequency w is sufficiently high such that only forcing fre- 
quencies 2 < w are of structural importance and will be con- 


sidered here. (To study cases 2 > w it would also be necessary 
to use compressible theory as the assumption 2 < c/Z would 
not be satisfied.) In the cases to be considered, Equation [2] 
has the following series solution satisfying all the boundary 
conditions, except those for z = L 


= 
E sin + >) Cn sinh 
n=1 


where the values u, are the absolute values of the roots of 
Equation [7]. 

The coefficients C’,, can be utilized to complete the problem 
by satisfying conditions at z = L caused by the motion of the 
bottom of the tank. For this purpose it is important to 
realize that the sequence of functions Jo, Jo defined by the 
roots of Equation [7] are the solutions of a well-known eigen- 
value problem,? and form a complete set of orthogonal func- 
tions. 

The principal point of this note is the observation that the 
first term of Equation [8] contains the major part of the 
response at points not too close to the bottom of the tank. 
This is due to the fact that the terms containing sinh (u,z/R) 
are of importance only if z is nearly equal to L, but become 
quite small* ifz<< ZL. Simple, but good approximations for 
fluid pressures and tank response, valid for z < L, can be 
obtained by using only the first term of Equation [8]; if 
2 < w, the root 4; is necessarily larger than 3.83 (which is the 
first root of J:(u) = 0) and the neglected terms are significant 
only for z >(L — R/2);i.e., very close to the bottom. 


Application to Tank With Rigid Flat Bottom 


Consider first the rather unusual case of a flat rigid bottom, 
Fig. 2. Let the vertical velocity ug of the tank bottom be 


where U is the amplitude of the imposed forced vibration. 


s=L 

for all values of the radius r. Utilizing the orthogonality 
property of the Bessel functions in the Series [8], the coeffi- 
cients C, can be determined in a routine manner; using 


R R?2 
Ho 


R R? 
r1o%(por/R)dr = — — 


* Courant-Hilbert, ““Methoden der mathematischen Physik,” 
first edition, vol. 1, p. 254, J. Springer-Verlag, Berlin, 1924. 

* In the corresponding problem of wave ropagation in a tube 
containing an incompressible fluid, only the root yo furnishes a 
propagating wave. 
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the required leading coefficient Co in Equation [8] becomes 


21; 1( Mo) 
12 
uo) — cos /R) [12] 
Similar expressions apply for the other coefficients C,, and can 
be obtained by replacing uo by tun. Except near the bottom 
of the tank, it is sufficient to use just the first term of the 
Series [8] 


@ Cosin (yoz/R) Io (wor /R)e*™........... [13] 
Application to Tank With Rigid Bottom of Arbitrary Shape 


In case of an arbitrarily shaped rigid bottom, it is possible 
to obtain at least an approximate value for Cy by considering 
the total amount of fluid passing through the plane z = L. 
As the bottom is rigid and the fluid incompressible, the veloc- 
ity uz given by Equation [9], while no longer the actual 
one, is still the average one. Using the Expression [13] for the 
potential, one can determine an approximate value of Cy 
from the condition that the average vertical velocity of the 
fluid must equal uz 

RU 


This result can also be applied to a flat bottom, but is, of 
course, cruder than Equation [12]. 


Application to Tanks With Elastic Bottoms 


The case of a tank being closed by a flat elastic plate could 
be treated easily, and expressions for all coefficients C,, can be 
obtained, because the Bessel functions in Equation [8] lend 
themselves to an expansion of the response of a circular plate. 
The case is however academic; practically important cases 
of other shapes can only be treated approximately based on 
Equation [13], as shown hereafter. 

To obtain a condition for the coefficient Cy in Equation 
[13], it is again specified that the average vertical velocity 
in the plane z = L in the cylindrical tank and in the bottom 
should match. The fluid being incompressible, the difference 
between the average vertical velocity in the plane z = L and 
the imposed velocity ug is equal to the change in volume of 
the tank bottom divided by 7R?. As a first approximation, 
valid only if the forcing frequency Q is sufficiently below the 
frequency of the tank bottom by itself, one can use the static 
relation between the pressure and the volume change of the 
bottom. This relation may be written 


where AV is the volume change, p the pressure change, and 
the value K depends on the geometry of the bottom; for a 
hemispherical bottom, Fig. 3, one finds easily 


where v is Poisson’s ratio. Other cases can also be handled 


Fig. 3 


without much difficulty. The average velocity @ and the 
rate of change of the pressure p, at z = L due to the potential 
¢, Equation [13], are 


008 (wo: Ia (wo 


sin (uoL/R) (wo)e*™ 


Ho 
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The continuity condition stated in the foregoing requires 


av 


and the coefficient Cy becomes finally 


RU 
2a cos(uol/R) 


Co = 
where 


Equations [14] and [19] differ only by the coefficient a 
which is unity for rigid tank bottoms. It is noted that a may 
become zero, indicating resonance. 

The use of the static value of K is only permissible if the 
forcing frequency @ is sufficiently smaller than the frequency 
wp of the tank bottom; the latter is of the order 


where M is the mass of the tank bottom including the fluid 
below z = L. Even if Equation [21] indicates that the use 
of the static value of K is not permissible, Equations [19] and 
[20] remain valid provided an appropriately corrected dynamic 
value of K is used. 


Effect of the Compressibility of the Fluid 


For sufficiently high frequencies, or large tanks, it may be 
necessary to allow for compressibility. Replacing Equation 
[2] by the appropriate equation A*g = (1/c?) g one can easily 
derive a refined general expression for the potential 


= 
sin Io (uor/R) + Cn sinh | eftt 


c c 


The retention of the term containing Q in the expression for 
jo introduces no complications at all. A more fundamental 
change occurs if the frequency 2 becomes so iarge that 
becomes imaginary; the hyperbolic function in Equation [22] 
changes in this case to a circular one, and one cannot draw 
the conclusion that the first term of Equation [22] is sufficient 
to express the response in the upper part of the tank. If 
w > Q the value uw, > 3.83, and this unpleasant situation 
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will not occur as long as 2 < 3.83 c/R; in practical cases only 
frequencies considerably higher than 200 cycles/sec will be 
affected. The correction to wo on the other hand may 
be necessary for large tanks where 2 ~ c/L; all results of this 
note apply in such a case provided the coefficient yo in all 
circular functions in Equations [12] to [20] is replaced by jo. 

Numerical Example. Fig. 4 shows the velocities of the 
tank wall for a rather large aluminum tank for a frequency 
2 = 500 rad/sec (about 80 cycles/sec) assuming a rigid 
bottom. Pertinent dimensions are also shown in the figure. 
The tank frequency w ~ 4000 rad/sec according to Equa- 
tion [6] is much larger than the forcing frequency 2. Equation 
[7] gives uo = 3.5; using c ~ 5000 ft/sec, the frequency Q is 
not small vs. c/Z, and one might suspect that the correction 
for compressibility, Equation [23], is appropriate. How- 
ever, computation gives jp = 3.52 ~ wo, showing that the 
correction is negligible. The amplitudes of the radial velocity 
w of the tank wall are 


2 coxuoL/R) Ue 1.76 sin (3.52/R) Ue™*. . . [24] 
It is seen that the velocities of the tank wall are of the same 
magnitude as the amplitude U of the imposed velocity. It is 
worth noting that Fig. 4 shows several nodes; this indicates 
that the selected frequency is considerably above the funda- 
mental frequency of the tank, which may surprise some 
readers, and resonant frequencies entailing large magnifica- 
tion exist close to the frequency considered. The assumption 
of a rigid tank bottom made in the foregoing is not realistic 
and in any actual case the correction according to Equation 
[19] would be indicated. 
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Jet Propulsion News 


Alfred J. Zaehringer, American Rocket Company, Associate Editor 


The Research Scene 


NEW type of thrust meter is being installed in the 
Boeing B-52 to give in-flight thrust readings. Made 
by Minneapolis-Honeywell and AiResearch, the engine pres- 
sure ratio indicating system takes the pressure at the front 
of the engine and from a pressure rake in the exhaust. These 
pressures are transmitted through tubing to an electrical in- 
strument where the ratio of one to the other is computed and 
transmitted electrically to an indicator in the cockpit. Each 
engine has its own system. 


@ Titanium research is in for expansion as the Department of 
Defense is to undertake a multimillion-dollar program devoted 
to titanium sheet rolling. Several years ago titanium sponge 
was $5/lb; it now costs betweeen $3.50 and $3.95/lb (fab- 
ricated parts may, however, be as high as $15-25/lb) and 
an eventual sponge price of $1-1.50 is being aimed at within 
the next few years. Major applications have been in military 
aircraft but potentialities for guided missiles are becoming 
greater. Titanium alloys promise to give properties of high 
strength steel at only 60 per cent of the weight of steel. In 
contrast to critical supplies of steel alloying materials, esti- 
mates place more than 100 million tons of titanium ore 
available in North America. Production is estimated at 
8000 tons in 1955 with a projected capacity of 22,500 tons in 
1957. Aircraft manufacturers alone have estimated that 
they could use 100,000 tons per year. Present high costs 
have been ascribed to the batch processes available, high 
scrap loss, and the difficult fabricating techniques. Some of 
the present titanium producers are: Cramet, Inc.; Dow 
Chemical; DuPont; Titanium Metals Corp. of America; 
and the Union Carbide and Carbon Corp. 


@ The ruggedness of the silver-zinc storage battery was dem- 
onstrated in the Viking 12 shot when such a battery sur- 
vived the 143-mile fall, was recovered intact, and is still in 
operating condition today. The battery was produced by 
the Yardney Electric Co. 


@ Two pilots have bailed out of jet aircraft at high speeds 
and have lived to tell about it. George Smith, U. S. 
test pilot, bailed out of his F-100A at 800 mph and at an 
altitude of 6500 ft. Subjected to 40 g’s, the pilot suffered 
from bloodshot eyes, hemorrhages, and intestinal wounds. 
His helmet, oxygen mask, shoes, and socks were ripped by 
the wind blast. His clothes were torn to ribbons. Another 
high speed bailout was made by an RAF pilot, Henry Hol- 
land, who jumped at 10,000 ft and a speed of 700 mph. Be- 
sides the traditional black eyes, the British pilot suffered 
fractures of the pelvis and arms. 


@ The use of welding in rocketry and nuclear reactions was 
the topic of a two-day meeting at the Armour Research 
Foundation of the Illinois Institute of Technology at Chicago 
on February 1-2. 


Air Conditioning Space Ships 


AS A result of present developments toward higher and 
higher human flights, Frederick H. Green of AiResearch 
(photo) has considered the fundamental characteristics of an 
air conditioning system for a space ship. His conclusions 
were: (1) oxygen may be supplied from thermos bottles of 
liquid oxygen; (2) carbon dioxide will be removed chemi- 
cally or by condensation; (3) other air impurities such as 
tobacco tars, cooking smokes (he doesn’t believe that space- 
men’s diet will be all pill), etc., could be condensed on cool 
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Garrett Corp. 


Giant high altitude pressure chamber can simulate altitudes 
beyond 65,000 ft 


surfaces and removed in a manner similar to a windshield 
wiper method; (4) to conserve energy, maximum use will be 
made of radiant energy transfer to and from external surfaces 
for heating and cooling; (5) automatic and coordinated 
temperature control for the entire ship will be needed; (6) 
forced circulation of air inside the ship will be needed; (7) 
automatic air pressure control of the entire ship is needed 
at all times; (8) increase or decrease of humidity will be 
controlled by a device which probably will operate in con- 
junction with a water economizer; (9) detection and repair 
of air leaks is essential; and (10) for certain flights, a vacuum 
cleaning apparatus may be required. 


Rockets and Guided Missiles 


@ Over 170 rockets will be fired by the U. 8. during the forth- 
coming IGY. These include AEROBEE, DEACON, and 
the ROCKOON. Items to be measured are: magnetic 
fields, ionosphere electrical currents, auroral particles, upper 
air composition, cosmic rays, solar x-rays, and ultraviolet 
radiation. 


@ The Netherlands’ new pilotless aircraft, AT-21, is called 
the AVIOLANDA. Powered by a French SNECMA pulse- 
jet of 190 lb thrust, use is slated as a radar tiainer and target 


‘craft. Launching is by rocket. 


@ A new auxiliary unit powered by liquid or solid propellants 
has been developed by AiResearch. Providing electrical 
and mechanical, or hydraulic power, to operate guidance and 
control system in guided missiles, the unit delivers 650 watts 
and can be up-rated to 1600 watts. The compact unit 
weighs only 29 lb and dimensions 13 X 10 X 7 inches (photo). 
Combustion of the rocket propellant furnishes gases which 
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New auxiliary power unit 


operate a turbo-generator rotating at 24,000 rpm. An out- 
put shaft and gearing can operate missile controls mechani- 
cally or turn a pump for hydraulic power up to 2.5 hp. 
AiResearch built one of the motors back in 1950 as a solid 
propellant gas generator to replace storage batteries in a 
missile. 

e@ Another gas generator uses isopropyl nitrate as a mono- 
propellant to start aircraft engines. The cartridge-type 
starter was developed by the British Thompson-Houston 
Co., Ltd., of Rugby and was recently demonstrated at the 
Farnborough Air Show. 


e A new rocket powered aircraft will probe the problems of 
manned flight at altitudes of up to 100 miles. North Ameri- 
ean Aviat‘on will work on the joint NACA-Navy-USAF 
project. 


@ Rockets were recently put to a novel—and effective—use. 
A five-ton military truck set a speed record of 107 mph with 
the aid of 6 RATO’s. The tests were made near Ipswich, 
England, to test out a new USAF nylon crash barrier devel- 
oped for jet fighters. 


Facilities 


OCKET engine development and production will be 
consolidated into Reaction Motors, Inc., new $4 million 
facility at Denville, N. J. (photo). The 200,000 sq ft area 
centers RMI’S administration offices, research, engineering, 
and manufacturing operations which were previously divided 
between Rockaway and Lake Denmark. 


e@ North American Aviation’s efforts in the electromechanical 
field will be vested in its Autonetics Division. From its 


RMI 
Aerial view of RMI facilities at Denville, N. J. 
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guided missile efforts over the past ten years, electromechanics 
now employs nearly 8000. Facilities are to be located at 
Downey. Another new North American group, the Rocket- 
dyne Division, was created for rocket work. A $3.475 million 
building for the group is nearing completion in the Canoga 
Park section of the San Fernando Valley. Headquarters 
will be near the NAA Field Laboratory in the Santa Susana 
mountains. Rocketdyne now employs 4000. 


@ Douglas Aircraft is separating its missile engineering func- 
tion from its aircraft engineering group. This was done in 
recognition of the increasing proportion of engineering effort 
being devoted to guided missiles. 


@ A 1575-mile long underwater phone cable, extending from 
Cape Canaveral, Fla., to Mayayuez, Puerto Rico, was com- 
pleted for the USAF’s guided missile test range. Called 
the world’s longest, the undersea cable will feed missile flight 
data to Florida. 


@ Soon to go into operation is a $41 million Navy laboratory 
at Trenton, N. J., for testing jet engines. Designed to com- 
plement the Air Force Arnold Engineering Center at Tulla- 
homa, Tenn., and the NACA Lewis Flight Propulsion Lab- 
oratory at Cleveland, Ohio, the installation features five 
huge concrete and steel test cells, heating and cooling equip- 
ment, and environmental equipment. To silence the turbo- 
jet and turboprop engines to be tested here, over $3 million 
was spent.” 


@ Palo Alto, Calif., is to get a new electronic research labora- 
tory to be devoted to the guided missile field. Admiral 
Corp’s Advanced Development Section will be housed here. 


@ Sylvania Electric Products has completed a new Avionics 
and Missile Systems Laboratory at Waltham, Mass. The 
120,000 sq ft facility will house 700 ciatl 
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U. S.Army 


Panoramic view of Picatinny Arsenal 


Rockets at Picatinny Arsenal 


LTHOUGH Picatinny Arsenal turned out cannon shot 
for Washington’s Continental Army, its rocket work 
started with the development of the 2.36-in. antitank rocket 
(BAZOOKA). Operated asa U.S. Army Ordnance facility, 
the Arsenal personnel totals more than 6500. Valued at over 
$76 million, Picatinny comprises a land area of over 5125 
acres with over 1200 buildings. Located in the rolling hill 
country of northern New Jersey, Picatinny is near Dover, 
some 40 miles west of New York City. 
Since the BAZOOKA, rocket activities have increased to 
a point where some recently developed rockets carry warheads 
of over 1200 lb (the BAZOOKA carried approximately 1 lb). 
After the BAZOOKA came the 4.5-in. Folding Fin Light 
Artillery Rocket (standardized as the M8). The rocket 
offered high rates of fire and absence of recoil and allowed one 
4.5-in. cluster launcher to produce a rate of fire equal to 3 


U. S. Army 
Double base powder propellant plant 


U. S. Army 


Firing of T50 Booster 


standard 105-mm howitzers. Another activity was concerned 
with the design and development of the sound and flash 
tracking system for the AEROBEE rocket and fuzing for 
ejection of instruments from high-altitude rockets. In the 
opening days of the Korean War, the Picatinny-developed 
3.5-in. BAZOOKA proved an effectivé antitank weapon. In 
the late 1940’s work was started on rocket boosters and cul- 
minated in the. development of an ATO unit and the T50 
Booster in 1953 for launching the B-61 Matador. In 1954, 
the T50 successfully launched a manned fighter plane from a 
zero length launcher. An unusual JATO for mine field 
clearance was also developed, called the SNAKE. 

Propellant research is also under progress at Picatinny. 
This is concerned with ignition, ballistic modifiers, stabilizers, 
storage, and the use of wood pulp to replace some of the 
cotton linters used in the manufacture of double base powder. 

Picatinny Arsenal, the cradle of American ammunition, is 
continuing to maintain the superiority of American rockets. 


PROPULSION. 
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Space Flight Notes 


Astronautics 


SPECTACULAR proposal for an inhabited satellite that 

ranks in boldness with the well-known earlier proposal 

of W. von Braun, was presented at the recent Annual Meeting 

of the ARS by Darell C. Romick of Goodyear Aircraft Corp. 

in a paper entitled ‘Preliminary Engineering Study of a 
Satellite Station Concept.” 

Mr. Romick described in great detail the construction of a 
giant manned space station. This is established in an orbit 
about the earth by large satellite ferry rockets which carry 
passengers and hardware to the steadily growing orbital sta- 
tion. Indeed, the first two ferry vessels themselves are not 
returned to earth, but form the embryo of the orbital station 
structure. It is estimated that about three years would be 
required to build a vehicle of 3-billion cu ft volume. The 
shape (photo) of this truly heroic vessel is cylindrical with a 
large rotating (1 rpm) wheel at one end. This wheel would 
provide living quarters for those passengers and crew who de- 
sire some feeling of ‘‘weight’’—which is obtained from the 
angular acceleration resulting from the wheel’s rotation. The 
nonrotating cylindrical body has living quarters for those who 
do not object to the weightlessness extant in an orbital vehicle. 


| 


ct 


This cylindrical segment is 1000 ft in diam and 3000 ft 
long; the wheel is 1500 ft in diam and 40 ft thick. The sta- 
tion has a scientific purpose for the study of solar energy, earth 
reconnaisance, and guidance of outgoing space ships; it also 
has a recreational purpose for tourists who wish to observe the 
splendors of the heavens or the earth from this vantage point. 


Comment 


Although this coming age of space flight may be expected to 
begin very modestly with a small instrumented satellite, it 
seems certain that larger unmanned and manned vehicles will 
follow in the future if the trend of growth in this field follows 
that of such other technical developments as aircraft and 
ships. Therefore a detailed speculation as in the discussed 
paper can be useful since it shows the magnitude of the engi- 
neering effort needed to build large inhabited space vehicles. 

Prognostication based on existing engineering knowledge 
can, of course, be misleading. The prophets of the middle and 
late 1880’s who could visualize the giant passenger sea vessels 
of the future, or great heavenly argosies of airships, did so in 
terms of their immediate engineering experiences. Thus the 
recorded descriptions of those times show that the great ocean 
vessels and airships of the future were “extrapolated”’ versions 
of existing vehicles such as the S.S. Great Western or the bal- 
loons in current use. 

Thus, when the age of great manned space vehicles arrives, 
the engineering experiences of our own days will surely seem 
crude. It is likely that methods of construction and propul- 
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Kurt R. Stehling, Naval Research Laboratory, Contributor 


sion in those distant times will be as advanced as the ideas of 
the giant spaceship itself. Indeed, if such vehicles are to be 
economically feasible, then novel and unorthodox methods may 
be necessary. 


REVIEW 


“Methods of Achieving Space Flight” (Wege zur Raum- 
schiffahrt), by Hermann Oberth. Published by R. Olden- 
bourg, Berlin, 1929. Now out of print. 

Part II—Chapters 9-11. (Continuation of review begun in 
November issue.) 

In Chapter 9, Oberth examined in great detail the phe- 
nomena of “Andruck” or “reactive force” or “reaction.” 
Andruck is a not easily translatable word which Oberth 
coined to describe the acceleration forces which can act on a 
body. Indeed, in one section of this long chapter he gives 
Andruck the same dimensions as acceleration. 

The physiological influence of acceleration, or deceleration, 
is treated; Oberth seemed greatly concerned with the effect of 
changes of velocity and position on the human organism. He 
illustrated with such homely examples as a carousel and a roll- 
ing or pitching ship the action of nonrectilinear motions upon 
the ear canal, the blood vessels, and internal organs. He was 
strangely prophetic in his assertion that motion sickness, ab- 
sense of weight, and sudden accelerations and decelerations 
would prejudice manned space flight. 

Large sums are being spent currently on various space medi- 

cine programs, including these very effects of rapid accelera- 
tion and decelerations and high “g” loadings on the human 
body. 
In Chapter 10, Oberth discussed ballistic trajectories and 
methods of overcoming the Earth’s gravitational attraction. 
Realizing that a minimum of 7 miles/sec would be needed to 
permit a body to escape the Earth, he postulated that three 
possible techniques would permit this: (1) The electromag- 
netic “‘cannon’’; (2) the electric ‘‘windmill” (described later 
in the book); (3) the rocket. 

He describes a horizontal cannon with a barrel 6000 miles 
long, consisting of segments of electromagnets, alternately N 
and S, which pull along a projectile to a muzzle velocity of 7 
miles/sec. It was evident that this enterprise as well as that 
of the electric windmill (to be described in a later review) 
were hopelessly impractical and that the application of rocket 
propulsion held out the only hope for escape from the Earth. 
The remainder of Chapter 10 then covers the computation of 
ballistic flight paths and the overcoming of the Earth’s gravity, 
using known and standard derivations of mechanics and 
ballistics. 

In Chapter 11 is outlined the effect of air drag on the rocket 
velocity. Oberth clearly realized that not only the shape but 
also the cross-sectional density (Ib/sq ft of missile area) were 
important in the determination of air drag losses. He showed 
that a vertical ascent by a ballistic vehicle is the most eco- 
nomical, when the velocity losses caused by air drag and grav- 
ity retardation are considered. 

In contrast, the useful application of air drag is shown for 
vehicles re-entering the atmosphere. Oberth mentioned air 
brakes and parachutes as likely devices for slowing down a 
rocket vehicle. He stated in a footnote that some of his air 
drag and parachute formulas have been “applied and proved” 
since the book was printed. A vague reference is given to a 
mail-carrying rocket idea which he developed. He shows also 
that an inclined or curved ascent rather than a vertical has 
merit only when a winged airborne rocket is used, since the lift 
which is naturally obtained can reduce the all-up weight of the 
vehicle even though it is designed to ultimately escape the 
Earth. 
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NEW 2-FOLD 


RETIREMENT PLAN 


for Aviation and Missile Engineers 


It’s natural that a long-time leader in aviation, like Republic— with many firsts to its credit — 
should be a leader, too, in providing for the welfare of its staff. Right in line with this forward- 
looking policy is the remarkable new Republic Retirement Income Plan. Here is how it works: 


PART 1 is a basic Retirement Income plan paid for entirely by Republic Aviation. 
PART 2 is a cooperative effort. It is completely optional. But if an engineer 
wishes to increase his retirement fund, by making a small monthly payment, 
Republic will MORE THAN MATCH his contribution. 


Take the case of a hypothetical engineer who joins Republic on January Ist, 1956, averages 
$8,000 a year for 15 years; then retires aged 65. If he elected PART 2 of the plan, he will have 
a total monthly retirement income of $225.80, including his social security. For this he himself 
will contribute only $8.50 a month to the Republic Retirement Income Plan. 


Of course, the MORE YOU EARN, the HIGHER your Retirement Income will be. And 
Republic pays a top salary scale in the industry. 


.» RETIREMENT PLAN JUST ONE OF MANY PLUS FACTORS POINTING TO A REPUBLIC CAREER 
FIRST —there’s the interest and prestige of the Research and Development Program has 
working for a pioneer in aircraft design, cre- just been announced. 
ater of such planes the P04 Thun- THIRD —an All-Expense Paid Relocation Plan 


derjet, the F-84F Thunderstreak, RF-84F 
Thunderfilash and XF-84H. (Soon to be fol- 
lowed on the production line by the new F-103, 
F-105 and RF-105.) 


SECOND —the company is expanding sharply, 
providing frequent opportunities for able men 
to advance. In fact a $12,000,000 increase in 


for qualified engineers living outside the New 
York City and Long Island areas, which makes 
it easy to move to Republic. Other liberal bene- 
fits: Life, Accident and Health Insurance; 
Hospital-Surgical Benefits for the whole fam- 
ily; educational aid covering %4 the cost of 
collegiate and graduate study. 


Serve your own best interests. Make full inquiries into the many advantages of joining 
Republic now, not the least of which is living on Long Island—the Playground of the East. 


Important engineering positions are now open at all levels: 


STAFF ENGINEERING SYSTEMS 
PRELIMINARY DESIGN MATHEMATICAL 
ELECTRONICS ANALYSES 
WEAPONS SYSTEMS AERODYNAMICS 
PROPOSALS TECHNICAL WRITING 
OPERATIONS RESEARCH DYNAMICS 


Please address complete resume 


AIRCRAFT 


Assistant Chief Engineer, Administration Mr. R. L. Bortner, 


Farmingdale, Long Island, N. Y. 


ELECTRO-MECHANICAL STRESS 
DESIGN AIRCRAFT DESIGN 
FLIGHT TEST WEIGHTS 
RESEARCH DESIGNERS: 
THERMODYNAMICS 
FLUTTER & VIBRATIONS Mechanical 


outlining details of your technical background to: 


AVIATION 


MISSILES 


Administrative Engineer, Mr. Robert R. Reissig 
Hicksville, Long Island, N. Y. 
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ARS News 


Los Angeles Conference on 


Satellites, Thermal Thicket, Air Cargo 


Three ARS Sessions Slated for March 14- 
16 Meeting with ASME Aviation Division 


A®S will cosponsor a three-day meeting 
with the ASME Aviation Div., March 
14-16, at the Hotel Statler in Los Angeles. 
ARS sessions, announced by Meeting 
Chairman David E. Shonerd of Aero- 
physics Development Corporation, will 
dea! with earth satellites, high tempera- 
ture materials, and design problems. 
ASME?’s Aviation Division will sponsor 
five sessions on Air Cargo and nine on the 
Thermal Thicket. Another meeting co- 
sponsor will be the ASME Instruments 
and Regulators Division which will stage 
three sessions on instrumentation. 
Discussed in Thermal Thicket sessions 
will be human factors, air frames, materials 
accessories, structures, and turbojets. 
Two luncheons will be held on March 
14 and 15 and a banquet on the evening 
of March 15. 
ARS sessions are listed below. The 
complete program including ASME ses- 
sions will be mailed to all members. 


WEDNESDAY, MARCH 14 

9:30 a.m. Earth Satellite 

Chairman: W. Bollay, Aerophysics De- 
velopment Corporation 

Vice-Chairman: T. Carney, Hughes Aircraft 
Scientific Aspects of the Satellite Program, 
W. H. Pickering, Jet Propulsion Lab. 

Project VANGUARD—The IGY Earth 
Satellite, Rear Admiral F. R. Furth, Naval 
Research Laboratory 


Feasibility of Orbiting Shaped Charge Parti- 
cles and Their Scientific Value, A. W. 
Ehler, Aerophysics Development Corp. 


THURSDAY, MARCH 15 


2:30 p.m. High-Temperature Application of 
Nonmetallic Materials 

Chairman: .F. Marion, Grand Central 
Rocket Company 

Vice-Chairman: W. Schaafma, Grand Cen- 
tral Rocket Company 

Use of Graphite in High Temperature 
Rocket Applications, L. D. Loch, Battelle 
Memorial Institute 

High-Temperature Short Lived Plastic Lami- 
nates, G. Dodson, H. I. Thompson Co. 

Pure Oxide Coating, W. M. Wheildon, R. 
J. Westerholm, Norton Co., and A. V. 
Levy, Marquardt Aircraft Co. 


FRIDAY, MARCH 16 


9:30 a.m. High Temperature Design Prob- 
lems 

Chairman: H. Schuerch, Aerophysics De- 
velopment Corporation 

Vice-Chairman: E. Knuth, Aerophysics 
Development Corporation 

Is There a Thermal Barrier for Supersonic 
Ramjet Tail Pipes? W. Unterberg, 
Marquardt Aircraft Company 

Structure Research at Heating Rates With 
Free Flight Rocket Vehicles, I. Zuckerman, 
Aerophysics Development Corporation 

Heat-Transfer Problems in Solid Propellant 
Rocket Design, T. T. Omari, Aerojet- 
General Corporation. 


Swanson Elected in Alabama, Succeeds Wiggins 


Conrad Swanson, third from left, new president of the Alabama Section of the American 
Rocket Society is congratulated by the outgoing president, J. W. Wiggins, Thiokol Chem- 
ical Corp., at the annual dinner meeting of the organization attended by some 125 per- 
sons Dec. 1 in Huntsville. Looking on at left is Cliff E. Fitton, who was installed as vice- 
president. All three men are employed at Redstone Arsenal, the Army Ordnance 
Guided Missile and Rocket Development Center. At right is Thomas L. Phillips, a 
guided missile expert from the Raytheon Manufacturing Company, Bedford, Mass., who 
was principal speaker for the meeting. Not in photo are the Chapter’s new secretary, 
William A. Davis, and new treasurer, Lawrence N. Nicastro, both of Redstone Arsenal. 


Frsruary 1956 


Report on Nuclear Science and 
Engineering Congress 
By ALEXANDER WEIR, JR. 
University of Michigan 


Ts Nuclear Science and Engineering 
Congress, held in Cleveland, Ohio, 
December 12-16, 1955, was sponsored by 
some 26 technical societies. Approximately 
5000 registrants attended this meeting, and 
over 350 papers were presented in the 50 
technical sessions. The session on Chemical 
Reactions and Catalysis was sponsored by 
the American Rocket Society and the 
American Institute of Chemical Engineers. 
The papers presented in this session are sum- 
marized below. 

The effect of gamma radiation on aqueous 
ethylene-oxygen solutions was discussed by 
Ernest J. Henley, Jane P. Schwartz, and 
Wesley P. Schiffries of Columbia University. 
Their results can be represented by the equa- 
tion 


CH.CH:; + 0: + H:O radiation 
CH;CHO + H,0 


In the ensuing discussion, Prof. Henley indi- 
cated that no evidence for the existence of 
ethylene oxide in the reaction products was 
found. D. E. Harmer, Joseph J. Martin, and 
Leigh C. Anderson of the University of Michi- 
gan presented chlorination rate data for aro- 
matic hydrocarbons under the influence of 
gamma radiation. Professor Martin indi- 
cated that, in the presence of gamma radia- 
tion, chlorine adds to the benzene ring to a 
much greater extent than it substitutes for 
hydrogen in the side chain. 

The results of experiments on the electro- 
lytic separation of metallic isotopes were re- 
ported by J. H. Andrews, Ivan Ceresna, 
F. A. Rohrman, and W. F. Utlaut of the 
University of Colorado. Their results indi- 
cated that, nickel-63 behaves more like a noble 
metal than other nickel isotopes. The use of 
nuclear radiation in combustion research was 
discussed by Martin E. Gluckstein, Richard 
B. Morrison, and John V. Nehemias of the 
University of Michigan. Mr. Gluckstein 
indicated that alpha and beta emitters would 
be expected to be more effective than gamma 
emitters, and he discussed the problems asso- 
ciated with the handling of such radioactive 
sources. The chairman of this ARS-AIChE 
session was Alexander Weir, Jr., of the Uni- 
versity of Michigan, while Hal F. Coats, 
Blaw-Knox Company, was co-chairman. 


Report on Second AGARD 
Combustion Colloquium 


By ALEXANDER WETR, JR. 
University of Michigan 


LARGE number of ARS members were 
among the 73 United States registrants 
for the Second AGARD Combustion Col- 
loquium, sponsored by the Advisory Group 
for Aeronautical Research and Development 
(AGARD) of the North Atlantic Treaty 
Organization (NATO). The colloquium was 
held December 5-9, 1955, at the Institute de 
Mécanique of the University of Liége, 
Belgium. Most of the other registrants were 
from the United Kingdom, France, Belgium, 
Western Germany, and the Netherlands, but 
Canada, Turkey, Norway, and Italy sent, 
delegations, and observers from Spain, 
Sweden, Egypt, Israel, and Australia were 
present. 
Presiding at the opening session, AGARD 
Chairman Theodore von Karman introduced 
A. Capetti, Ecole Polytechnique de Turin, 
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TO THE FINE ENGINEERING MIND 


SEEKING THE CHALLENGING PROJECTS IN 


ROCKET PROPULSION ENGINEERS are offered unusual career opportunities now 
at Convair in beautiful, San Diego, California, including: Design Engineers for 
design and analysis of advanced high performance rocket engine systems and 
components including propellant systems, lubrication systems, contro! systems, 
mounting structure, and auxiliary power plants; Development Engineers for 
liaison with Engineering Test Laboratories and Test Stations in the planning,. 
analysis, and coordination of rocket engine system and component tests; Devel- 


opment Engineers for coordination with Rocket Engine Manufacturers in the . 


installation design, performance analysis, and development tests in conjunction 
with Convair missile programs. Professional engineering experience in rocket 
missiles and aircraft propulsion system development will qualify you for an 
exceptional opportunity. 

CONVAIR offers you an imaginative, explorative, energetic engineering depart- 
ment... truly the “engineer’s” engineering department to challenge your mind, 
your skills, your abilities in solving the complex problems of vital, new, long- 
range programs. You will find salaries, facilities, engineering policies, educa- 
tional opportunities and personal advantages excellent. 


Generous travel allowances to engineers who are accepted. Write at once 
enclosing full resume to: 


H. T. Brooks, Engineering Personnel, Dept. 1414 


A Division of General Dynamics Corporation 
3302 PACIFIC HIGHWAY SAN DIEGO, CALIFORNIA 


SMOG-FREE SAN DIEGO, lovely, sunny city on the coast of Southern California, 
offers you and your family a wonderful, new way of life ...a way of life judged 
by most as the Nation's finest for climate, natural beauty and easy (indoor- 
outdoor) living. Housing is plentiful and reasonable. 


who summarized factors affecting aircraft 
engine performance at Mach numbers up to 
4. Discussing this paper, R. A. Tyler, 
National Aeronautical Establishment, (t- 
tawa, indicated that engines operating on a 
reheated turbojet cycle should have immedi- 
ate practical interest in the Mach 2 ranye. 
R. R. Hibbard and H. C. Barnett, NACA, 
presented the results of studies concerned 
with the effect of fuel volatility and composi- 
tion on both air frame and engine require- 
ments. 

B. P. Mullins, AGARD Combustion Panel 
Chairman, presided at the session on ignition 
where papers by B. Lewis and G. von Elbe, 
Combustion and Explosion Research, L. D. 
Wigg, National Gas Turbine Establishment 
(NGTE), M. Barrére and A. Moutet, Office 
National d’Etudes et de Recherches Aéro- 
nautiques (ONERA), and R. 8S. Brokaw, 
Princeton, were presented. Considera)le 
discussion of these papers by J. W. Linnett, 
Oxford; E. A. Watson, Joseph Lucas, Ltd.; 
A. Van Tiggelen, Louvain; W. F. Jost, 
Géttingen; M. Vernet-Lozet, Société d’Etude 
de la Propulsion par Reaction; and Jr. 
Mullins resulted, primarily concerning the 
excess enthalpy concept proposed by Ir. 
Lewis. 

The chairman of the session on transport 
properties was S. S. Penner, California Insti- 
tute of Technology. Dr. von Karman’s paper 
on aerothermochemistry was discussed by 
Professors Kirkwood (Yale) and Herzfeld 
(Catholic University). E. A. Brun, Faculté 
des Science, Paris, presented the results of a 
study of convection at high temperatures, and 
J. Hilsenrath, Bureau of Standards, sum- 
marized and correlated the transport coeffi- 
cients found in some 400 references. At the 
conclusion of his paper, Mr. Hilsenrath pre- 
sented a nine-volume library of thermody- 
namic tables to J. Ducarme, Deputy Chair- 
man of the AGARD Combustion Panel, for 
use at the University of Liége. 

High-altitude combustion problems were 
discussed in a session whose chairman was 
Louis Deffet, Institute Belge des Hautes 
Pressions. Papers by J. Soissons, Centre 
d’Essais en Vol; 8S. L. Bragg and J. B. 
Holliday, Rolls-Royce; S. Way, Westing- 
house; and H. G. Wolfhard, NGTE, were 
presented. The problems encountered in 
scaling turbojet, ramjet, and rocket combus- 
tion chambers were discussed by A. E. Weller, 
Battelle; D. G. Stewart, NGTE; B. Salmon 
and H. Vigne, Société National d’Etude et de 
Construction de Moteurs d’Aviation; and 
C. C. Ross, Aerojet. C Zwicker (Nether- 
lands) was the chairman of the latter session. 

In addition to the technical papers, movies 
were shown by A. O. Tischler, NACA; E. F. 
Zukowski and F. E. Marble, California Insti- 
tute of Technology; H. G. Wolfhard, NGTE; 
and J. Surugue, ONERA. The proceedings 
of this colloquium are to be published in a 
book edited by Prof. M. W. Thring, Uni- 
versity of Sheffield, Sheffield, England. 


Reagan of Beech Heads 
Wichita Section 


AMES F. REAGAN, manager of mis- 

sile engineering for Beech Aircraft 

Corporation, has been elected president of 
the Wichita Section. 

During a recent meeting of the newly 
formed group, other officers elected were: 
Lawrence J. McMurtrey, power plant staff 
engineer, Boeing Airplane Co., vice-presi- 
dent; Wilbur F. Pro, rocket propulsion 
systems engineer, Beech Aircraft, secretary; 
and George H. Ward, aircraft engines 
specialist, General Electric Co., treasurer. 

In addition to the new officers, those 
named to a nine-member board of directors 
include: 
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Henry F. Waring, chief military aircraft 
engineer, Cessna Aircraft Co.; Maj. 
Clarence H. Mills, wing safety officer, 
McConnell Air Force Base; Richard B. 
Holloway, aerodynamics staff engineer, 
Boeing; Alex N. Petroff, director of re- 
search, Cessna; and David Henderson, a 
student at Wichita University. 


First Meeting of 
“California Desert”’ 
Section 


NEARLY 200 attended the first meeting 
of the proposed ‘California Desert’’ 
Section at Mojave, November 30. 
Richard Schmidt, asst. chief of the 
Rocket Engine Test Lab., Edwards AFB, 
introduced G. P. Sutton, National Mem- 
bership Chairman of ARS, who talked on 
“Rockets Behind the Iron Curtain.”’ 


Section Forming 
At Columbus 


RS members and interested individuals 

from North American Aviation, Bat- 
telle Memorial Institute, and Ohio State 
University met on Nov. 17 to organize a 
section of the Society in Columbus. 
Rudolph Edse, project engineer, Rocket 
Research Laboratory (Ohio State Univer- 
sity) spoke on “Problems Related to 
Rocket Research.”’ A second meeting was 
held at Battelle on December 13. James 
J. Harford, Executive Secretary, outlined 
the objectives of ARS and spoke briefly of 
its history. William Beans, power plant 
test and development engineer (North 
American) talked on ‘‘A Method of 
Evaluating and Selecting Power Plants.”’ 
M. W. Jack Bell was elected temporary 
chairman and David B. Haswell, secre- 
tary-treasurer (both North American). 
The Executive Committee consists of 
Stephen Davis, Kirsey Smith, Donald 
Wertz, and William Welch (North Ameri- 
ean); Dick Bowman and Alex Lemon 
(Battelle), and Rudolph Edse (Ohio State 
University ). 


San Diego Section 
Under Way 


AN organizational meeting of a pro- 
posed San Diego Section was held on 
Dec. 20. 
Those interested in the Section’s activi- 
ties should contact the temporary chair- 
man, William E. Clark of Convair. 


Section Doings 


Arizona. A brief business meeting was 
held December 6, followed by a lecture on 
airfoil design by Richard Thompson. A 
film was shown and refreshments served to 
the 32 members and guests present. At a 
Jan. 10 dinner meeting, R. Sumner lec- 
tured and showed a film on rocket sleds. 

Indiana. J. Donald Haas (Reaction 
Motors) spoke at a Purdue University 
meeting, November 29, on the status of the 
missile program. Color movies were 
shown, one entitled “Exploring Space’’; 
the other was a film on NACA experi- 
mental rocket planes. 

National Capital. Refreshments and 
door prizes preceded a dinner meeting on 
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SANBORN Pe_mapapet 


is a Sanborn feature that 
is especially valuable in 
multi-channel recording. 
The user of a 2-, 4-, 6- 
or 8-channel “150” sys- 
tem can record interre- 
lated events simultane- 
ously on one strip of 
record paper, and then 
correlate them against 
time quickly and accu- 
rately, even during the 
recording. 


PROVIDE THESE 


SPECIFIC 
ADVANTAGES 


The three foremost 
operating advantages of 
a Sanborn “150” are + 
evident in the record, - 
and are given at the f--== 
right. Equally popular | 
with industrial users is 
the wide versatility of ©~" 
these systems. Each one 
comprises a basic assem- ~ 
bly to which the user | 
adds his choice of an 
interchangeable, plug-in 
type preamplifier for 
each channel. Presently 
available different pre- 
amplifiers now number - 
eleven. 


The nichrome ribbon 
tip “150” heated stylus 
removes the white 
opaque surface of San- ; 
born recording paper 

(Permapaper) permit- 
=) ting the black under- 
ayy coating to show. Trac- 
ings are permanent, will 
not fade or smudge. 


| 


Over-all Linearity 


; resulting from current 
Jee, feedback design of 
Driver Amplifiers and 
new shorted coil frame, 
high torque (200,000 
<j dyne cm.) galvanome- 
ters. Maximum error 
i over middle 4 cm. of 
chart: 0.25 mm. . 

over entire 5 cm. chart 
width: 0.5 mm. 


Rass INDUSTRIAL DIVISION 


SANBORN CO. 


Preamplifier (A) 195 Massachusetts Ave., Cambridge Mass. 


with Power Supply (B) provides complete de- 
which are normally + scriptive data, specifi- 


cations and prices on 
all 150 Systems, Acces- 
sories and Cabinets. A 


already in place in the 
Basic Cabinet Assembly tH 


2 * copy will be sent on 
request, 
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December 14 at the National Aviation 
Club. Andrew G. Haley introduced 1955 
National President Richard W. Porter of 
General Electric Co. His topic: ‘The 
Artificial Satellite and the ARS.’’ 

Officers announced for 1956: president, 
Jack Gilchrist (Aerojet-General);  vice- 
president, Erik Bergaust (Aero Digest); 
secretary, Mason Comer (NRL); and 
treasurer, Richard Snodgrass (NRL). 

New Mexico-West Texas. Subject of a 
talk by John W. Townsend, Jr., at the 
December 8 meeting was “The Aerobee 
Rocket and Upper Atmosphere Re- 
search.’’ Dr. Townsend is a physicist at 
the Rocket-Sonde Research Branch of the 
NRL, and is head of the Aerobee-Hi devel- 
opment program. He described the de- 
velopment of the Aerobee from the original 
or standard versions to the new Aerobee- 
Hi which is expected to-carry a payload of 
150 Ib to a height of 180 miles (Navy 
version) or to 135 miles (Air Force ver- 
sion). 

At the annual meeting the following offi- 
cers were elected: 

President, Russell K. Sherburne (re- 
elected); vice-president, G. Harry Stine; 
secretary, Keith E. Hennigh; treasurer, 
Lawrence W. Gerdenhire; Board of Direc- 
tors: Nathan Wagner, Levering Smith, 
Dudley M. Gottler, Gilbert L. Meredith, 
William L. Hancok, Allen Niles, Herbert 
L. Karsch, and Henry F. Flamm. 

New York. In the December elections, 
Directors chosen to serve three years are: 
John D. Clark, head of the Chemistry 
Branch, NARTS; Robert A. Gross, chief 
research engineer, Fairchild Engine Div.; 


and Norman V. Peterson, Missile Section, 
Sperry Gyroscope Co. 

North Texas. Walter Dornberger, 
guided missile specialist for Bell Aircraft 
at Buffalo since 1950, addressed the newly 
formed Section at Arlington State College 
on December 7. An overflow audience of 
more than 200 heard him speak on con- 
temporary rockets as well as on the future 
of rocketry. 

Pacific Northwest. Officers elected for 
1956 are: president, Jim C. Drury 
(Boeing); vice-president, James H. Fisher 
(Univ. of Washington); secretary, Joseph 
A. Brousseau (Boeing); and (reasurer, 
George W. Hettrick (Boeing). Board of 
Directors: G. Truxton Ringe; William L. 
Gray; William G. Ramroth; W. Emmett 
Coon; Robert Plath; Henry L. Turner, 
and Robert M. Bridgforth. Al! Directors 
are with Boeing. 

Southern California. [Elected to the 
1956 Board of Directors are: president, 
Richard D. Geckler (Aerojet); vice-presi- 
dent, James A. Broadston (Rocketdyne); 
secretary, William J. Flaherty (Hughes); 
treasurer, Robert G. MacEwan (Douglas). 
Directors (2 years): C. E. Bartley (Grand 
Central); W. E. Campbell (Aerojet); and 
R. B. Canright (Rocketdyne). For 1 year: 
T. F. Dixon (Rocketdyne); R. M. 
Gompertz (ERTS); W. A. Hawkins 
(Lockheed); and W. C. House (Aerojet). 


Indiana Hears Stehling 


URT R. STEHLING of Bell Air- 
craft Corporation spoke to approxi- 
mately 95 persons at a meeting of the 


Indiana Section held at Purdue Univer- 
sity, November 17. 

Mr. Stehling discussed the chain of 
events leading to the announcement of 
the satellite program. Slides were shown 
illustrating the various methods that have 
been proposed for launching a satellite 
vehicle as well as some of the suggested 
designs of the satellite. 

Movies were shown of the Bell X-2 
experimental rocket plane and the Bell 
verticle take-off jet plane. 


Student Groups Forming 
in New York 


(THERE are three student groups in the 

process of formation within the New 
York Section. New York University: 
Contact is Ned Greenburg, New York 
University, College of Engineering, Bronx 
53, N. Y. Academy of Aeronautics: Con- 
tact is Ralph Hautau, Instructor, Academy 
of Aeronautics, La Guardia Airport, 
Flushing 51, N. Y. Brooklyn Polytechnic 
Institute: Contacts are Don Morin, 156 
Dean Street, Brooklyn 17, N. Y., and 
Mario Cordullo, 651 De Graw Street, 
Brooklyn, N. Y. 


Sacramento Section 
Being Initiated 


Ps are being made to initiate a 

Sacramento Section of ARS. Mem- 
bers in that region can obtain further in- 
formation by getting in touch with Frank 
A. Coss at 2605 Midland Way, Carmichael, 
Calif. 


AIMS AND PURPOSES: 


detail the careers they offer. 


—combustion 
—guidance 

—rocket development 
—instrumentation 
—telemetry 
—ballistics 


Announcing a Special Supplement on 


to be published with the 


May 1956 issue of Jer PRoPULSION 


—heat transfer 

—high temperature materials 
—operations research 
—aerodynamics 

—structures 

—electronics 

—computers 


Careers in the Rocket and Missile Industry 


@ To offer companies, as v. ell as government and university research and testing centers, an opportunity to describe in 


@ Toreach an audience which includes, in addition to ARS members, students who are interested in, and potentially 
qualified for, the rocket and missile field. Distribution to students will be made through the cooperation of their professors. 
@ Toinclude the facts on careers in such areas as: 


—propellants 

—missile range testing 
—missile tracking 

—upper atmosphere research 
—controls 

—boundary layer 

—vibration and shock 


@ To help the student relate the work he is doing in the university to its ultimate application in industry. 
@ To give him a perspective view of the entire rocket and missile field. 


Invitations to reserve space (at.regular rates for JET PROPULSION advertising) and submit copy for this section have already 
been mailed. If your organization has not received one, and wishes to participate, please contact Jet Proputsion, 500 Fifth 
Avenue, New York 36, N. Y., or one of the advertising representatives listed on page 65. Deadline for material is April 2. 
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Newell Speaks on Upper 
Air Research 


OMER E. NEWELL, Jr., discussed 
“Upper Air Research With Rockets” 
at the November 17 Dinner Meeting 
of the Cleveland-Akron Section in Akron. 
Dr. Newell is head of the Rocket- 
Sonde Branch, Naval Research Labora- 
tory, and has played a leading role in our 
country’s upper air research program. He 
related past accomplishments in the iono- 
sphere, solar physics, cosmic rays, and 
measurements of pressure, temperature, 
and density. 
The meeting was preceded by a tour 
of the Goodyear Aircraft facilities and a 
soci:l hour sponsored by Goodyear. 


Rocket Films Shown for 
New York Section 


A FILM program on applications of 

rocket power was featured at the 
New York Section’s November 18 meet- 
ing, held at the Engineering Societies 
Building. Three films were shown—a 
captured German film on rocket flight, 
launching and test facilities at NAMTC 
Point Mugu, Calif. (introduced by Paul 
Terlizzi of NARTS, Lake Denmark, 
N. J.), and a film describing Reaction 
Motors’ application of rocket power to 
helicopters (rockets-on-rotor). Joe Witko 
of RMI introduced the ROR film. 


Activities Carnival and 
Meeting at Purdue 


MAN Y spectators were attracted to 

the display presented by the Indiana 
Section at the Activities Carnival held at 
Purdue University on October 4. The 
annual event is held to acquaint students 
and faculty with various student organi- 
zations active on the campus. 

Exhibits included a continuous strip 
movie showing the firing of an experimen- 
tal rocket motor at the Purdue University 
Jet Propulsion Center. Accompanying 
the movie was a tape recording of a test 
firing. Also of interest to the public was 
the Purdue Rocket Laboratory Demon- 
stration Rocket Engine. Other items on 
display included an Aerojet ALDW-1500 
Rocket Assisted Take-Off Unit, a Walter 
“Cold” Rato Unit manufactured by 
Walter Motor Works, Kiel, Germany, 
and a cutaway view of the motor from 
an Aerojet Jato unit. 

Approximately 75 persons attended a 
meeting on October 25. P. M. Diamond, 
Section president, introduced the guest 
speaker, M. J. Zucrow, Professor of Gas 
Turbines and Jet Propulsion at Purdue 
and Director of the Purdue University 
Jet Propulsion Center. Dr. Zucrow spoke 
on “Research Problems in Rocketry.” 
Following the talk, two North American 
Aviation films were shown. 

C. F. Warner, Professor of Mechanical 
Engineering at Purdue University, re- 
cently spoke to the Indiana Section on 
“The Fundamentals of Rocketry.” The 
lecture was the first of a series of space 
flight seminars presented for the benefit 
of the members of the Section. Dr. 
Warner discussed some of the important 
parameters of both solid and liquid rockets. 
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Engineers 
whose minds 
travel 7 miles 


find a rewarding future in rocketry at 
REACTION MOTORS 


Here are futures not bounded by the narrow limits of specialization 
or regimented production. At Reaction Motors, YOU are the engineer. 


The problems provoked by man’s determination to conquer time and 
space, to establish service satellites, to master the force of gravity 

and to develop precision in rockets will challenge your creative abilities. 
And the solutions add immeasurably to your professional stature. 


Here at RMI you have at your disposal the tremendous research facilities 
and the accumulated technical skill gathered through 14 pioneering years 
of rocket research and development. All the important speed and 

altitude records for planes and single stage rockets are held by aircraft 

or missiles powered by engines developed by Reaction Motors. 


If you ARE an engineer whose mind travels 7 miles per second, you'll 

reap a richer harvest at RMI. And you’ll enjoy working in a completely 
equipped, brand_new four million dollar plant in the heart of New Jersey’s 
lake and mountain resort area, just 50 minutes from New York City. 


We Have immediate Openings 


For Men With Experience or Training 


in The Following Fields: 


@ Stress Analysts—Aircraft © Solid Propellant Grain 


@ Formulation of Control 


Engines or Airframes Design for Seat Catapult System Problems, 

e@ Environmental Testing Ejection Systems Computer Set-Up and 
of Components © Thrust Chamber and Solution Analysis 

@ Aircraft Power Plant Injector Design and @ Rotating Machinery — 
Installation Combustion Stability Design and Consultation 

@ General Project © Design and Development @ Rocket Engine Testing 
Engineering Work of Hydraulic, Pneumatic and Instrumentation 

e Engine Controls and and Mechanical @ Internal Combustion 
Controls Systems Components and for Jet or Rocket 

© Design of Airframe System Synthesis Engine Applications 


or Aircraft Engine- 
Type Parts 


@ Servo-Mechanisms 


REACTION MOTORS, INC. 


“First in Rocket Power” 
FORD ROAD 


@ Analysis and Consultation @ Preliminary Design or 
on Rocket Engines and 
Power Plant Systems 


Proposal Experience in 
Aircraft Power Plant Field 


Send Complete 
Resume Including 
Experience & Salary 
Requirements 

to Employment 
Manager 


DENVILLE, N.J. 
AFFILIATED WITH THE OLIN MATHIESON CHEMICAL CORP. 


*Earth’s Escape Velocity 
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positions now open 
at all levels 
(Do not involve board work) 


Aerodynamic Design and Analysis 
Mechanical Development and Design 
Controls and Accessories 

Product Evaluation 

Test Facilities Design 

Manufacturing Engineers 


Preliminary Design 


Please Write To: 

MR. MARK PETERS 
Technical Recruiting, Bldg. 100 
Aircraft Gas Turbine Division 


ocr 


Work where 
making engineering history 
is a habit... 


and creative thinking 
is the way of doing it 


-When General Electric developed America’s first jet engine, it was only 


the beginning of a long line of engineering “firsts” (and “mosts”) destined 
to come from the department now known as the Aircraft Gas Turbine 
Division. 

From that start came eventually the nation’s most-produced jet engine, 
the J47, the powerplant backbone of the Air Force today, as well as the 
J73, J79 and other notable GE models. Advanced concepts in rocket 
engines are also developed and built in this GE department. Work on 
atomic aircraft engine applications is in progress. 


The moving force behind these accomplishments has been the creative 
thinking of engineers who work in a stimulating atmosphere of freedom 
... who are ready to take on new responsibilities as they arise ... who know 
the jet engine field from every aspect because they’ve worked on every 
aspect...and who have confidence in what the future holds for them at GE. 


There is room — and opportunity —for more men like them to take the jet 
field into more advanced stages and to advance their own careers as they 
do so. As much as opportunity, these openings also offer the outstanding 
advantages of association with GE—in terms of employee benefits, 
educational opportunities, stability and exceptional working conditions. 


If this is the kind of environment, the kind of opportunity that interests 
you, why not learn more about these openings? 


GENERAL ELECTRIC 


CINCINNATI 15, OHIO 
JET PROPULSION 
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New Equipment and Processes 


Equipment 


Electrical, Electronic 


Ultrasonic Flowmeter. Contains two 
transducers which measure mass or volu- 
metric flow at almost zero pressure drop. 
Fischer & Porter Co., Hatboro 35, Pa. 


Dust & Explosion Proof Panelboards. 
Protects against vapors and dusts for 120 
volt service. Crouse-Hinds Co., Wolf & 
Seventh St., Syracuse 1, N. Y 


Small Selenium Diodes. Types 6U1, 
7U1, and 8U1 are for bias service at temps 
of —80 to +100 C. Pigtail leads. Inter- 
national Rectifier Corp., 1521 E. Grand 
Ave., El Segundo, Calif. 


High-Speed Recorder. Electronic re- 
corder allows 0.25-sec span for 11-in. wide 
chart. Used with DC output transducers. 
Basic chart speeds of 2 in/sec or 4 in/min. 
Minneapolis-Honeywell Regulator Co., 
Industrial Div., Wayne & Windrin Ave., 
Philadelphia 44, Pa. 


8-Channel Oscilloscope. Features 17- 
in. rectangular CRT with 2% display 
accuracy. Each channel has '/2 sec scan. 
A-c or d-c operation. Electromec, Inc., 
Strain Gage Oscilloscope Dept., 5121 San 
Fernando Road, Los Angeles 39, Calif. 


Subminiature Modules. RF, IF, and 
power supply modules weigh about 2 lb 
and measure 1” X 3” X 9”. Available 
frequencies from 40 to 235 me. Land-Air, 
Inc., Instrument & Electronics Div., 
Bldg. 711, Oakland International Airport, 
Oakland 14, Calif. 


Sun Battery Cells. Areas of 0.14-10.5 
sq. in. and 0.1-15 milliwatts output. 
Area of 8 sq. in. can drive 5-stage transistor 
radio. Encapsulated for long life. Inter- 
national Rectifier Corp., Product Informa- 
tion Dept., 1521 E. Grand Ave., El 
Segundo, Calif. 


Ceramic Micro Miniature Electron 
Tube. New tubes are as small as transis- 
tors. First tube available will be 6 by 4. 
Temperature service up to 500 C. 5/1¢-in. 
diam., !/,-in. long. General Electric Co., 
Schenectady, 


Microjet Pneumatic Control. For tur- 
bojet engine afterburner light-off and blow- 
out detection, turbine temperature con- 
trol, overspeed control, and afterburner 
fuel control. Also for ramjet temperature 
control, cruise matching, compressor surge, 
and shock wave position. Solar Aircraft 
Pacific Highway, San Diego 12, 

‘alif. 

Silicon Power Rectifier. A-c to d-c is 
rated at 200 volts peak and 15 amperes, 
operates at —55 to +170 C, has negligible 
voltage drop and is thimble-sized. Ray- 
theon Mfg. Co., Waltham 54, Mass. 
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Shock-Resistant Galvanometer. With- 
stands 25G, compact galvanometer has 
sensitivity up to 0.105 microamp/mm div. 
Completely housed models. Howell In- 
., 1106 Norwood, Fort Worth, 

ex 


Miniature Pressure Transducer. Gages 
absolute, gage, and differential pressures. 
Features low shift and can operate at —65 
to +250 F. Statham Laboratories, Inc., 
a W. Olympic Blvd., Los Angeles 65, 

if. 


Oxygen and Combustibles Measuring. 
Exhaust gases are measured and displayed 
on meters in this portable (18 X 12 x 18 
in.) precision instrument. Instrument 
55 Halleck St., Newark, N 


Germanium Rectifier. A-c to d-c has 
rating of 125 volts and 800 amperes, yet is 
about the size of a dime. Requires water 
cooling for power lines. General Elec- 
tric Co., Schenectady, N. Y. 


Analog-to-Digital Converter. 0.1% ac- 
curacy, resolves 30 mv, converts at rate of 
100,000 per sec. Can be used with trans- 
ducers, tapes, readers, etc. J. B. Rea, Inc., 
rite A Cloverfield Blvd., Santa Monica, 

alif. 


Data Reduction Filter Set. Analyzes 
vibrations in missiles from tape input. 
23 band-pass filters give output to panel 
jacks. Hycor, Inc., 11423 Vanowen St., 
North Hollywood, Calif. 


Bridge Balance. Four channel, Type 
8-110 controls output from transducers to 
recorders, meters, or conversion systems. 
Consolidated Engineering Corp., 300 N. 
Sierra Madre Villa, Pasadena 15, Calif. 


Doppler Data Digitizer. Missile data 
from radar tracker is converted to digital 
form and stored on tape. Potter Instru- 
ment Co., Great Neck, L. I., N. Y. 


Pressurized Connector for Missiles. 
Provides for 40 contacts with fast discon- 
nect, Series GA. DeJur-Amsco Corp., 
83 Northern Blvd., Long Island City 1, 


aN. 


Recording Potentiometer. Weighs less 
than 25 lb for use with iron-constantin, 
copper-constantin, and chromel-alumel 
thermocouples. 5 in. wide chart. Pen 
speed is 1 sec full scale. Accuracy 0.5%. 
Westronics, Inc., 3605 McCart St., Fort 
Worth, Tex. 


Mechanical 


Hi-Pressure Pumps. Plunger-type can 
handle up to 50,000 psi. Flows of 2 cc/hr 
to 15,000 gph. Handles acids, oxidants, 
etc.’ Philadelphia Pump & Machinery Co., 
1513 Race St., Philadelphia, Pa. 


Emergency Field Shower. Self-con- 
tained, mobile shower is trailer mounted. 
—— Emergency Showers, Glendale, 

alif. 


AN Check Valves. Approved '/,” Al 
alloy for 3000 psi. Operating range of 
—65 to +160 F. Can be supplied for 275 
F operation. Aircraft Products Co., 300 
Church Rd., Bridgeport, Pa 


Lightweight Flexible Hose. Stainless 
steel braid with Teflon inner lining resists 
solvents, corrosive fluids from —100 
to +450 F. Titeflex, Inc., 500 Freling- 
huysen Ave., Newark 5, N. J. 


Missile Fuze Tester. Vari-G unit can 
test two 1-lb units up to 300 G or two 10-Ib 
units up to 100 G. Rocket launcher condi- 
tions can be simulated in this cabinet 
centrifuge and console. Hycon Manu- 
facturing Co., 2961 E. Colorado St., 
Pasadena 8, Calif. 


REMOTE ACTUATED COUPLING 


WIGGINS REMOTE COUPLING 
IN COUPLED POSITION 
ON MISSILE 


ELECTRICAL CABLE 


Remote Actuated Coupling. Uses stand- 
ard AN connections. Lanyard-type re- 
lease. E. B. Wiggins Oil Tool Co., 3424 
E. Olympic Blvd., Los Angeles 23, Calif. 


Test 


Test Accelerator. Objects up to 5-in. 
cube (max 3 Ib) can be tested at accelera- 
tions of up to 250G. Controlled speed 
range up to 1000 rpm. Equipped with 
control console. Shaevitz Machine Works, 
Camden 1, N. J. 


Sequence Camera. 70-mm takes photos 
at rates of 5-20 per sec. Weighs 12'/2 lb. 
Exposures of 1/25 to '/2s% sec. Lenses of 
3-80 in focal length. 12 volt, 1 roll gives 
470 pictures, single frame rates up to 5 
per sec. Charles A. Hulcher Co., 40 
Manteo Ave., Hampton, Va. 


Materials 


Metals 


Fine, Flat Wire. Commercially availa- 
ble i in most alloys in sizes down to 0.003 
in. wide X 0.0007 in. thick. Tolerances of 
+0.00001 in. thick and 0.001 in. wide avail- 
able. Hamilton Watch Co., Allied Prod- 
ucts Div., 904 Wheatland Ave., Lancaster, 
Pa. 


Platinum-Rhodium Thermocouple. 
Measures temps to 1800 C. Thermo 
Electric Co., Inc., Saddle River Township, 
Rochelle Park P. '0., N. J. 


Pressure-Sealed Thermocouples. 
Available in iron-constantin, or copper- 
constantin with 1/;” pipe thread connec- 
tion. Conax Corp., 7811 Sheridan Dr., 
Buffalo 21, N. Y. 


Titanium Shear Bolts. Shear strengths 
of 95,000 psi for NAS-333 Hi-Ti 100° flush 
head bolts. Diameters of 0.25 to 0.5 in. 
Standard Pressed Steel Co., Jenkintown, 
Pa. 
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Plastics 

Single Pa e Epoxy Adhesive. Pot 
life of over rp gives up to 3000 psi 
shear strengths on metal-to-metal bonds. 
Miracle Adhesives Corp., Box 466, New 
Philadelphia, Ohio. 

Single Package Thixotropic Adhesive. 
Bondmaster M620 is based on Ciba Aral- 
dite AN-100 epoxy resin. Various curing 
cycles applicable. Rubber & Asbestos 
Corp., 225 Belleville Ave., Bloomfield, N. J. 


Low Viscosity Epoxy. Liquid Epon 815 
is in the 500-900 centipoise range. Makes 

ssible glass fiber contents of 70% in 
He mes tes. Shell Chemical Co., 50 West 
50th St., New York 20, N. Y. 


Foamed-in-Place Insulation. Two- 
package polyurethane foam features resin 
and catalyst. Chase Chemical Corp., 
3527 8 St., Pittsburgh 1, Pa. 


Resistant Coatings. For missile skins, 
test stands has high strengths, chemic 
resistance. Clear or pigmented. Amchem 
Corp., 283 Tuxedo, Hi d Park 3, Mich. 


Product Literature’ 


Working With Stainless Steel. Fabri- 
cator’s 162-page handbook covers cutting, 
heat treating, forming, joining, and ma- 
chining. Crucible Steel Co., Bldg., 
Pittsburgh 30, Pa. 

Solid Propellant Stud Driver. How to 
use powder actuated fasteners. Industrial 
Sales Div., Remington Arms Co., Bridge- 
port 2, Conn. 

Pressure Transducers. Bulletin MPT- 
1, 130, 131, and 132 covers transducers and 
adaptors. Statham Laboratories, Inc., 


1 Unless otherwise specified, literature is 
free and available by request on company 
letterhead. 


12401 W. Olympic Blvd., Los Angeles 64, 
Calif. 


Patents. 19-page booklet tells what an 
invention is and steps in por in | it. CGS 
Laboratories, Inc., 391 Ludlow St., Stam- 
ford, Conn. 


Titanium. 2-page folder gives list of 
available titanium alloys, processing meth- 
ods, prices, etc. Arthur D. Little, Inc., 
30 Memorial Drive, Cambridge 42, Mass. 


Steel Tubing. 8-page catalog lists car- 
bon and alloy steel tubing for mechanical, 
ressure, air-frame, and aircraft use. Ohio 
ess Tube Div., Copperweld Steel 

Co., Shelby, Ohio. 


Radiation Pyrometer. Bulletin 9320 
describes Radiametric detector. Response 
is less than 0.5 sec up to 2600 F. Indus- 
trial Div., Minneapolis-Honeywell Regu- 
lator Corp., Philadelphia 44, Pa. 

Pressure Conversion Chart. Bulletin 
CEC 1557 presents charts and tables for 
conversion. Consolidated Engineering 
Corp., 300 N. Sierra Madre Villa, Pasa- 
dena 15, Calif. 


Missile Vibration Testing. 8-page book 
tells of problems in developing procedures. 
MB Manufacturing Co., Inc., 1060 State 
St., New Haven, Conn. 


Fluocarbon Plastics. 20-page brochure 
Bulletin No. IN-554 on Teflon, Kel-F, 
and Bakelite Fluorothene. United States 
Gasket Co., Box 93, Camden, N. J. 


Ceramics. Bulletin CHE-R5 describes 
porcelain, stoneware, coatings, graphite, 
and alumina for heat and corrosion resist- 
09 General Ceramics Corp., Keasbey, 


Humidity Control for Solid Propellant 
Development. [Illustrated pamphlet K- 
1955-C describes constant humidity con- 
trol equipment for propellant processing of 


hygroscopic materials for Phillips Petro. 
leum Co., Bartlesville, Okla. Thomag 


Engineering Co., Tulsa, Okla. 


Missile Program. Bulletin R55A0519 
reviews ten-year program. General Elec. 
tric Co., Special Defense Projects Div,, 
Schenectady 5, N. Y. 


Moly Stainless Steel. Technical data 
on AM350 stainless steel is given in this 
10-page booklet. Allegheny Ludlum Steel 
ai 2020 Oliver Bldg., Pittsburgh 22, 

‘a. 


Steel Alloys. Brochure describes high 
strength, corrosion, and temperature re- 
sistant alloys for jets, rockets, and missiles, 
'yclops Steel Corp., Bridge 

e, Pa. 


Asbestos Textiles. Cloth, cord, lap, 
rope, roving, tape, tubing, wick, and yarn 
asbestos textiles for high temperatures are 
described in this 16-page manual. Johns- 
—— 22 E. 40th St., New York 16, 


Phenolic Tooling Plastics. Product 
Data No. 107, a 6-page booklet gives prop- 
erties, handling, and applications of R- 
728 phenolic tooling. Rezolin, Inc., 5736 
W. 96th St., Los Angeles 45, Calif. 


Plaskon Polyester Resins. Gives prop- 
erties, applications, of a line of flexible, 
rigid polyesters. Barrett Div., Allied 
Chemical & Dye Corp., 40 Rector St., 
New York 6, N. Y. 


Grinders. Catalog of a line of Pulva- 
Sizers, feeders, and auxili: equipment. 
— Corp., 550 High St., Perth Amboy, 


Vacuum Metallizing. Catalog 551 de- 
scribes vacuum metallizing process and de- 
tails experimental and production units. 
High Vacuum Equipment Corp., 349 
Lincoln St., Hingham, Mass. 


SCINFLEX WATERPROOF PLUGS 


FOR USE WITH MULTI-CONDUCTOR CABLES! 


THESE BUILT-IN FEATURES ASSURE TOP PROTECTION 
AGAINST CIRCUIT FAILURE: 
Shock and Vibration Resistant * Die Cast Alu- 


These new Bendix*-Scinflex waterproof plugs 
are a modification of our standard AN type 
“E” (environment resistant) connector. They 
are designed to meet all “E” performance re- 
quirements when used with multi-conductor 
cables. Each plug includes a modified AN- 
3057B cable clamp which provides inward 
radial compression on multi-conductor cables. 
This unique feature completely eliminates 
cable strain—a common source of circuit 
trouble. 
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In addition, there are gaskets at all mating 
surfaces and an accessory sleeve is available 
to accommodate an extreme range of cable 
sizes. A descriptive folder is available on 
request. 


*TRADE-MARK 


minum Shell * Cadmium Plate—Olive Drab 
Finish Moisture-Proof, Pressurized High 
Arc Resistance, High Dielectric Strength ° 

Silver-Plated Contacts * Resilient Inserts 


SCINTILLA DIVISION OF 


SIDNEY, NEW YORK 


AVIATION CORPORATION 


Export Sales and Service: Bendix International Division, 205 E. 42nd St., New York 17, N. Y. 
Factory Branch Offices: 117 E. Providencia Avenue, Burbank, California * Brouwer Building, 176 W. Wisconsin 
Avenue, Milwaukee, Wisconsin * Stephenson Building, 6560 Cass Avenue, Detroit 2, Michigan © 512 West 
Avenue, Jenkintown, Pennsylvania * 8401 Cedar Springs Rd., Dallas 19, Texas * American Building, 4 South 
Main Street, Dayton 2, Ohio 
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ING FOR THE 


... taxes your knowledge of the many com- 
plexities of rocket fueling — energy, weight, 
supply, cost, conductivity, ignition qualities. 


It will pay you to consider the characteristics 
of Pennsalt Fluorine and Pennsalt Chlorine 
Trifluoride. Both are now available in un- 
limited quantities. 


Whether yours is a research or a purchasing 
job, Pennsalt’s wide experience in developing 
and supplying interhalogen compounds and 
fluorinated oxidants will be of value to you. 
Plan now to discuss your liquid rocket-pro- 
pellant problems with a Pennsalt representa- 
tive. Call or write Technical Services Depart- 
ment 264, Industrial Chemicals Division, 
Pennsylvania Salt Manufacturing Company, 
Three Penn Center Plaza, Philadelphia 2, Pa. 


Pennsalt 


Chemicals 


PENNSYLVANIA SALT MANUFACTURING COMPANY 
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New Patents 


Fuel scheduling control system for gas 
turbine engines (2,720,751). Wm. J. 
Kung, Jr., South Bend, Ind., assignor 
to Bendix Aviation Corp. 

An electric servomotor for varying rate 
of fuel feed to a burner. Means for 
automatically limiting the magnitude of 
the deceleration signal as a function of 
density of air flowing to the engine. 


Gas turbine apparatus (2,724,546). Law- 
rence D. Barrett, Jr., Glenn Mills, and 
John yi Prospect Park, Pa., as- 
signors to Westinghouse Electric Corp. 

Casing structure with longitudinal 
guides and a series of integral annular 
diaphragms with recesses for slidably en- 
gaging the guides. Diaphragms carry 
radial stationary vanes in stacked relation 
within the shell and between the end abut- 
ment members mounted in opposite ends 
of the shell. 


Automatic speed control for aircraft 
(2,724,564). John C. Newton, Roslyn 
Heights, and Max Irwin Beers, Great 
—— N. Y., assignors to Sperry Rand 


rp. 

Air speed reference device and baro- 
metric altitude measuring device which 
limit the maximum air speed to a pre- 
determined Mach number by control of 
the supply of fuel to the engine. 


Pressure differential control for variable 
jet nozzle (2,724,947). Chas. A. Meyer, 

rexel Hill, Pa., assignor to the U. S. 
Navy. 

Electro-responsive variable area jet 
nozzle and afterburner, operated inter- 
mittently to agument the thrust at the 
nozzle. Switches open or close operating 
circuits in response to sudden variations 
in pressure incident to ignition or extinc- 
tion of gases in the afterburner. 


High-speed landcraft (2,724,966). John 


K. Northrop, Claude A. Brosterhous, 
Leigh E. Dunn, John M. Ellison and Dick 
R. Herman, Los Angeles, Calif., assignors 
to Northrop Aircraft, Inc. 

Sled for testing airfoils at transonic 
speeds in open air propulsion. Legs 
extending from the body connect to slip- 

rs mounted to slide on track rails. A 

m extends forward to hold the airfoil 
being tested, and a rocket motor at the 
rear supplies the motive power. 


Resilient mounting arrangement in gas- 
turbine engine installations (2,724,948). 
Geo. H. Hiscock and Robert Garbutt, 
Bristol, England, assignors to Bristol 
Aeroplane Co. 


A compartment surrounding the com- 
pressor section, and a second surroundin 
the combustion equipment, separa 
by a wall to prevent the spread of fire 
between compartments. Resilient means 
connect a mounting unit within the first 
compartment, and air is passed to the 
second compartment for ventilation. 


2,726,507 


Jet power plant controls (2,726,507). 
ohn H. Baker, Washington, D. C., and 
Wendell E. Reed, Chula Vista, Calif., 
assignors to Solar Aircraft Co. 

Actuating means within the tailpi 
for applying vernier adjustment to the 
nozzle in either of its basic positions, by 
varying the area of the nozzle in response 
to temperature varitions. 


2,726,508 


Aircraft propulsion units of the gas- 
turbine jet type (2,726,508). F. B. 
Halford, E.S. Moult, and Wm. H. Arscott, 
H nden, England, assignors to The 
DeHavilland Engine Co., Ltd. 

Impeller in the passage in advance of 
the compressor and driven in the same 
direction, but at a lower speed than the 
centrifugal impeller. The impeller has 
vanes of greater diameter than the entry 
eye of the compressor, and the entry pass- 
age diameter progressively diminishes from 
the impeller vanes to the entry eye of 
the compressor. 


Head shield for airplane occupants 
(2,726,054). Robert F. Lesley, Cuyahoga 
Falls, Ohio, and Robert C. Tomkins, 
Los Angeles, Calif., assignors to Goodyear 
Aircraft Corp. 

A fixed head-rest box, open at front and 
bottom, for covering the rear half of the 
head and restraining the head from lateral 
movement. A rigid one-piece shield forms 
an open sector pivotally mounted on the 
box for forward movement. When ejection 
takes place, the face shield is tilted against 
the chest, completely covering the head. 


Method of operating an aircraft propulsion 
system (2,725,108). Robert A. Grossel- 
finger, Washington, D. C., assignor to 
Grovar, Inc. 

Steps for increasing fuel flow to the 
turbine and simultaneously increasing the 
rotor pitch. These steps decrease turbine 
speed simultaneously with power increase, 
and the kinetic energy stored in the 
turbine-rotor system is released for 
propulsion of the aircraft. 


George F. McLaughlin, Contributor 


Ignition system for jet motors and the like 
(2,725,718). Herman E. Sheets, Howard 
R. Hegbar, John C. Feldscher, and Carl A. 
Meneley, Cuyahoga Falls, Ohio, Assignors 
to Goodyear Aircraft Corp. 

A pulsating oscillator of high frequency 
to establish corona discharge between a 
flame arrester and an electrode. [lec- 
trical means connecting low-frequency 
power in arcing relation with the flame 
arrester in the region of the corona dis- 
a between the arrester and the elec- 
trode. 


Fuel system for gas-turbine and 
the like (2,725,932). David H. Ballan- 
tyne and Albert Jubb, Buttershaw, 
Bradford, England, assignors to Rolls 
Royce, Ltd. 

Centrifugal pump delivering engine 
fuel, and a fixed capacity positive-dis- 
placement pump supplying fuel require- 
ments at low rotational speeds delivering 
fuel to the fuel injector in parallel with 
the centrifugal pump. 


2,726,510 
Flight-control apparatus involving steering 
combustion chambers (2,726,510). Robert 
H. Goddard (deceased), by Esther C. 
Goddard, executrix, Worcester, Mass., 
assignor of one half to the The Daniel 
and Florence Guggenheim Foundation. 

A pair of combustion chambers and 
discharge nozzles at the front of the casing, 
and separate means for supplying two 
liquid combustion agents to both cham- 
bers. Valves vary the relative flow of 
both liquids and the relative thrust of each 
combustion chamber. 


= 


Afterburners (2,726,511). Paul A. Pitt, 
San Diego, Calif., assignor to Solar 
Aircraft Co. 

A first grid of fuel distribution pipes 
and a second grid having a pair of inter- 
connected carburetion and combustion 
gutters. An intermediate interconnected 
combustion gutter is attached down- 
stream to the second grid. 


Catapult for aircraft personnel and other 
uses (2,726,055). Walton Musser, 
Philadelphia, Pa., assignor to the U. S§. 
Army. 

An explosive charge within gunlike 
tubes furnishing a force for propelling the 
inner tube and attached seat out of the 
aircraft. Axial movement of remaining 
tubes is prevented by a lock which may be 
released to permit axial movement with 
respect to the outermost tube. 


Eprror’s Nore: The patents listed above were selected from recent issues of the Official Gazette of the U. S. Patent Office. 
Printed copies of patents may be obtained at a cost of 25 cents each, from the Commissioner of Patents, Washington 25, D. C. 


126 


JET PROPULSION 


| 
mer | | 
274947 * 
N | 
4 
| 
4 
] 


haw, 
Rolls 


igine 
uire- 


‘ring 
with 


Book Reviews 


Servomechanisms and Regulating Sys- 
tem Design, vol. II, by Harold Chestnut 
and Robert W. Mayer, John Wiley & 
Sons, New York, 1955. $8.50. 

Reviewed by A. RosENBLOOM 
Ramo-Wooldridge Corp. 


The authors state in their preface that 
this book was written for the practicing 
servomechanisms designer and the ad- 
vanced graduate student. This reviewer 
only half concurs. The book contains a 
wealth of practical material, but the level 
of presentation limits its usefulness as an 
advanced textbook in servomechanisms, 
as does the absence of problems. 

The range of topics covered is diverse, 
ranging from a detailed chapter on d-c 
amplifier design to an inadequate section 
on random noise. Analysis of nonlinear 
control systems is presented using the 
describing function method and small 
signal perturbation theory. No phase 
plane techniques are discussed. There 
are many examples worked out in detail. 

There are nine major sections: 1. 
Measurement techniques; 2. Influence 
of input characteristics on control system 
design; 3. Selection of controlled sys- 
tem power element with proper rating; 
4. Networks for obtaining desired at- 
tenuation-frequency characteristics; 5. 
Amplifier design; 6. All alternating- 
current servomechanism operation; 7. 
Linearization of nonlinear elements for 
small departures; 8. Linearization of 
nonlinear elements for large departures; 
9. Application of nonlinear elements to 
control systems. However, no mention is 
made of such important topics as root 
locus techniques, Wiener filters, modern 
techniques of network synthesis, or the 
effect of noise in nonlinear servos. 

The style of this book is somewhat 
verbose, and the presentation is occa- 
sionally marred by loose statements, such 
as those in Section 7.4, where the impres- 
sion is given that the presence of time 
variable parameters in a system implies 
that the system is nonlinear. 

This book is of value to those concerned 
with the detailed design of control 
systems. 


Introduction to Nuclear Engineering, by 
Raymond L. Murray, Prentice-Hall, 
Inc., New York, 1955, 418 pp. $9.35. 

Reviewed by A. R. Crocker 
Aircraft Nuclear Propulsion Dept. 
General Electric Co. 


This book comes closer to being a real 
introduction to the new and fascinating 
field of nuclear engineering than the few 
books now published on this subject. It 
is intended as a text for junior and senior 
year students. There are twenty chap- 
ters which may be arbitrarily grouped in 
the following topical categories: Basic 
Information—Atomic and nuclear physics; 
Neutrons; Fission and the chain reaction; 
Nuclear reactor principles; Shielding; 
Neutron experiments. Reactor Design— 
The water boiler; Design of gas-cooled 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


enriched uranium reactor; Design of a 
liquid metal natural uranium reactor. 
Applied Problems—Reactor start up and 
operation; Principles of heat transfer and 
fluid flow; Radioactive waste disposal; 
Detectors and control instruments; Sep- 
aration of isotopes; Plutonium produc- 
tion; Radiation hazards. Application— 
Uses of radioactive and stable isotopes; 
Nuclear propulsion of aircraft submarines 
and rockets; Electrical power from 
fissionable material. 

The only apparent omission in subject 
matter which should have been included 
to cover the field completely is that of 
remote systems engineering; i.e., the 
application of special apparatus, optical 
and mechanical, for the handling of 
highly radioactive materials while per- 
forming chemical, metallurgical and me- 
chanical operations and examinations. 

One of the very useful features of the 
book are the two chapters which give the 
procedure and calculations for the nuclear 
design of two types of reactors. 

Of interest to jet propulsion engineers 
are the sixteen pages devoted to the 
nuclear propulsion of aircraft and rockets 
which provide some unclassified infor- 
mation which can serve as a starting point 
for consideration of the particular prob- 
lems associated with nuclear propulsion. 

Nuclear reactor theory is always a 
stumbling block for the average engineer 
or engineering student. In this book, 
one short chapter provides an introduction 
with a descriptive explanation of the terms 
and principles with a minimum of mathe- 
matics. It should provide a picture of 
the processes such that it will evoke 
interest in further study rather than dis- 
heartening the student at the outset with 
“cold” mathematics. 

Since there are an average of 19 pages 
per chapter, the book cannot exhaustively 
treat any subject (a book could probably 
be written on each chapter), but it does 
provide a well-balanced introduction to 
nuclear engineering for the average 
engineer or engineering student. 

While the technician and amateur may 
have difficulty with the mathematical 
treatment, they should nevertheless find 
the discussions and applications of con- 
siderable value to them. Chapters 5 
through 7, covering amplifiers, oscillators, 
and other applications, will give the tech- 
nician and amateur a good working knowl- 
edge of transistors. 


In my opinion, this book is best suited | 


to the practicing engineer who has no 
transistor background as stated above. 
A bibliography would be a helpful and 
almost necessary addition to this book. 


Thermodynamics, by J. F. Lee and F. W. 
Sears, Addison-Wesley Publishing Co. 
Inc., Cambridge, Mass., 1955, 543 pp., 
$7.50 

Reviewed by C. F. WARNER 
Purdue University 


The two authors, Lee and Sears, emi- 
nently successful authors in their separate 


fields of engineering and physics, have 
combined their talents to produce a most 
worth-while contribution to the thermo- 
dynamics literature. The authors have 
skillfully blended the basic physical ther- 
modynamics and the applied thermody- 
namics of engineering to form a coherent 
treatise. The blending of the classical and 
the applied thermodynamics is evident 
throughout the book and is not limited to 
the so-called basic chapters as is the case in 
some previously published textbooks in 
this field. 

As stated in the preface, “The purpose of 
the book is to provide a mature approach 
to the basic principles of thermodynamics 
for all engineering and applied science stu- 
dents regardless of their ultimate profes- 
sional objectives. Emphasis is placed 
on basic principles, rather than on de- 
tailed engineering applications involving 
descriptions of equipment and design 
techniques, in order to encourage a fuller 
understanding of thermodynamics as a 
basic science. ..” It is the opinion of the 
reviewer that the authors have accom- 
plished their stated purpose; however, the 
book may find wider acceptance as a 
graduate level text and a reference book 
than as an undergraduate text. 

Extensive use is made of differential 
equations in the presentation of the phys- 
ical properties of matter, so that the stu- 
dent acquires the concept of the univer- 
sality of the physical property function and 
is thus able to deal with substances whose 
properties are not completely known. 
Numerous examples of the experimental 
measurement of physical properties are 
presented to augment the mathematical 
treatment. The many three-dimensional 
figures of the p-v-t surfaces for real sub- 
stances aid in explaining the physical 
significance of some of the partial differen- 
tial coefficients, and in showing the differ- 
ence in the variation of these properties 
for different substances. 

Emphasis is placed upon the selection 
of a system and then its surroundings as 
the first step to be taken in the analysis 
of- any thermodynamic _ problems. 
Throughout the book various systems 
of an electrical, chemical, mechanical, 
etc., nature are used to illustrate the uni- 
versality of this approach. 

The second law presentation is one of 
the outstanding features of the book. 
The use of heat reservoirs or heat storage 
systems in the development of the second 
law concept is common practice. In 
addition to using heat reservoirs, the au- 
thors introduce work reservoirs in explain- 
ing the difference between heat and work 
flow processes. For example, a block 
attached to a spring may be considered 
as two work reservoirs. When the system 
is set in motion on a level frictionless sur- 
face, there is a continual interchange of 
energy between the two reservoirs; how- 
ever, the total energy remains constant. 
At the extreme limits of the motion the 
energy is all potential and at the mid- 
point of the motion it is all kinetic. The 
two reservoirs, the spring and the block, 
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are characterized by the fact that energy 
interchange between them can take place 
in either direction with no change in the 
amount of energy available for doing work 
on another system. Consider two heat 
reservoirs of the same material but at 


different temperatures. When the two 
reservoirs are brought in contact, heat 
will flow from the one at the higher tem- 
perature to the one at the lower tempera- 
ture and the energy flow out of one is 
equal to the energy flow into the other. 
No “oscillations” of temperature are ever 
observed; that is, the temperature of the 
reservoirs do not “overshoot” the equilib- 
rium temperature, with a subsequent re- 
versal of the direction of heat flow. This 
example illustrates the fundamental dif- 
ference between “heat flow” and “work 
flow” and the one-directional nature of 
heat flow. Additional use is made of 
the work reservoir concept throughout the 
discussion of the second law and entropy. 

Fluid mechanics is combined with ther- 
modynamics in treating steady-state fluid 
flow. The major portion of the material 
presented is for compressible fluids and 
extensive use is made of Mach number and 
the common fluid dynamic equations in- 
volving Mach number. 

In addition to the subjects previously 
discussed, chapters dealing with heat 
transfer, single-phase and two-phase mix- 
tures, power cycles, and thermodynamics 
of reactive systems (combustion) are pre- 
sented by the authors. Of these, the one 
dealing with power cycles is the most ex- 
tensive. 

In the opinion of the reviewer, some defi- 
ciencies of the book are the omission of the 


consideration of the nonsteady-state proc- 
esses, the limited discussion of combustion 
processes, and the limitation of the jet 
propulsion discussion to the turbojet en- 
gine. Notwithstanding these criticisms, 
this is one of the best thermodynamics 
textbooks to appear in some time. 


Fluid Dynamics of Jets, by Shih-I Pai, 
D. Van Nostrand Co., Inc., New York, 
1954, 227 pp. $5. 

.Reviewed by J. S. IsenBere 


\\ ? Bell Aircraft Corporation 


In the words of the author, ‘““The pur- 
pose of this book is to present the experi- 
mental and theoretical facts /about the 


Flow of a Jet of Viscous Fluid, Turbulent 
Flow of a Jet of Incompressible Fluid, 
Turbulent Flow of a Jet of Compressible 
Fluid, Jet Mixing of Gases of Different 
Kinds, The Stability of a Jet of Incompres- 
sible Fluid, Stability of a Jet of Compres- 
sible Fluid, Fundamental Equations of 
Fluid Dynamics. 

Unfortunately, no attempt has been 
made to consider some important practical 
problems such as the impingement of jets 
on plates, the jet syphon, ejectors, etc., 
and the discussion of a topic such as thrust 
augmentation is very inadequate. 

Most chapters start with a combined 
general discussion of the subject and his- 
torical sketch of the main contributions. 


steady flow of a jet of gas issuing from a )_/ The bulk of the material is concerned with 


nozzle into the surrounding stream at/ 
rest or in uniform motion.”/ It is the 
opinion of the reviewer thaf the stated 
purpose has by and large been achieved, 
although an improvement in the treat- 
ment of the subject matter and an exten- 
sion of its scope would make a much more 
readable and useful volume. This book is. 


recommended /primarily for those serious \\/ 


students of fluid mechanics who already | 
have a grounding in the fundament 

of viscous and of incompressible flow along 
with an understanding of the pertinent 
mathematical techniques. Engineers look- 
ing for quick solutions to their problems 
will not find them here. 

The logical organization of the book and 
its scope is shown by the following list of 
chapter headings: Flow of a Jet of In- 
viscid and Incompressible Fluid, Subsonic 
Flow of a Jet of Inviscid Fluid, Supersonic 
Flow of a Jet of Inviscid Fluid, Laminar 


the mathematical analysis of idealized jet 
flow problems, and the author has admira- 
bly tried to give the essence of the method 
without becoming involved in details. A 
limited number of experimental results are 
quoted. 


Book Notices 


Thermal Conductivity of Gases, 
‘Measurement of the Thermal Conduc- 
tivity of Gases at High Temperatures,” 
Report PB 111571, $1. OTS, Ss. Dent 
of Commerce, Washington 25, D 


Conversion Factors and Tables, by 
O. T. Zimmerman and Irvin Lavine, 
ee Research Service, Inc., Dover, 

N. H., 2nd ed., 1955, 501 pp., $5. This 
pocket-sized book contains an alphabetical 

ting of U. S. and metric units with 
conversion factors. Tables of conversion 
factors include data on pressure, viscosity, 


COMBUSTION 
RESEARCH 


Engineers needed with up to ten years’ experience for re- 
search in high speed combustion, fluid dynamics pertaining 
to jet propulsion, and in combustion in high compression ratio 
reciprocating engines. Work broadly covers the field from 
fundamental research to engine design. Opportunity for a 
career with a leading petroleum research company. 


Give full and specific details of education, desired salary, 
availability date, and references. All! inquiries will be con- 
sidered promptly and held confidential. 


(Formerly Standard Oil Development Co.) 


Personnel Division 


P. O. Box 51 
Linden, N. J. 


energy, hardness, wire gages, etc. 
design engineer 


ME OR AE 


Power Plants for 
Nuclear Engines and Reactors 


A career opportunity in the newest 
phase of aeronautical advance- 
ment is now open with this major 
company, located in the Midwest 
and West. It is a position which 
requires excellent judgment, sound 
engineering fundamentals, and 
the ability to work with others. 
In turn, it offers the chance for 
rewarding professional experience, 
an unusually fine working environ- 
ment, and a wide range of personal 
benefits. 


The engineer who qualifies for this 
job will be responsible for the 
design and proper functioning of 
the mechanical portions of power 
plant equipment associated with 
both engines and reactors. This 
means he must have 3 to 8 years’ 
experience in the mechanical 
design of light weight, high tem- 
perature machinery, preferably 
- experience in turbojet engines. 
An advanced degree is desirable. 


Address replies to: 


Box H 
American Rocket Society 
500 Fifth Ave. 

New York 36, N. Y. 
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America’s first successful long-range pilotless 
bomber—the TM-61 Martin Matador—is a major 
weapons system conceived and developed in peace- 
time. And it is helping to keep the peace in an im- 
portant way. 

For the Matador is in quantity production, and 
already widely deployed. As the first aircraft having 
complete interchangeability of parts, it is transport- 
able by air and can be delivered unassembled for 
stockpiling in critical outpost areas. 

The new engineering concept that produced the 


THE FORTUNES OF PEACE 


Matador has already revised many design and pro- 
duction standards in the aircraft industry. For this 
versatile weapons system is being produced at the 
lowest cost-per-pound for comparable production, de- 
spite performance requirements more severe than 
those for most piloted aircraft. 

Today, the new Martin engineering concept is re- 
vising the calendar on some of the most advanced 
flight and weapons systems projects now in the re- 
search and development stages. 

These are among the Fortunes of Peace! 
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Technical Literature Digest 


M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Jet Propulsion Engines 


Preliminary Performance Data of Sev- 
eral Tail - Pipe - Cascade - Type Model 
Thrust Reversers, by James G. Hensel, 


Jr., and Jack G. McCardle, NACA RM 
E55F09, Aug. 1955, 48 pp. 

Definitions of the Thrust of a Jet Engine 
and of the Internal Drag of a Ducted 
Body, by a special panel of the Perform- 
ance and Tagine Aerodynamics Sub- 
committees, Aeronautical Research Coun- 
cil, London, J. Roy. Aeron. Soc., vol. 59, 
Aug. 1955, pp. 517-526. . 

The Significance of Diffuser Losses for 
the Design of Supersonic Jet Engines, 
by O. Schrenk, Zeitschr. fir Flugwiss., 
vol. 3, June 1955, pp. 164-165 (in 
German). 

Effect of Three Flameholder Configura- 
tions on Subsonic Flight Performance of 
Rectangular Ram Jet Over Range of 
Altitudes, by Dugald O. Black and Wesley 
E. Messing, NACA RM ES8I01, Nov. 
1948, 28 pp. (declassified from Confiden- 
tial Sept. 15, 1955). 

Effect of Exhaust-Nozzle Ejectors on 
Turbojet Noise Generation, by Warren 
J. North and Willard D. Coles, NACA TN 
3573, Oct. 1955, 26 pp. 

Design Analysis of Allison 501, by David 
A. Anderton, Aviation Week, vol. 63, 
Nov. 7, 1955, pp. 32-33, 35, 37, 41, 43. 

Generalized Turbojet Performance 
Equations for IBM Calculation, by Victor 
P. Kovacik, Wright Air Dev. Center, Power 
Plant Lab., TN WCLP 55-39, March 3, 
1955, 23 pp. 

Pulsejet Helicopter Power Plant Control 
System Development, Final Report, by 
R. W. MecJones, Fairchild Eng. and Air- 
plane Corp., Amer. Helicopter Div., Rep. 
195-J-1, Dec. 15, 1954, 56 pp. 


Rocket Propulsion Engines 


ATO Rocket Burns Liquid Propellants, 
Aviation Age, vol. 24, Nov. 1955, pp. 
142-145. 

Rocket Engine Drives Missile Accesso- 
ries, by R. L. Zimmerman, Aviation Age, 
vol. 24, Nov. 1955, pp. 166-171. 

Effect of Heat Transfer on Rocket Mo- 
tor-Tube Strength, by Frank W. Jandl, 
Navord Rep. 3410; NOTS Rep. 1004, 
March 24, 1955, 18 pp. 

Research on Direct Measurement of 
Rocket Exhaust Velocities, Ohio State 
Univ. Res. Found., Quart. Prog. Rep. 
5, Sept. 1—-Nov. 30, 1955, 22 pp. 


Heat Transfer and 
Fluid Flow 


The Determination of Coefficients of 
Heat Transfer to a Rocket Motor Nozzle 
by a Transient Method, Part II, by W. 8. 
Long, Gt. Brit. Roy. Aircr. Estab. TN 
R.P.D. 114, Dec. 1954, 14 pp. 

a age of Stall in a Compressor 
Blade Row, by Frank E. Marble, J. Aeron. 
Sci., vol. 22, Aug. 1955, pp. 541-554. 

__ Simultaneous Measurements of Veloc- 
ity and Temperature Profiles for Adiabatic 
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Supersonic Flow of Air in a Tube, by 
Joseph Kaye, George A. Brown, Julius 
C. Wesiantaleed, and Edward F. Kurtz, 
Jr., J. Aeron. Sct., vol. 22, Aug. 1955, 
pp. 575-577. 

The Calculation of Diffusion in Centrif- 
ugal Turbocompressors, by G. L. Val- 
denazzi, Aerotecnica, vol. 35, April 1955, 
pp. 79-86 (in Italian). 

The Atomization of Liquids. Part IL 
by H. A. Troesch, Engrs. Digest, vol. 16, 
Aug. 1955, pp. 367-370. 

An Electro-Chemical Method of Visual- 
izing the Laminar Boundary Layer in 
Water, by Edward W. Geller, Inst. 
Aeron. Sci., Student Branch Competition, 
First Award Papers, 1954, pp. 28-32. 

Heat-Transfer Measurements in an 
Inexpensive Supersonic Wind Tunnel, 
he Seuehs Kaye, J. H. Keenan, G. A. 
Brown, and R. H. Shoulberg, J. Appl. 
Mech., vol. 22, Sept. 1955, pp. 289-310. 

The Theories of Turbulence, by L. 
Agostini and J. Bass, NACA 7M 1377, 
Oct. 1955, 163 pp. (translated from 
France, Min. de l’ Air, Pub. Sci. et Tech., 
no. 237, 1950). 

Production of High Speed Liquid Jets 
by Means of a Propellant, by Brian 
Dunne, Jr., Herbert Gass, and Benedict 
Cassen, Atomic Energy Comm., UCLA- 
275, Jan. 1954, 11 pp. 

Recent References to Thermodynamic 
Data, by Wendell M. Latimer, Atomic 
Energy Comm., UCRL-1680, Aug. 1954, 
27 pp. 

Thermal Conductivity of Gases at High 
Temperatures, by Albert Joel Rothman, 
Atomic Energy Comm., UCRL-2339, Jan. 
1954, 115 pp. 

Thermal Conductivity and Its Variabil- 
ity With Temperature and Pressure, by 
L. 8. Kowalezyk, Trans. ASME, vol. 77, 
Oct. 1955, pp. 1021-1035. 

Secondary Flow in Axial Turbomachin- 
ery, by L. H. Smith, Trans. ASME, 
vol. 77, Oct. 1955, pp. 1065-1076. 


Heat Transfer Problems in Rocket 
Motors, by H. Ziebland, /. Brit. Interplan. 
i vol. 14, Sept.-Oct. 1955, pp. 249- 


Heat Transfer Conditions at the Flame 
Tube Walls of an Aero Gas-Turbine Com- 
bustion Chamber, by E. F. Winter, Fuel, 
vol. 34, Oct. 1955, pp. 409-428. 


On the Stability of the Axially Sym- 
metric Laminar Jet, by H. G. Lew, Quart. 
a. Math., vol. 13, Oct. 1955, pp. 310- 

14. 


The Flow of Gases of Varying Specific 
Heat—Thermodynamic Properties of Air 
and Combustion Products of Hydrocarbon 
Fuels, Part II, by D. Fielding, J. E. C. 
Topps, and W. R. Thompson, Gt. Brit. 
Nat. Gas Turbine Esitab., Rep. R. 160, 
June 1954, 24 pp. 


Flow Phenomena in Starting a Hyper- 
sonic Shock Tunnel, by H. S. Glick, A. 
Hertzberg, and W. E. Smith, Cornell 
Aeron. Lab., Rep. Ad-789-A-3, March 
1955, 75 pp. 


The Effect of Stretching a Vortex Core, 
by John R. Weske, Maryland Univ. 
Inst. Fluid Mech. and Appl. Math., TN 
BN-57 (Off. Sci. Res. TN-55-282) Aug. 
1955, 15 pp. 


Flow Equations for Multicomponent 
Fluid Systems, Part I, General Equations, 
by Newman A. Hall, Minn. Univ. Mech, 
Engng. Dept. TR 2, Aug. 1955, 33 pp. 

Nozzle Flow With Heat Transfer and 
Friction, by W. Gordon Watson, -\irer. 
Engng., vol. 27, Oct. 1955, pp. 339-345. 

The Rate of Dissipation of Energy and 
the Energy Spectrum in a Low-S 
Turbulent Jet, by Wan-cheng Chiu and 
Louis N. Rib, New York Univ. Coll. 
Engng. Sci. Rep. 2, April 1955, 39 pp. 

Interference Due to Viscous Wakes 
Between Stationary and Rotating Blades 
in Turbomachines, by Rudolf X. Meyer, 
The Johns Hopkins Uni. Inst. for Coop, 
Res. Rep. I-24, May 1955, 102 pp. 


Combustion 


Investigations of Limits of Flame Sta- 
bility of Disk Flameholders in a Free 
Stream, by Gustave G. Kutzko, Battelle 
Mem. Inst. TR 15033-4, Jan. 1955, 15 pp. 

On Rendering Engine Exhaust Gases 
Nonexplosible, by B. P. Mullins and 
J. M. Hawkins, Gt. Brit. Nat. Gas Turbine 
Estab., Memo. Rep. M. 225, Oct. 1954, 
17 pp. 

Mechanism and Kinetics of the Reac- 
tion Between Fuming Nitric Acid and/or 
Its Decomposition Products and Gaseous 
Hydrocarbons, by Albert L. Myerson and 
Francis R. Taylor, The Franklin Inst., Labs. 
for Res. and Dev. Quart. Prog. Rep. Q- 
2452-2, Sept. 11, 1954-Jan. 15, 1955, 14 


p. 

Infrared Study of the Kinetics of Homo- 
geneous Gas Reactions, Terminal Report, 
Univ. of Calif., Berkeley, Oct. 15, 1954, 
21 pp. 

Flame Propagation Rates Along a Ni- 
chrome Surface, by Arthur L. Thomas and 
R. H. Wilhelm, Princeton Univ. Chem. 
Kinetics Project, TN 22 (Off. Sci. Res. 
TN-55-361), Sept. 1955, 36 pp. 

The Effect of Particle Size on the Com- 
bustion of Uniform Suspensions, by 
J. A. Browning, T. L. Tyler, and W. G. 
Krall, Project Squid, TR DART-6-P, 
Sept. 1955, 16 pp. (ASTIA AD 70,871; 
available only on microcard). 

Burning Velocities of Various Premixed 
Turbulent Propane Flames on Open Burn- 
ers, by Paul Wagner, NACA TN 3575, 
Oct. 1955, 32 pp. 

Study ot Screeching Combustion in a 
6-In. Simulated Afterburner, by Perry 
L. Blackshear, Warren D. Rayle, and 
Leonard K. Tower, NACA TN 3567, 
Oct. 1955, 55 pp. 

Integration of a Simplified Kinetic 
Model for a Hydrogen Bromide Flame, 
by Edwin 8. Campbell, Wisconsin Univ. 
Naval Res. Lab. CM-849, July 1955, 82 pp. 

Spark Ignition on Flowing Gases. V. 
Application of Fuel-Air-Ratio and Initial- 
Temperature Data to Ignition Theory, 
by Clyde C. Swett, Jr.. NACA RM 
E55116, Nov. 1955, 19 pp. 

Applications of High Speed Computing 
to Chemical Problems, by Joseph 0. 
Hirschfelder, Wisconsin Univ., Naval 
Res. Lab., Rep. WIS-ONR-16, Aug. 
1955, 12 pp. 

Studies of the Diffusion Coefficients of 
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pinary Gas Systems, by R. E. Bunde, 
Wisconsin Univ., Naval Res. Lab., CM- 
950, Aug. 1955, 149 pp. 

Interferometric Studies of Refractive 
Indices and Diffusion Coefficients for 
Selected Gases, by W. F. Rumpel, 
Wisconsin Univ., Naval Res. Lab., CM- 
851, Aug. 1955, 108 pp. 

The Thermal Decomposition of Tri-n- 
Butylborane, by Louis Rosenblum, J. 
Amer. Chem. Soc., vol. 77, Oct. 5, 1955, 
pp. 5016-5017. 

The Premixed Hydrogen-Fluorine Flame 
and Its Burning Velocity, by A. V. Grosse 
and A. D. Kirshenbaum, J. Amer. Chem. 
Sec., vol. 77, Oct. 5, 1955, pp. 5012-5013. 

The Influence of Methyl Bromide on 
Flames, Part 1, Premixed Flames, by 
R. F. Simmons and H. G. Wolfhard, 
Trans. Faraday Soc., vol. 51, Sept. 1955, 
pp. 1211-1217. 

The Behavior of Diborane-Air and of 
Propane-Diborane-Air Flames on a Vortex 
Burner, by Philip F. Kurz, Battelle Mem. 
Inst., 7 R 15036-24, March 4, 1955, 16 pp. 

Calculated Weak Limit Flame Tempera- 
tures of Hydrocarbon-Air-Diluent Mix- 
tures, by B. P. Mullins and J. M. Marley, 
Gt. Brit. Nat. Gas Turbine Estab., Mem. 
M. 235, March 1955, 13 pp. 

Effects of Additives on Flame Propaga- 
tion in Acetylene, III, by M.S.B. Munson 
and Robbin C. Anderson, Texas Univ., 
TN 17, Jan. 1955, 14 pp. 


Fuels, Propellants, and 
Materials 


Experience Teaches How to Form Ti- 
tanium, by J. M. Ohlson, SAE J., vol. 63, 
Aug. 1955, pp. 57-58. 

The Influence of the Geometric Forms 
of Powder Grains on Their Burning Rates, 
by George Seitz (translated from Spreng- 
technik no. 12, 1952, pp. 221-227, and 
Explosivstoff no. 1/2, 1953, pp. 15-18), 
Picatinny Arsenal, Trans!. 1, July 1955, 
29 pp. 

The Application of Subsonic Vibrations 
During Solidification of Castings with Par- 
ticular Reference to a Material for Gas 
Turbine Blades, by S. Hinchliff and Josiah 
W. Jones, Cranfield College of Aero. Rep. 
89, April 1955, 23 pp. 

Kinetics of the Thermal Decomposi- 
tion of Nitric Acid Vapor. III. Low 
Pressure Results, by Harold S. Johnston, 
Louise Foering, and James R. White, J. 
Amer. Chem. Soc., vol. 77, Aug. 1955, pp. 
4208-4212. 

Measurement and Correlation of Vapor 
Pressure Data for High Boiling Hydro- 
carbons, by H. S. Myers and M. R. 
Fenske, Indust. Engng. Chem., vol. 47, Aug. 
1955, pp. 1652-1658. 

Vapor Pressure Chart for All Types of 
Low Boiling Hydrocarbons, by H. S. 
Myers, Indust. Engng. Chem., vol. 47, Aug. 
1955, pp. 1659-1660. 

Properties of High Molecular Weight 
Hydrocarbons, by R. W. Schiessler and 
F. C. Whitmore, Indust. Engng. Chem., 
vol. 47, Aug. 1955, pp. 1660-1665. 

Heat Transfer Design Data—Aqueous 
Solutions of Nitric Acid and of Sulfuric 
Acid, by Thomas R. Bump and Wilmer 
L. Sibbitt, Indust. Engng. Chem., vol. 47, 
Aug. 1955, pp. 1665-1670. 

Phase Diagram of Liquid Ozone-Oxygen 
System, by Callaway Brown, Abraham 
W. Berger, and Charles K. Hersh, J. 
hey Phys., vol. 23, July 1955, pp. 1340- 
343. 

Structure of the Gas Phase Combustion 
Region of a Solid Double Base Propellant, 
by C. A. Heller and Alvin S. Gordon, 
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J. Phys. Chem., vol. 59, Aug. 1955, pp. 
773-774. 

The Thermal Decomposition of Nitrous 
Oxide. I. Secondary Catalytic and Surface 
Effects. II. Influence of Added Gases 
and a Theory of the Kinetic M i 
by F. J. Lindars and Cyril Hinshelwood, 
Proc. Roy. Soc., vol. (A)231, Aug. 22, 
1955, pp. 162-197. 

Transport Properties of Gaseous Hydro- 
carbons, by J. D. Lambert, K. J. Cotton, 
M. W. Pailthorpe, A. M. Robinson, J. 
Scrivins, W. R. F. Vale, and R. M. 
Young. Proc. Roy. Soc., vol. (A)231, 
Aug. 1955, pp. 230-290. 

Estimated Solubility of Helium in Red 
Fuming Acid, With Application to Gas- 
Pressurized Rocket Systems, by J. S. 
Gordon, Wright Air Dev. Center, TN 
WCLP 55-6, March 1955, 14 pp. 

The Liquid Ammonia Chemistry of 
Boron, by William L. Jolly, Atomic Energy 
Comm., UCRL-4332, April 1954, 17 pp. 

A Survey Report on Lithium Hydride, 
by Thomas R. P. Gibb, Jr., and C. E. 
Messer, Atomic Energy Comm., NYO- 
3957, May 1954, 36 pp. 


Instrumentation and Ex- 
perimental Techniques 


Research on Direct Measurement of 
Rocket Exhaust Velocities, Oliio State 
Univ. Res. Found. Quart. Prog. Rep. 5, 
Sept. 1-Nov. 30, 1954, 22 pp. 

A Method for the Analysis of Gas Mix- 
tures Containing Oxides of Carbon and 
Nitrogen, by R. H. J. Gercke and T. E. 
Hicks, Atomic Energy Comm., LRL-87, 
Feb. 1954, 6 pp. 

A Tensile Test for Rocket Propellants 
Extruded in the 2.75-In. Internal- 
Burning Grain Form, by F. J. Worcester 
and W. H. Gough, Naval Powder Factory, 
Mem. Rep. 97, April 7, 1955, 14 pp. 

Averaging of Periodic Pressure Pulsa- 
tions by a Total-Pressure Probe, by R. C. 
Johnson, NACA TN 3568, Oct. 1955, 


30 pp. 


Terrestrial Flight, Vehicle 
Design 


Stability of Spin-Stabilized Rockets 
With Very Large Initial Yaws, by B. A. 
Hodes and E. M. DePue, Aberdeen Prov- 
ing Ground, Ballistic Res. Lab., Rep. 926, 
Jan. 1955, 16 pp. 

Self-Aligning Rockets, by Serge J. 
Zarrodny, Aberdeen Proving Ground, Bal- 
listic Res. Lab., Mem. Rep. 884, Jan. 1955, 
21 pp. 

The Dynamic Motion of a Missile Des- 
cending Through the Atmosphere, by 
Hans R. Friedrich and Frank J. Dore, 
J. Aeron. Sci., vol. 22, Sept. 1955, pp. 
628-632, 638. 

Guided Missile Problems, by W. H. 
Pickering, Ordnance, vol. 40, Sept.-Oct. 
1955, pp. 241-243. 

Prediction of Missile Reliability, by 
Maurice J. Kirby and Harry R. Powell, 
Sperry Engng. Rev., vol. 8, July-Aug. 
1955, pp. 12-20. 

Optimum Flight Paths for Jet Aircraft, 
by Angelo Miele, NACA TM 1389, Sept. 
1955, 47 pp. (translated from Aero- 
tecnica, vol. 32, 1952, pp. 206-219). 


INSTRUMENTATIO. 


MODEL P146 
PRESSURE TRANSDUCER 
RANGES: 0-100 TO 
0-5,000 PSIG 


Model P146 Pressure Transducers 
employing the Statham biradial un- 
bonded strain gage were selected 
for use at this new facility because 
of their high accuracy, interchange- 
ability and capacity to perform dur- 
ing the grueling tests required in the 
development of rocket engines. 


Please request Bulletin No. P146. 


LABORATORIES 
LOS ANGELES 64, CALIFORNIA | 
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new Baldwin instrumentation system 
measures engine thrust electronically 


In the new test cells at an Air Force base in California, jet engine 
thrust and consumption of lubricating oil are measured by a Baldwin 
instrumentation system that is instantaneously responsive and accu- 
rate to 

At the heart of the system are the Baldwin SR-4® transducers 
which sense changes in load and convert them into electrical quan- 
tities. Load cells are hermetically sealed and have no moving parts, 
so the system is always reliable and virtually maintenance-free. Signals 
are transmitted from the SR-4 load cells, amplified, and indicated 
on specially modified 5-range potentiometers which are calibrated 
in pounds for direct reading. Every Baldwin system is pre-loaded and 
calibrated to achieve the guaranteed accuracy. Signals may also be 
fed into a digital computer. 

Baldwin designs and builds instrumentation systems for any appli- 
cation involving force, weight, pressure, torque, flow, tension—from 
simple measuring and indicating systems to complete feedback 
control. Next time you have a control problem, why not let our 
engineers help you? In the meantime, for a copy of our illustrated 
bulletin 4309, write to Electronics & Instrumentation Division, 
Baldwin-Lima-Hamilton Corporation, 806 Massachusetts Avenue, 
Cambridge, Mass. 


£. > BALDWIN-LIMA-HAMILTON 


Divisions: Austin-Western * Eddystone * 
Electronics & Instrumentation * Hamilton * 
Lima * Loewy-Hydropress * Madsen ° Pelton 
* Standard Steel Works 


Space Flight, Astro- 
physics, Aerophysics 


An Electrical Machine for Use in Extra. 
Terrestrial Environment, by G. S. Bro. 
san, J. Brit. Interplan. Soc., vol. 14, 
Sept.-Oct. 1955, pp. 270-273. 

instrumentation of the Rocket-Grenade 
Experiment for Measuring Atmospheric 
Temperatures and Winds, by W. G. 
Stroud, E. A. Terhune, J. H. Venner, 
J. R. Walsh, and 8S. Weiland, Rev. Sci. 
Instrum., vol. 26, May 1955, pp. 427-432. 

Projecting the Earth, by Otto Weil, 
Sperry Engng. Rev., vol. 8, May-June 
1955, pp. 16-20. 

Dynamic Problems of Interplanetary 
Flight, by Derek F. Lawden, Aeron. 
Quart., vol. 6, Aug. 1955, pp. 165-18). 


Atomic Energy 


Challenge to the Aircraft Propulsion 
Engineer: The Nuclear Reactor, by Wil- 
liam J. Graff, Jr., Aviation Age, vol. 24, 
Oct. 1955, pp. 18-25. 

Reactor Fuels, by J. P. Howe, Chem. 
Engng. Progress, vol. 51, Oct. 1955, pp 
482-484. 

The Safety of Nuclear Reactors, by C. 
Rogers McCullough, Mark M. Miills, 
and Edward Teller, Chem. Engng. Prog- 
ress, vol. 51, Oct. 1955, pp. 446-450. 

Critical Assemblies at Los Alamos, by 
H. C. Paxton, Nucleonics, vol. 13, Oct. 
1955, pp. 48-50. 


NORTHAM 


Miniature Magnetic 
Tape Recorder 


Height: 5’’—Diameter: 4’ 


Weight: 4% Ibs. 


Eight channels of information can be applied 
simultaneously to a %-inch tape in the 
Model MR-1 Recorder (illustrated). This 
self-contained unit includes transistor tim- 
ing oscillator and battery power source and 
operates under axial accelerations up to 
500 g. A pre-record carrier system permits 
direct coupling with D.C. devices such as 
thermocouples and strain gauge transducers. 


MODEL MR-1 SPECIFICATIONS: 


Sensitivity :....28 miHivolts D.C. full scale 
across 40 ohms impedance 

Frequency Bandwidth:.......... 0-250 cps 

Timing: ...... 1000 cps transistor oscillator 

Applications: ....Running time, ber of 


channels and recording characteristics are 
adaptable to meet requirements of appli- 
cations—including rocket and aircraft 
flight testing, oil well logging and 
atmospheric studies. 


WRITE FOR BULLETIN NO. ...JP-104 
NORTHAM PRODUCTS INCLUDE... 


Transd s for p " leration and 
disp! t ent and auxiliary 
electronic equip t for complete systems. 


NORTH AMERICAN INSTRUMENTS, INC. 


2420 North Lake Avenue * Altadena, Calif 
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HIGHER 


SPEEDS 


MODELS 65F005 
PAD AND 10260* 


Rated Capacity: 0.5 gpm 2 1500 rpm. 
Max. Continuous Speed: 10,000 rpm. 
Continuous Pressures: 1500 psi and 3000 psi. 
Weight: 31 oz. 

MODELS 65F010 


PAD PER AND 10260°* 


Rated Capacity: 1.0 gpm £3 1500 Sg 

Max. Continuous Speed: 1 

Continuous psi and 3000 psi. 
Weight: 3.2 I 


MODELS 65F015 
PAD PER AND 10260" 


5 44 


Rated Capacity: 1.5 gpm @ 1500 a? 

Max. Continuous Speed: 10,000 rp 

Continuous Pressures: 1500 psi psi. 
Weight: 4.0 Ib. 


MODELS 65F020 
PER AND 10261" 


Rated Capacity: 2.0 gpm 
Max. Continuous Speed: 7500 
Continuous Pressures: 1500 one 3000 psi. 
Weight: 5.5 Ib. 


MODELS 65F030 
PAD PER AND 1026!* 


Rated Capacity: 3.0 gpm @ 1500 rpm. 

Max. Continuous Speed: 7500 rpm. 
Continuous Pressures: 1500 psi and 3000 psi. 
Weight: 6.0 Ib. 


*Other mounting flanges and port arrangements available. 


WATERTOWN 


SMALLER 
ENVELOPES 


A New Series of Fixed Displacement Pumps 
That Check Out on the Points You Want Most! 


LIGHTER HIGHER 


WEIGHTS 


CEILINGS 


STRATOPOWER 


HYDRAULIC PUMPS 


“\ The 65F Series STRATOPOWER Hydraulic 


Pumps pack more power than ever into ex- 
tremely small envelopes . . . check the draw- 
ings! These Fixed Displacement Pumps operate 
at higher speeds .. . check the figures! 


The 65F Series Pumps feature a piston actua- 
tion mechanism which is unaffected by fluid tem- 


peratures or inlet pressures. They will operate . 


with reservoirs pressurized to over 80 psi abso- 
lute, yet do not require pressurization for altitude 


operation. 
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THE NEW YORK AIR BRAKE (ff) 


STARBUCK AVENUE e 
INTERNATIONAL SALES OFFICE, 90 WEST ST, NEW YORK 6, N. Y. 


WATERTOWN-:WN Y. 
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Data 
Screw Torque Data 
Adapter Problems 

General 


frying to male 
Faced with “‘impossible’’ problems 
of design diminution, 
weight and size reduction, 
power and component capsulation? x 
Send for complete story on how MPB's 


BALL BEARINGS ACTUAL SIZE 


such as these are helping to make things 


....$maller and smaller and smaiier 


MINIATURE PRECISION BEARINGS, INC. 


20 Precision Park, Keene, N. H. 


ADEL has the answer to the need for 
a rugged, heat and abrasion-resistant 
clamp that can withstand the attack of 
the stratosphere’s sub-zero cold and 
the intense heat of today’s aircraft, mis- 
sile and rocket power plants. TEMPERA- 
TURE CHARACTERISTICS: Far above 
+550° F. to well below —90° F. Per- 
formance far beyond specifications of 
any other clamp. 


SERVICE-FITTED 
SERVICE-TESTED 
SERVICE-APPROVED 


~ SILICONE IMPREGNATED 


FIBERGLAS CUSHIONED 


CLAMPS 


Withstand Jet Engine Heat or Stratospheric Cold... 


Easily installed on “bundle” or removed, 
and they retain their firm, cushioned 
grip despite severe vibration and extreme 
temperatures. Cut maintenance and 
replacement costs...cushions won’t tear 
during engine ‘harness assemblies and 
if material chars, due to extremely high 
temperatures, there still remains a 
safe, non-conductive layer—A positive 
“design-fix” for jet engines. 


SAMPLES & SPECIFICATIONS: Aircraft, Component & Guided Missile Manufacturers mail 
coupon today for samples and detailed information on how CLAMPS fit your specific needs. 
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speeds supersonic ramjet progress 


Putting the modern ‘‘weapons system”’ into 
action...that’s the purpose of the OMAR 
Joint Technical Committee. 


Here undivided responsibility in the devel- 
opment of supersonic aircraft and propelled 
missiles is a reality—from theory and draw- 
ing board, through design and engineering, 
to testing and full-scale production. This 
team activity, coordinating for the first time 
the chemical and mechanical experience of 
Marquardt Aircraft, Olin Mathieson Chem- 
ical and Reaction Motors, promises: con- 
stant advancement in the development and 
production of ramjets, rockets, and special 
liquid and solid propellants of a practical 
nature for supersonic flight. 


Laboratory testing of a Marquardt supersonic ramjet engine 
capable of producing 10,000-plus horsepower at high altitudes. 
Here testing facilities can simulate altitudes to 100,000 feet and 
air flows equivalent to three times the speed of sound. Marquardt 
today ranks as the West's largest jet engine research and 
development center for ramjet power plants. Testing is just one 
facet in the coordinated activities of Marquardt Aircraft Com- 
pany, Olin Mathieson Chemical Corporation and Reaction Motors, 
Inc. dedicated to the improvement of ramjet and rocket engines 
and special fuels for supersonic propulsion. 


RMI liquid propellant rocket engines 
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Marquardt Aircraft Company Olin Mathi Chemical Corporation Reaction Motors, Inc. 
Van Nuys, California New York, New York Denville, New Jersey 
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THE M. W. KELLOGG COMPANY, NEW YORK 7, N.Y. 9 
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SUBSIDIARIES OF PULLMAN INCORPORATED CHEMICAL 


PLANTS 


136 JET PROPULSION 


| 
Outline o emperature-Pressure Problems 
‘Any rocket engine looks) le electric power plants, the Kellogg assign: 
ly problems include Selection'or development 
= ne M. W. ogg Con roper ‘to provide correct 
n. to corrosion or erosion, but also the welding 
tech S, ferials, and skillsto fabricate 
the c opment 2 d pre thea an ‘se or intricate 
: ‘or example, have much: concerm dy rocket rogram, andto 
| latter, nozzle temperature those responsible for process equipment in 
4 
CHIMNEYS 


